M Khaska

Négrel Ph

Tracking selenium in the Chalk aquifer of northern France: Sr isotope constraints

Groundwater at the southern and eastern edges of France's Paris Basin has a selenium content that at times exceeds the European Framework Directive's drinking-water limit value of 10 µg/L. To better understand the dynamics of the Chalk groundwater being tapped to supply the city of Lille and the Se origins, we used a combination of geochemical and isotopic tools. Strontium isotopes, coupled with Ca/Sr, Mg/Sr and Se/Sr ratios, were used to identify the main groundwater bodies and their mixings, with the Mg/Sr ratio constraining a ternary system. Groundwater in the agricultural aquifer-recharge zone represents a first endmember and displays the youngest water ages of the catchment along with the highest strontium isotopic signature (0.7084) and low Se contents. Anaerobic groundwater constitutes a second major end-member affected by water-rock interactions over a long residence time, with the lowest strontium isotopic signature (0.7079) and the lowest Se content, its low SF6 content confirming the contribution of old water. Se-rich groundwater containing up to 30 µg/L of Se represents a third major end-member, with an intermediate Sr isotopic ratio, and is mainly constrained by the clayey Se-rich formation overlying the Chalk aquifer. The spatial and temporal Se variability in the groundwater is clearly linked to the presence of this formation identified as Tertiary and also to the hydrological conditions; saturation of the Se-rich clays by oxygenated groundwater enhances Se mobility and also Sr adsorption onto the clays. This multi-tool study including Sr isotopes successfully identified the Se origins in the aquifer and has led to a better understanding of the regional mixing and processes affecting the Chalk groundwater.

1.

Introduction

High concentrations of selenium (Se) in both surface and groundwaters over the past 30 years have emerged as a major issue world-wide. Se, which occurs naturally as a trace constituent in geologic formations (e.g. Cretaceous black shale in the USA [START_REF] Kulp | Speciation and weathering of selenium in upper Cretaceous chalk and shale from South Dakota and Wyoming, USA[END_REF]Tuttle et al., 2014a) and associated soils, is currently contaminating watersheds in mining areas [START_REF] Chuncai | Mobility behavior and environmental implications of trace elements associated with coal gangue: A case study at the Huainan Coalfield in China[END_REF][START_REF] Liu | Selenium Speciation in coal ash spilled at the Tennessee Valley Authority Kingston Site[END_REF][START_REF] Vengosh | Isotopic imprints of mountaintop mining contaminants[END_REF], irrigated alluvial valleys [START_REF] Bailey | Simulating variably-saturated reactive transport of selenium and nitrogen in agricultural groundwater systems[END_REF][START_REF] Bailey | The Influence of nitrate on selenium in irrigated agricultural groundwater systems[END_REF][START_REF] Bajaj | Hazardous concentrations of selenium in soil and groundwater in North-West India[END_REF][START_REF] Deverel | Relation of salinity and selenium in shallow groundwater to hydrologic and geochemical processes, Western San Joaquin Valley, California[END_REF][START_REF] Dhillon | Quality of underground water and its contribution towards selenium enrichment of the soil-plant system for seleniferous region of Northwest India[END_REF][START_REF] Gates | Assessing selenium contamination in the irrigated stream-aquifer system of the Arkansas River, Colorado[END_REF]Tuttle et al., 2014b), wetlands [START_REF] Zawislanski | Selenium distribution and fluxes in intertidal wetlands, San Francisco Bay, California[END_REF][START_REF] Zhang | Isotopic and microbiological signatures of pyrite-driven denitrification in a sandy aquifer[END_REF], and sediments [START_REF] Beisner | Selenium and trace element mobility affected by periodic displacement of stratification in the Great Salt Lake, Utah[END_REF][START_REF] Guan | Selenium distribution in the Rhône delta and the Gulf of Lions[END_REF][START_REF] Masscheleyn | Transformations of selenium as affected by sediment oxidation-reduction potential and pH[END_REF][START_REF] Oremland | Selenate reduction to elemental selenium by anaerobic bacteria in sediments and culture: biogeochemical significance of a novel, sulfate-independent respiration[END_REF]…

Se is a non-metallic element of both nutritional and toxicological importance, but the window between its beneficial contribution to health and its toxic contribution (from excess to deficiency) is very narrow and the safety margin is unknown [START_REF] Vinceti | Risk of chronic low-dose selenium overexposure in humans: insights from epidemiology and biochemistry[END_REF]. [START_REF] Barron | Selenium exposure in subjects living in areas with high selenium concentrated drinking water: Results of a French integrated exposure assessment survey[END_REF] showed that the main human Se inputs come from drinking water in regions where the groundwater has a high Se level. For the purpose of health control, the Water Framework Directive (2000/60/EC) defines an Se threshold value of 10 µg/L for human consumption.

Nevertheless, there is no consensus on what is an appropriate standard of Se intake and [START_REF] Vinceti | The need for a reassessment of the safe upper limit of selenium in drinking water[END_REF] have reaffirmed the need for a reassessment of the safe upper limit of selenium in drinking water.

Analyses of some groundwater bodies at the southern and eastern edges of France's Paris Basin over the past decade have shown a high Se content (Fig. 1) which raises concern as regards groundwater and tap-water quality. The high Se content in the Paris Basin is related mainly to the presence of Eocene, mainly Ypresian, and Oligocene formations [START_REF] Chabart | Origine des anomalies en sélénium dans les captages d'AEP du département de la Marne. Première approche[END_REF]Gourcy et al., 2011). The Ypresian formations were deposited according to various marine transgressions and regressions, and are made of limestones, clays, carbonated sands and fluvial deposits rich in organic matter. The Oligocene formations are essentially made of sands. In the Lille area (northeastern France), a high Se spatial variability has been observed within the well fields since 2000; in 2012, the Se range was between 0.04 to 30 µg/L. This has resulted in water operators no longer exploiting the most Se-contaminated wells.

The aim of the present study is to understand the interactions between different water bodies and to better assess the origins of the Se, given the constraint of the high pumping rates in the four studied well fields located in an agricultural, urban, and former industrial region southwest of Lille (Fig. 2). The temporal and spatial heterogeneities of the groundwater chemistry of four well fields have been documented along with a regular monthly monitoring.

The results were examined in terms of major-element, Sr and Se concentrations, H, O and Sr isotopes, and chlorofluorocarbon (CFC) and sulphur hexafluoride (SF6) measurements.

Our first objective was to understand the spatial disparities in the Sr isotopic signatures and Sr and Se contents, both in the regional Chalk aquifer and also more locally at well scale.

Our second objective was to investigate anthropogenic and lithological influences on the groundwater chemistry. Groundwater sources of Se from recharge sources and water pathways within the basin were traced using Sr isotopic signatures. This approach leads to a better understanding of the global functioning of the Chalk aquifer and confirms that Sr isotopes can be valuably used for tracing the origins and mobilization of Se in the Chalk aquifer.

Geological and hydrogeological context

In the study area, the stratigraphic succession of the Paris Basin begins with a highly faulted, more-or-less dolomitic, several hundred meter thick sequence of Lower Carboniferous limestone and calcareous shale. This is overlain by marly Cenomanian chalk, locally called the "Dièves", resulting from the Cretaceous transgression (Fig. 2). The "Dièves" is followed by a thick deposit of Turonian marly limestone that evolves upward to white and grey flintbearing Turonian and Senonian Chalk with phosphate beds. The unconformably overlying Tertiary formations comprise the Louvil Clay (Lower Landenian, 9-12 m thick), the Ostricourt Sand (Upper Landenian, 25-30 m thick), and the Flandres Clay (Ypresian, 10-15 m thick). It is the Tertiary formations throughout the Paris-London Basin that confine the underlying Chalk aquifer, inducing a redox boundary [START_REF] Edmunds | Residence time indicators in groundwater: the East Midlands Triassic sandstone aquifer[END_REF][START_REF] Elliot | Hydrochemical trends, palaeorecharge and groundwater ages in the fissured Chalk aquifer of the London and Berkshire Basins, UK[END_REF][START_REF] Kloppmann | Redox zones in the chalk aquifer of the Paris and North German basins[END_REF][START_REF] Smedley | Redox patterns and trace-element behaviour in the East Midlands Triassic Sandstone Aquifer[END_REF]. Quaternary loess (pre-Weichselian and Weichselian) covers practically the entire basin.

Given the structural context, a new hypothesis is suggested for the age of the clays present at 5-15 m depth in the Wells A11 to A15 area (Fig. 2), newly interpreted as Early Tertiary (Thanetian) instead of Quaternary as described in the geological database. These clays were protected from erosion through fault activation (Lutetian, Bartonian and Oligocene) prior to the Quaternary. The last 5 m are constituted by the quaternary alluvium of the Deûle River, mainly comprising heterogeneous sandy clayey deposits.

The Dieves marl forms the impermeable base of the Chalk aquifer (Fig. 2). The Turonian-Senonian Chalk constitutes the Chalk aquifer encountered at 10-50 m depth, favored by the presence of highly permeable horizons with open fissures mainly found at 10-15 m depth [START_REF] Caous | Impact de l'activité industrielle sur la qualité de la nappe de la craie (champs captants du Sud de Lille -59). Etat initial des connaissances géologiques et hydrogéologiques[END_REF]. The fissures and stratification joints in the Chalk close with depth, resulting in a significant decrease in rock permeability. The inter-annual average rainfall is 689 mm with an effective rainfall of 177 mm, and the main recharge area is close to Emmerin [START_REF] Caous | Impact de l'activité industrielle sur la qualité de la nappe de la craie (champs captants du Sud de Lille -59). Etat initial des connaissances géologiques et hydrogéologiques[END_REF]. The groundwater of the Chalk aquifer circulates in the fault network, which is highly developed in the valleys where the well fields are sited.

Drinking water for the urban community of Lille (more than 1 million inhabitants) is pumped through four well fields tapping the Chalk aquifer to the southwest of Lille (Fig. 2 from the Les Ansereuilles well fields has drained this source [START_REF] Talbot | Champ captants de Ansereuilles, d'Emmerin et d'Houplin-Ancoisne. Modélisation des écoulements souterrains et des transferts de masse. Acquisition de données[END_REF], Bourg et al., 1985).

The Emmerin and Houplin-Ancoisne well fields will be referenced hereafter as E-H.A. The Emmerin well field is constituted by 8 wells located at the top of a WNW-ESE anticline, the thickness of the Cretaceous beds is reduced as shown in the cross section (Fig. 2). The Houplin-Ancoisne well field has 10 wells (H1 to H10) aligned along the Seclin canal, which follows the trajectory of the major NW-SE-trending Seclin Fault, and one offset well (H11) (Fig. 2). The Seclin Fault affects the depth of the Chalk (9 m in Well H4; 6 m in Well H11).

The Chalk groundwater is locally confined below the clayey deposits (present in Well H11).

The Les Ansereuilles well field has 27 wells and is parallel to the highly polluted Deûle River.

The groundwater, which is unconfined to the southwest of the Deûle, flows northwestward to become confined under the alluvial cover whose thickness locally reaches 15 to 20 m. It is also locally confined in the A11-A15 well area due to the presence of more than 10 m of clay deposited above the Chalk (Fig. 2). In the Salomé well field (17 wells), the Chalk is covered with a layer of sandy clayey silt, 1 to 10 m thick, making the aquifer unconfined, with the groundwater circulating SW to NE with a gradient of 1.5‰. Chemical separation of the Sr was done using an ion-exchange column (Sr-Spec), with a total blank of 0.5 ng for the entire chemical procedure. After chemical separation, 1/5 of the sample was loaded onto a tungsten filament and analysed with a Finnigan MAT 262 multiple collector mass spectrometer. The 87 Sr/ 86 Sr ratios were normalised to an 86 Sr/ 88 Sr ratio of 0.1194. The reproducibility of the 87 Sr/ 86 Sr ratio measurements was tested through duplicate analyses of the NBS 987 standard and the mean value was 0.710243 ± 12.10 -5 (2σ, N=18).

Material and methods

A

The 2 H and 18 O isotopes of the water molecule were measured on raw water using a

Finnigan MAT 252 mass spectrometer with a precision of 0.1‰ for δ 18 O and of 0.8‰ for δ 2 H.

The isotopic compositions are plotted on the usual -scale in ‰ with reference to V-SMOW according to the equation δ sample (‰)=(R sample /R standard )*1000, where R is the 2 H/ 1 H and 18 O/ 16 O atomic ratios.

Chlorofluorocarbons (CFC-11, CFC-12, CFC-13) and sulphur hexafluoride (SF6) were analysed on raw water sampled in 500 mL glass bottles using gas chromatography in the Spurenstofflabor laboratory (Germany). The detection limit was close to 10 -4 pmol, allowing the measurement of CFC concentrations down to 0.01 pmol/L and SF6 concentrations to 0.1 fmol/L. The reproducibility was about ±5% for water samples. The calculations follow the programs developed by the International Atomic Energy Agency (IAEA, 2006). The CFC-11 and CFC-12 were used mainly for refrigeration and air-conditioning, while CFC-113 was used as a solvent and SF6 for its electrical and thermal insulation properties. CFC and SF6

tools are increasingly being used to date recent groundwater components allowed the apparent groundwater ages to be estimated [START_REF] Darling | The practicalities of using CFCs and SF6 for groundwater dating and tracing[END_REF][START_REF] Gooddy | Using chlorofluorocarbons (CFCs) and sulphur hexafluoride (SF6) to characterise groundwater movement and residence time in a lowland Chalk catchment[END_REF]Gourcy and Brenot, 2011).

For the solid samples, flint and white Chalk were sampled from the Emmerin quarry, and the Tertiary deposits (Landenian Ostricourt Sand and Ypresian Flandres Clay) and Quaternary silt were sampled from the Templeuve quarry. Cuttings of sediments were done by drilling with a helical auger without fluid addition close to wells A13, H4, H11 and S3 till reaching the chalk. After total rock dissolution by classical tri-acid mineralization, the chemical composition of the major and trace elements was measured by ICP-AES. The strontium isotope composition of solid samples (except sediments) was obtained after simple leaching with deionized water and analysed with a Finnigan MAT 262 multiple collector mass spectrometer [START_REF] Négrel | Strontium isotopes as tracers of groundwater-induced floods: the Somme case study (France)[END_REF].

Results

The concentrations of the major elements, Se and Sr, δ 2 H and δ 18 O, and 87 Sr/ 86 Sr are given under Supporting Information in Table 1. The compositions of the rocks are presented under Supporting information in Table 2.

Groundwater chemistry

The groundwater in the E-H.A area showed the lowest Na, Mg, and HCO 3 concentrations, with Sr close to the average concentration in the Deûle River (Supporting Information Table 1). Conversely, the groundwater in Well A26 and its neighboring wells A25, A27 and A28 in the Les Ansereuilles area, stood out in having the highest concentrations in Cl (1.9 to 3 mmol/L), Na, Mg, HCO 3 and Sr, and the lowest Ca content similar to that measured in the Deûle River (around 3 mmol/L). The Sr concentration in Well A26 was 10 times higher than

CFCs and SF6 contamination

All the groundwater samples were contaminated by CFCs and SF6 at various levels. All the groundwater samples were contaminated by CFCs and SF6 at various levels. The groundwater in wells A3, A4, H4, H6, and S3 was found to have higher CFC-12 concentrations than modern air-equilibrated water (2.97 pmol/L of CFC-12) [START_REF] Gooddy | Using chlorofluorocarbons (CFCs) and sulphur hexafluoride (SF6) to characterise groundwater movement and residence time in a lowland Chalk catchment[END_REF], making CFC-12 ineffective as a groundwater dating tool. Moreover, CFCs may possibly be reduced within the confined parts of aquifers. Therefore, only SF6 data were used for dating.

The initial program showed that SF6 concentrations in the groundwater varied between 0.7±0.1 fmol/L (Well A14) and 3.6±0.4 fmol/L (Well S8), the latter being above modern values (2.3 fmol/L, [START_REF] Gooddy | Using chlorofluorocarbons (CFCs) and sulphur hexafluoride (SF6) to characterise groundwater movement and residence time in a lowland Chalk catchment[END_REF] (Supporting information, Table 1). The group most contaminated by SF6 comprises wells E10, H4, H6 and H9, with concentrations equaling that of modern air-saturated water. Samples showing a higher SF6 concentration than air suggest contamination; Well A23 was probably contaminated by the nearby treatment plant. The wells located in the depressed area in the northeastern part of the Les Ansereuilles catchment displayed a relatively low contamination (0.7 and 0.9±0.1 fmol/L respectively for A14 and A12. The lowest concentrations were found in the confined Well A26 with 0.2±0.1 fmol/L for SF6.

SF6 concentrations measured during the monitoring program were high in Well A24 (from 1.5 to 2.6 fmol/kg), and even higher in Well H4 (from 1.7 to 3.3 fmol/kg, mean of 2.1 fmol/kg).

Lower, but varying, SF6 concentrations were measured in Well A13 (from 0.9 to 1.7 fmol/kg, with a mean of 1.2 fmol/kg). Well S3, however, displayed very high SF6 concentrations varying between 1.8 and 4.7 fmol/kg, indicating possible contamination.

Discussion

The spatial disparities observed in the element concentrations and 87 Sr/ 86 Sr values for the wells could reflect varied water pathways, varied recharge of the groundwater bodies, and varied land uses in the catchment.

Aquifer recharge, surface water / groundwater interactions, and spatial variability

The stable isotope variations in the groundwater are due mainly to (i) natural variations in the isotopic composition of the rainfall, (ii) the mixing with pre-existing waters, and (iii) the influence of evaporation during percolation through the soil or the unsaturated zone. CFCs from the atmosphere dissolve in water according to their Henry's Law solubilities at the recharge temperature [START_REF] Plummer | Chlorofluorocarbons[END_REF] which, in the study area, can be estimated at 12 °C based on the regional mean annual air temperature and can be used as age tracers. Knowing that the average altitude of the sampling points is around 20 m.a.s.l., the error in the estimation of apparent ages due to the uncertainty of point elevation is very low. The uncertainty due to an error in estimating the recharge temperature is 5% using the binary mixing model.

The δ 18 O and δ 2 H values of the groundwater samples collected in the four well fields lie within their 95% uncertainty envelope between the global meteoric water line (GMWL) and the meteoric line at Trier (Germany), , which is the closest locality with sufficient data from the Global Network of Isotopes in Precipitation (GNIP) [START_REF] Craig | Isotopic variations in meteoric waters[END_REF] (Fig. 3). It supports the idea that the groundwater most probably comes from present-day precipitation; the groundwater and surface-water samples showed only a narrow range of δ 18 O and δ 2 H signatures, especially when compared to the rain input throughout France in which δ 18 O may vary from -16‰ to 0‰ and δ 2 H from -120‰ to 0‰ [START_REF] Millot | Multi-isotopic composition (δ7Liδ11B-δD-δ18O) of rainwaters in France: Origin and spatio-temporal characterization[END_REF]. This narrow range reflects groundwater mixing over at least one hydrologic cycle with buffering of the different water inputs [START_REF] Huddard | dD and d18O evidence for inputs to groundwater at a wetland coastal boundary in the southern Great Lakes region of Canada[END_REF].

The Emmerin groundwater samples plot at the intersection between the local Trier meteoric line and a 4.5 slope line, and displayed the most depleted values of the unconfined aquifer.

The Les Ansereuilles groundwater signatures showed a spatial isotopic variability. Except for Well A26, this variability corresponds exactly to the range of the summer and winter isotopic means of the Deûle River water. Knowing that evaporation lines have a slope between 3.5 and 6.0 as function of local atmospheric humidity [START_REF] Gat | Isotopic fractionation[END_REF], this means that both the surface water and groundwater include slightly evaporated water. All tributaries to the Deûle River receive surface runoff from agricultural and industrial-commercial areas. The river water samples display a mixing with both rain water and stagnant water that has evaporated within the remaining wetlands where water renewal is limited; this is clear in the water from the Seclin canal which displays the highest enrichment. The fractionation trend in the groundwater can also be explained by partial evaporation of the water before infiltration into the aquifer, as has already been showed for the Chalk aquifer in Champagne [START_REF] Vachier | Mouvements de l'eau dans la zone non saturée et alimentation de la nappe de la craie de Champagne (France)[END_REF]. Interestingly, data acquired in the 1980's from the same Chalk aquifer in the Lille region were close to the GMWL and did not display an evaporative contribution [START_REF] Bernard | Acquisition et rassemblement des données géothermiques disponibles et nouvelles dans la zone Franco-Belge de Saint-Ghislain à Saint-Amand-les-Eaux (Nord)[END_REF][START_REF] Garnier | Un modèle de l'estimation de l'activité initiale du carbone dissous et des vitesses d'écoulement en régime captif[END_REF][START_REF] Simon | Apports de la biogéochimie isotopique à la connaissance de la dénitrification dans les eaux souterraines[END_REF]. The appearance of this rather recent evaporative contribution may be linked to a change in the recharge contribution, i.e. the intense pumping rates of the last 30 years may have triggered either (i) leaching of the entire unsaturated profile, thus soliciting a higher contribution of recharge water, or (ii) a contribution of Deûle River water, both integrating an evaporative component. The groundwater of Well A26 was the most depleted in δ 18 O (-7.2‰) and plotted on the local meteoric line. It can be considered as having the lowest influence from surface and/or alluvial water.

The absence of variation in the stable isotopic composition of the water over time from the monitored wells suggests that the groundwater can buffer the seasonal variations in rainwater input ( [START_REF] Gooddy | Using chlorofluorocarbons (CFCs) and sulphur hexafluoride (SF6) to characterise groundwater movement and residence time in a lowland Chalk catchment[END_REF]Gourcy and Brenot, 2011). In some wells (S8 and H9), SF6 contamination made the apparent age estimates doubtful (Supporting Information Table 1). At site scale, the apparently youngest groundwater (mean modern fraction of 91% for Emmerin and 99% for wells H4, H6 and H9) thus reflect the local role of recharge in the area which can be highly contaminated by present or past activities. The modern fraction of Well H11 (72%) was significantly lower, which may be explained by the regional Seclin Fault, used to advantage by water operators, possibly giving rise to mixing with pre-tracer water flowing along the fault.

In the Les Ansereuilles area, the modern water fraction for Well A26 is only 10%, thus confirming a large contribution of old water. SF6 data for the groundwater of Well A23, close to a sewage plant, displayed a very recent contribution of recharge water (99%). Surprisingly, two neighboring wells (A17 and A18) displayed very different modern fractions of 54% and 92%, respectively. The presence of a fault having a screen role in the Chalk between wells A17 and A22 may induce mixing with old water in Well A17. The wells in the A11-A15 area of Les Ansereuilles displayed a smaller modern fraction (35% for A12 and 65% for A14).

Although the Deûle River water was not analysed, we can assume that its SF6 concentrations are close to or above modern air-saturated water [START_REF] Gooddy | Using chlorofluorocarbons (CFCs) and sulphur hexafluoride (SF6) to characterise groundwater movement and residence time in a lowland Chalk catchment[END_REF].

With the monitoring program, the modern fraction of Well A24 groundwater was almost constant (between 66% and 74%, average of 70%); this would be equivalent to approximately 3 parts modern water to 1 part pre-tracer water when considering a binary mixing model. For Well A13, the modern fraction was higher during the winter months, varying from 34% in May 2012 to 65% in February 2012 (average of 46%), apart from the July measurement which gave 53%. A seasonal variation was also noted in Well H4, with a higher modern fraction in winter, i.e. 65% in June 2012 to 9-99% in the 2012-2013 winter.

Well S3 displayed the highest contribution of modern water (from 75% in May 2012 to 100% in the 2012-2013 winter). Knowing that the Tertiary and Quaternary sediments (T-Q) can be rather thin in Salomé, fast transfers could dominate.

Groundwater types and origins

The chemical facies of the groundwater samples were of a relatively homogeneous Ca-Mg-HCO 3 type. Carbonate dissolution in the Chalk aquifer leads to the release of Ca, Mg and Sr into the water, justifying the use of the bivalent Ca/Sr and Mg/Sr ratios as indicators of waterrock interaction; these are directly related to the groundwater residence time [START_REF] Edmunds | Baseline geochemical conditions in the Chalk aquifer, Berkshire, UK: a basis for groundwater quality management[END_REF]. The chemistry of both the groundwater and the solids are shown in the following figures (4 to 8).

The groundwater type evolves between two facies: (1) the E-H.A type (excluding Well H11) with low Na, Mg, and Sr concentrations indicating a relatively short residence time with little water-rock interaction, and (2) the confined type, represented by the groundwater of Well A26 with its high Na, Mg, and Sr concentrations indicating a long residence time with greater water-rock interaction [START_REF] Edmunds | Baseline geochemical conditions in the Chalk aquifer, Berkshire, UK: a basis for groundwater quality management[END_REF] (Fig. 4). The groundwater showed higher Ca, Mg, and Sr concentrations than would be expected from an ideal mixing with seawater and local rainwater. The dots in all the diagrams plot above the seawater dilution line, except for some E-H.A wells that showed a slight Na depletion. Since no evaporites have been described in the study region, the groundwater Cl can only originate from rainfall recharge and human activities, mainly through the application of fertilizers and wastewater inputs [START_REF] Widory | Nitrate in groundwater: an isotopic multi-tracer approach[END_REF]. The abundance of Ca due to carbonate dissolution masks the Ca excess that could be combined with the Na deficit due to cationic exchange with the clays in the alluvial sediment above the Chalk. Sr and Mg gave similar variations when plotted against Cl, confirming the similarity between the behavior of Mg 2+ and Sr 2+ [START_REF] Dogramaci | Strontium and carbon isotope constraints on carbonate-solution interactions and inter-aquifer mixing in groundwaters of the semi-arid Murray Basin, Australia[END_REF].

Sr isotope systematics for defining lithological controls on groundwater chemistry

The primary sources of Sr in groundwater are atmospheric input, dissolution of Sr-bearing minerals, and anthropogenic input. The concentration of Sr across the aquifer is usually closely related to the congruent and incongruent reactions of calcite. The biogenic calcite of coccoliths, which constitute the main phase of the Chalk, dissolves and brings Mg and Sr as trace elements to the groundwater, whilst a purer calcite re-precipitates [START_REF] Edmunds | Baseline geochemical conditions in the Chalk aquifer, Berkshire, UK: a basis for groundwater quality management[END_REF]. The highest Sr concentration found in the confined part of the Chalk aquifer (53.8 µmol/L) made water just below strontianite and celestite saturation (SI = -0.75 and -0.97).

Strontium-isotope ratios vary in nature because one of the strontium isotopes ( 87 Sr) is formed by the radioactive decay of the naturally occurring element rubidium ( 87 Rb). The 87 Sr/ 86 Sr ratio is classically used as a tracer of water-rock interaction [START_REF] Blum | 87Sr/86Sr ratios of Sierra Nevada stream waters: implications for relative mineral weathering rates[END_REF][START_REF] Négrel | Strontium isotope systematics used to decipher the origin of groundwaters sampled from granitoids: the Vienne Case (France)[END_REF]. Strontium isotopic signatures in groundwater are not dependent on dilution processes (e.g. hydrological conditions) and so enable sources of dissolved Sr to be traced. Here normalized-Sr ratios, the Ca/Sr against the Mg/Sr ratios, are first investigated for both solids and the groundwater (Fig. 5), as well as the 87 Sr/ 86 Sr ratio against the Ca/Sr, Mg/Sr and Se/Sr ratios (Fig. 6, 7,and 8) where mixing is represented by straight lines [START_REF] Faure | Principles of Isotope Geology[END_REF].

The Ca/Sr and Mg/Sr ratios indicate a mixing of groundwater with different chemistries linked to lithological variations, especially between the carbonate and silicate poles, the presence or absence of dolomite, and variations in residence time. In the Ca/Sr vs. Mg/Sr diagram for solids (Fig. 5a), samples of the regional geological formations constitute two contrasting end- In the Ca/Sr vs. Mg/Sr diagram for water (Fig. 5b), the Ca/Sr and Mg/Sr ratios of the studied groundwater were in the range of the Chalk groundwater of the London Basin in Berkshire (respectively 17-5300 and 14-130; [START_REF] Edmunds | Baseline geochemical conditions in the Chalk aquifer, Berkshire, UK: a basis for groundwater quality management[END_REF], but displayed a larger range than for the Chalk aquifer of the Somme region (respectively 700-1200 and 30-110; [START_REF] Négrel | Strontium isotopes as tracers of groundwater-induced floods: the Somme case study (France)[END_REF] or at La Bassée (respectively 700-1300 and 20-50; (Gourcy and Brenot, 2011). At the study site, three main end-members are indicated to explain the chemical composition of the groundwater:

(1) The E-H.A group, with the highest Ca/Sr ratios that could result from a mixing between the two lithological poles mentioned above.

(2) The groundwater of the confined area represented by Well A26 displayed the lowest Ca/Sr and Mg/Sr ratio values. The ratios were similar to those measured in the confined parts of the Chalk aquifer in Berkshire [START_REF] Edmunds | Baseline geochemical conditions in the Chalk aquifer, Berkshire, UK: a basis for groundwater quality management[END_REF] and also to those measured in the confined Carboniferous groundwater sampled nearby [START_REF] Gourcy | Apport des outils chimiques et isotopiques dans la connaissance de l'aquifère du calcaire carbonifère[END_REF].

(3) The groundwater from Well A12 with the highest Mg/Sr ratios plotted correctly on a mixing line between the Chalk and the T-Q pole. Well A17 displayed a similar composition than that of Well A12.

All the groundwater samples showed a trend from low (10) to high (1000) Mg/Sr ratios between the Well A26 and Well A12 type end-members. They clearly displayed a varied contribution of confined groundwater, which is consistent with the proximity of the Tertiary formations as well as with the piezometric fluxes. The Mg/Sr ratio successfully discriminates mixing of groundwater from Les Ansereuilles with old groundwater with a long residence time enhancing the water-rock interactions.

Surface water mixing with the groundwater is also possible, given the similar value of the Ca/Sr ratio in the groundwater of the southwestern part of the Les Ansereuilles catchment area (wells A23, A24 [annual mean], A17, A18 and A3) and the surface water of the Deûle River. The Deûle River may indeed represent an end-member of the system due to its contribution (900 m 3 /d) to the groundwater.

Significant monthly variations were recorded for the Ca/Sr and Mg/Sr molar ratios in monitored wells A13 and A24 whereas they were slightly less for wells H4 and S3

(Supporting Information Table 1). This indicates a change in the water chemistry of the four well fields during the period of monitoring resulting from degrees of mixing between groundwater originating from the confined and unconfined parts of the aquifer.

In the 87 Sr/ 86 Sr vs. Ca/Sr diagram (Fig. 6), a global trend is seen between two contrasted poles with a large Ca/Sr and 87 Sr/ 86 Sr variation: (1) the E-H.A groundwater (E10, E5, E7, H4, H9) with the high 87 Sr/ 86 Sr values (0.70840-0.70845), and (2) the A26 confined groundwater with the lowest strontium isotopic signature (0.70789). The Sr isotopic compositions of these two end-members correspond to lithological signatures that have been identified in a ternary system at regional scale: (1) Cretaceous Chalk leachates (0.70735) related to the oxic conditions, and (2) the Well A26 type in anoxic conditions, and (3) the T-Q pole with the very high 87 Sr/ 86 Sr ratios of Ypresian (0.70977) and Quaternary (0.71-0.73; [START_REF] Drouet | Strontium isotope composition as a tracer of calcium sources in two forest ecosystems in Belgium[END_REF])) rock leachates.

The leachate signature of the Chalk is similar to that of the UK Turonian Chalk (0.7074) [START_REF] Mcarthur | Strontium isotope stratigraphy for the Late Cretaceous: a new curve, based on the English Chalk, High Resolution Stratigraphy[END_REF]. The base of the Ostricourt Sand displayed an intermediate 87 Sr/ 86 Sr ratio (0.70886) which could result from the infiltration of leachate water from the clayey Upper

Ypresian and Quaternary formations. The strontium isotopic values of the groundwater are similar of those of La Bassée (Gourcy and Brenot, 2011).

A deviation towards high 87 Sr/ 86 Sr is observed and the strontium isotopic signature of Well A26 groundwater was slightly more radiogenic than both the local Chalk rock and the Chalk groundwater in Berkshire; a deviation that had already been identified in the Chalk aquifer [START_REF] Brenot | Interaction between different water bodies in a small catchment in the Paris basin (Brévilles, France): Tracing of multiple Sr sources through Sr isotopes coupled with Mg/Sr and Ca/Sr ratios[END_REF]Gourcy and Brenot, 2011;[START_REF] Négrel | Strontium isotopes as tracers of groundwater-induced floods: the Somme case study (France)[END_REF]. The strontium isotopic signature of the local Chalk aquifer can be explained by (1) the mixing of the various Tertiary weathering products with the Chalk component that presents a low Sr isotopic signature, and (2) the major impact of carbonate dissolution on the isotopic signature of the Chalk aquifer due to the dissolution kinetics of carbonates being higher than that of clay or K-and Rb-bearing minerals (illite and muscovite with high isotopic ratio).

The confined groundwater of Well A26 may integrate a slight high 87 Sr/ 86 Sr component from the overlying Sr-radiogenic T-Q silicate end-member. Similarly, in Berkshire, the strontium isotopic variation measured in the Chalk aquifer (0.7071 to 0.7074) could also be explained by leaching of the overlying Eocene Reading Beds (sand, loam, clay...) that may also present more radiogenic values.

The 87 Sr/ 86 Sr ratios of E-H.A groundwater, being higher than those of Well A26, may well be explained by a contribution of high 87 Sr/ 86 Sr from fertilizers [START_REF] Brenot | Interaction between different water bodies in a small catchment in the Paris basin (Brévilles, France): Tracing of multiple Sr sources through Sr isotopes coupled with Mg/Sr and Ca/Sr ratios[END_REF][START_REF] Vitoria | Fertilizer characterization: Isotopic data (N, S, O, C, and Sr)[END_REF]. Correlating NO 3 and Sr shows that supplementary Sr with a more radiogenic signature was provided to the system in the E-H.A recharge area where a high nitrate content from the use of fertilizers was measured [START_REF] Cary | Multi isotopic tools to understand selenium origins in groundwaters of the Chalk aquifer in Northern France[END_REF]. The E-H.A pole is in line with the A26-A12 correlation.

In the fault depression of the northeastern part of the Les Ansereuilles catchment area, wells A14 and A12 yielded Ca/Sr and strontium isotopic ratios that represented the entire variation range of Les Ansereuilles (apart from Well A26). They constitute two clear local endmembers among which all points of Les Ansereuilles (minus Well A26) were well distributed.

In the Salomé well field, the groundwater of wells S4b and S8 diverged from the mixing line with lower Sr values and lower Sr isotopic signatures that are closer to the isotopic signature of the Chalk aquifer of the Somme [START_REF] Négrel | Strontium isotopes as tracers of groundwater-induced floods: the Somme case study (France)[END_REF] or La Bassée (Gourcy and Brenot, 2011). This isotopic variation can be explained by a contribution of slightly different groundwater closer to the local Chalk end-member solid leachates. According to the piezometric maps, the Salomé groundwater incorporates groundwater from the southwest and may also receive surface water from the Aire canal that could have a composition distinct from that of the Deûle River. In the clayey area in the northeast of the Les Ansereuilles, the difference in the 87 Sr/ 86 Sr ratios of wells A12 (0.70826), A13 (0.70811) and A14 (0.70806), which are 200 m apart, suggests (i) the contribution of a more radiogenic component through desorption of Sr from the upper radiogenic T-Q end-member linked to the hydrological conditions, and (ii) a decreasing contribution of confined groundwater from wells A14 to A12. Given the geological context, the input of water from a different origin flowing along the Seclin Fault can also be expected and would have a greater impact on the groundwater of Well A14 (Fig. 2).

However, if these fluxes had a chemical composition similar to the groundwater of Well H11, then they would bring a very low Sr content (0.01 µmol/L) with a strontium isotopic signature (0.70812) similar to that of Well A13 and would be of minor consequence. The wells in the central (A3-A4-A16) and southern (A17 to A24) parts of Les Ansereuilles displayed a lower variability than for the northern part and also a decreasing influence of the confined groundwaters of A26 type compensated by a higher contribution of the E-H.A end-member in coherence with the piezometric fluxes (Fig. 2).

Interestingly, data from the Upper Ypresian Cuisian aquifer in the Paris Basin (Brévilles catchment) displayed similar and higher Sr isotopic values and Ca/Sr ratios [START_REF] Brenot | Interaction between different water bodies in a small catchment in the Paris basin (Brévilles, France): Tracing of multiple Sr sources through Sr isotopes coupled with Mg/Sr and Ca/Sr ratios[END_REF]. The Cuisian groundwater also plots on the mixing line between Well A26 and the T-Q end-member. This relationship can be explained by the fact that the Cuisian aquifer is bound by the clayey Ypresian (Sparnacian), leaching of which could supply radiogenic Sr to the groundwater.

Sr isotopes to constrain Se origins

The significant variations of the Sr content and the 87 Sr/ 86 Sr ratio in groundwater during the monitoring year suggest various water-body pathways and mixings. The Sr isotopes were then used to better elucidate the Se origins.

The 87 Sr/ 86 Sr vs. Se/Sr diagram (Fig. 8) again reveals the three end-members of the ternary system: (1) The E-H.A groundwater, with the exception of Well H4 due to its higher Se content;

(2) the A26-type confined groundwater which displayed an Se content of less than the quantification limit due to reductive conditions, but still represents an important endmember of the system; and (3), the A12 end-member completed with wells H11 and H4

representing the Se-rich groundwater in the depressed faulted area due to its high Se content and high Se/Sr ratio (0.045).

The T-Q end-member displayed high 87 Sr/ 86 Sr and relatively high Se concentrations. More locally, the Se-rich end-member is related to the groundwater of wells A12, H4, and H11. The variations in the Se/Sr ratio and in the Sr isotopic compositions in the neighboring wells A12, A13 and A14 can be explained by various mixings with groundwater of the A26 or A14 type and with various contributions of radiogenic Sr and Se from the upper clayey formation.

Interestingly, the high Se contents in Well A13 in June and October (respectively 0.33 and 0.23 µmol/L) corresponded to high and similar piezometric levels and were associated with the almost constant Sr concentrations (respectively 15.3 and 17.3 µmol/L) and the Sr isotopic ratios that were exactly equal for the two dates (0.7081). The A13 groundwater evolved along the mixing line between Wells A14 and A12 according to the different kinetics of the main processes affecting Sr and Se mobility, i.e. dissolution/precipitation and sorption.

The Houplin-Ancoisne wells showed distinct strontium isotopic ratios, but similar Se/Sr ratios (0.03) only slightly inferior to that of Well A12. Well H11 thus plots directly on a mixing line between the Well A26 and Well A12 end-members, confirming the influence of confined groundwater. Conversely, the composition of H4 can be explained by a mixing between groundwater from the E-H.A and the A12 end-members. As already discussed, this organization is supported by the presence of the main Seclin Fault which forms a hydraulic screen and provides high groundwater flows along preferential flow paths typically used for exploitation. This pattern is affected by the fact that Se is a redox reactive element that is immobile under reductive conditions and mobile under oxidizing conditions, and these are controlled by the confinement state of the Chalk aquifer and so can be linked to the Mg and Sr concentrations. Thus the very low Se/Sr ratios in the groundwater of wells A26 and A14 are explained by the low redox potential and the confinement conditions.

Focus on Sr and Se mobility in wells A13 and H4: evidence of piezometric level influence

The monitoring data have been completed with data acquired since 2000. Only the chemical monitoring in wells A13 and H4, which showed significant variations, is given here.

Taking the case of Well A13 in Les Ansereuilles, it is worth noting that (1) the Se concentrations have exceeded the standard for drinking water (0.12 µmol/L) from 2005 to 2013, and (2) the static piezometric groundwater level has been located in the sediments overlying the Chalk (Fig. 9). The previous period of high piezometric levels corresponds to the rather famous 2001-2003 Somme floods [START_REF] Négrel | Strontium isotopes as tracers of groundwater-induced floods: the Somme case study (France)[END_REF]. The monitoring data showed that the variability of Se concentrations correlated well with variations in piezometric level and dissolved oxygen, except after 09/2012 when dissolved oxygen remained high (4-5 mg/L) while the Se content decreased; the high oxygen level could have been a consequence of the anomalous 2012 spring-summer recharge. The correlation between the Se content and the piezometric level is explained by the presence of Se in the clayey overlying formation [START_REF] Cary | Tracking selenium behaviour in chalk aquifer (northern France): Sr and 34S-sulphates isotopes constraints[END_REF]. Furthermore, the presence of Se in the A11-A15 particular geological context is a second element that corroborates the hypothesis of a Tertiary age for the clayey formation in the WellA11-A15 area at 5-15 m depth (Fig 2). The strontium content showed slight variations opposed to those of Se (-0.64) and the piezometric levels, with the lowest Sr values recorded during the anomalous period of high piezometric level in the summer (due to a very rainy spring and early summer).

During this period of microporal saturation, oxygen and nitrates caused oxidation of the Sebearing phases present in the clayey deposits [START_REF] Bailey | Simulating variably-saturated reactive transport of selenium and nitrogen in agricultural groundwater systems[END_REF][START_REF] Bailey | The Influence of nitrate on selenium in irrigated agricultural groundwater systems[END_REF] and that the clays favored Sr adsorption. Conversely, during periods of lower piezometric level localized in the Chalk aquifer, the clayey formation was not leached by groundwater, thus inducing less Se mobilization and a higher Sr content in water.

The SF6 concentrations similarly showed a reverse correlation with the piezometric levels and the Se content (-0.73), except for the higher July value in Well A13. The highest SF6 contents were measured during winter where very low piezometric levels were recorded.

These low piezometric levels can have triggered infiltration of recharge water, especially surface water from the Deûle River. The BMM correctly represented this evolution with a higher fraction of modern water in the winter and a mean 35% of modern water during the spring and summer seasons, except for the contaminated measurement in July 2012 (Fig. 9).

The end-members identified with the isotopic strontium typology are again discriminated when plotting the percentage of modern water against the Se concentrations (Fig. 10).

Because a time lag can be expected for the diffusive transport of SF6 through the unsaturated zone of a deep poorly porous aquifer, the apparent tracer age will be older than the true age of the recharge [START_REF] Gooddy | Using chlorofluorocarbons (CFCs) and sulphur hexafluoride (SF6) to characterise groundwater movement and residence time in a lowland Chalk catchment[END_REF]. The old water almost certainly originated from a mixing with old confined water of the A26 type and with leachates from the overlying clayey formation.

Taking the case of Well H4 in Houplin-Ancoisne, we again see an apparent increase of Se appeared between Se and the piezometric level nor with SF6, probably because the overlying deposits were constantly saturated. The Se content of the overlying clayey sediments is inferior to the one of the clays of Well A13 (Cary et al., in prep).

Finally, the drilling operation which requires daily stops and ignition of pumping impacts the variations of the dynamic level and thus alternating saturation or desaturation of the overlying clays. Strong and daily variations of the dynamic level of the water most surely contribute secondarily to mobilize the Se sources in the saturated zone. In order to limit Se contamination of the groundwater, the water operators have already decreased the pumping rates in the depressed area of wells A11-A15. A complete study of the spatial distribution of the Tertiary and Quaternary formations in the study area should enable new wells to be drilled away from them.

Conclusions

The specific geochemical signature of the Chalk groundwater south of Lille revealed by Sr isotopes, in addition to the major-and trace-element concentrations demonstrated a mixing of groundwater bodies where the main poles could be represented by: (1) the upstream groundwater in the recharge area with a higher 87 Sr/ 86 Sr ratio, probably controlled by fertilizer inputs, and SF6 contamination;

(2) the confined groundwater with a high Sr concentration and the lowest 87 Sr/ 86 Sr, Ca/Sr and Mg/Sr ratios of the study area, along with a higher contribution of pre-tracer water and a very low Se content;

(3) the Se-rich groundwater. The four monitored wells tapping the Chalk aquifer showed strong temporal and spatial disparities as regards Sr and Se concentrations during a 21-month study.

The spatial and temporal disparities of the Se concentrations in this dynamic system can be explained by the geological and hydrogeological constraints, i.e. the presence of a Se-rich clayey formation identified as Tertiary overlying the Chalk aquifer, and variations of the piezometric level induced by natural recharge and pumping. Se mobility is enhanced by saturation and oxidation of the reductive clayey formation that locally confines the Chalk aquifer by SF6-contaminated groundwater. Strontium could, at the same time, be favorably adsorbed on clays. Pumping of the groundwater may also favor Se transfer towards piezometrically depleted areas.

In order to limit Se contamination in the groundwater, the geology of the study area should be revisited, so that new wells should be drilled in areas devoid of Se-rich clayey deposits, and connections with Se-rich wells should be as limited as possible.

As a first approach, this multi-tool study successfully identified Se origins in the Chalk aquifer. There is now a need to explore Se-bearing phases in the Tertiary and Quaternary sediments and to better constraint the oxido-reductive conditions that control Se mobility and Se speciation, and especially the role of denitrification. This will help to a better exploitation of the groundwater resource.

Figure caption

Figure 1: Se content (µg/L) in groundwater (scale 1:1,000,000) affecting the southern and eastern parts of the Paris Basin due mainly to the Tertiary formations (data from the ADES database www.ades.eaufrance.fr). [START_REF] Bernard | Acquisition et rassemblement des données géothermiques disponibles et nouvelles dans la zone Franco-Belge de Saint-Ghislain à Saint-Amand-les-Eaux (Nord)[END_REF][START_REF] Simon | Apports de la biogéochimie isotopique à la connaissance de la dénitrification dans les eaux souterraines[END_REF][START_REF] Vachier | Mouvements de l'eau dans la zone non saturée et alimentation de la nappe de la craie de Champagne (France)[END_REF] and locally [START_REF] Gourcy | Apport des outils chimiques et isotopiques dans la connaissance de l'aquifère du calcaire carbonifère[END_REF], plus local Carboniferous data, are also plotted.

Figure 4: Diagrams showing element concentrations in the groundwater and surface water of the study area, including results from the initial sampling and analysis and the monthly monitoring programs (the Sr vs. Cl and Se vs. Cl graphs both give all the monitoring values for wells A13 and H4). Data from ADES (*) are given in black crosses for confined groundwaters and grey crosses for A26 data when available.

Figure 5: a) Bedrock end-members: Ca/Sr vs. Mg/Sr (mol/mol) in the solids sampled in the study region and at depth close to the wells. The Ypresian sample is from Montreuil-sur-Epte (Gourcy and Brenot, 2011) and the carbonates given by [START_REF] Meybeck | Composition chimique des ruisseaux non pollués de France[END_REF]. b) Groundwater end-members: Ca/Sr vs. Mg/Sr (mol/mol) in groundwater sampled from the study site. The groundwater data are from Gourcy and Crastes de Paulet (2012) for the Carboniferous and local Chalk groundwater, [START_REF] Négrel | Strontium isotopes as tracers of groundwater-induced floods: the Somme case study (France)[END_REF] for the Somme catchment Gourcy and Brenot (2011) for la Bassée, [START_REF] Brenot | Interaction between different water bodies in a small catchment in the Paris basin (Brévilles, France): Tracing of multiple Sr sources through Sr isotopes coupled with Mg/Sr and Ca/Sr ratios[END_REF] for the Cuise Sand, [START_REF] Roy | Geochemistry of dissolved and suspended loads of the Seine river, France: Anthropogenic impact, carbonate and silicate weathering[END_REF] for the Oise River, and [START_REF] Edmunds | Baseline geochemical conditions in the Chalk aquifer, Berkshire, UK: a basis for groundwater quality management[END_REF] for the Chalk aquifer in the London Basin. Errors are within the points. Mean values were calculated for the 2011, 2012 and 2013 samples.

Figure 6: Ca/Sr (mol/mol) vs. strontium isotopes in the sampled groundwater of the study area. Groundwater data from the bibliography were also plotted for the Chalk aquifer in Berkshire (London Basin [START_REF] Edmunds | Baseline geochemical conditions in the Chalk aquifer, Berkshire, UK: a basis for groundwater quality management[END_REF], the Brévilles catchment [START_REF] Brenot | Interaction between different water bodies in a small catchment in the Paris basin (Brévilles, France): Tracing of multiple Sr sources through Sr isotopes coupled with Mg/Sr and Ca/Sr ratios[END_REF], La Bassée (Gourcy and Brenot, 2011) and the Somme catchment [START_REF] Négrel | Strontium isotopes as tracers of groundwater-induced floods: the Somme case study (France)[END_REF]. The strontium isotopic signature of local atmospheric inputs and Quaternary silt, together with the Ca/Sr variation of atmospheric inputs, are from [START_REF] Drouet | Strontium isotope composition as a tracer of calcium sources in two forest ecosystems in Belgium[END_REF]. The groundwater plots on a mixing line between two poles in the Chalk aquifer and were inscribed in a regional ternary system defined by the silicate and carbonate poles, the latter evolving between free and confined conditions that imply varying degrees of water-rock interaction.

Figure 7: Mg/Sr (mol/mol) vs. strontium isotopes in the sampled groundwater of the study site. The regional system is defined by a very radiogenic silicate end-member and by the low isotopic composition of the solid Chalk. The diagram shows three poles of mixing that better define the Sr isotopic composition of the Chalk aquifer groundwater: (1) the Emmerin (E-H.A) end-member pole, (2) the A26 confined pole, and (3) the A12 end-member pole that is clearly different from the other Chalk poles. Data for the Somme catchment [START_REF] Négrel | Strontium isotopes as tracers of groundwater-induced floods: the Somme case study (France)[END_REF], the Cuise Sand in the Brévilles catchment [START_REF] Brenot | Interaction between different water bodies in a small catchment in the Paris basin (Brévilles, France): Tracing of multiple Sr sources through Sr isotopes coupled with Mg/Sr and Ca/Sr ratios[END_REF], and the Chalk aquifer in the London Basin [START_REF] Edmunds | Baseline geochemical conditions in the Chalk aquifer, Berkshire, UK: a basis for groundwater quality management[END_REF] are indicated. Errors are within the points.

  ), i.e. the Emmerin, Houplin-Ancoisne, Les Ansereuilles and Salomé fields. The Deûle River and other small tributaries (Seclin canal, Tortue and Naviette) drain the area and quantitative and qualitative pressures of industrial, urban and agriculture origins occur. Exploitation of the aquifer in the study area decreased from 40-50 Mm 3 in the 1960-1980's to 30-40 Mm 3 due to decreased industrial pumping. The potentiometric surface of the Chalk layer is distorted by a cone of depression centered on Les Ansereuilles (Fig. 2). In 1985, around 20% of the water for Lille came from the Deûle River via the Quaternary alluvium, but the high pumping rate

  first initial sampling and analysis program was carried out in June 2011 in which the geochemistry and strontium isotope ratios were determined for 20 wells all acting as drinking water wells. A second program was a monthly monitoring program over a period of 14 months(2012)(2013). It aimed at studying four wells (A13, A24, H4 and S3), yielding a range of Se contents corresponding to the high variability observed in the well fields, and the monitoring of the Deûle River and of its tributaries. Measurements of the physical and chemical parameters (temperature, conductivity, pH, Eh, and dissolved oxygen) were made in situ. The groundwater samples were then filtered in the field at 0.45 µm. Samples for cations and Sr isotopes analyses were acidified by HNO 3 . Samples for anions were stored non acidified.Analyses of Ca, Mg, Na, K were carried out by inductively coupled plasma-atomic emission spectroscopy (ICP-AES, Ultima-2 model, Jobin Yvon), and Cl, SO 4 , and NO 3 by ion chromatography (model DX120, Dionex) in the BRGM laboratory. The accuracy of both techniques was around 5-10% depending upon the concentration. Selenium and strontium concentrations were measured by inductively coupled plasma mass spectrometry (ICP-MS) (uncertainty 5%). Alkalinity was determined by HCl titration and Gran's method on the sampling day.

  members: (1) the local Chalk end-member with high Ca/Sr (1500-1620 mol/mol) and low Mg/Sr (18 mol/mol) ratios, and (2) the T-Q silicate end-member with low Ca/Sr and high Mg/Sr ratios. The Mg/Sr ratio of the Ostricourt sand evolved from 209 at the bottom to 1260 mol/mol at the top, showing enrichment in Mg probably due to leachates from the Upper Ypresian clay.Within the four studied wells, the chemical composition of the solids evolved clearly between these two lithological end-members. With similar Mg/Sr values, the sediment of Well H4 showed a higher Ca/Sr ratio than of Well A13, the latter being almost dominated by clay whereas the former comprises successive horizons of sand and clay, in which Chalk grains can be found, inherited from alternating Quaternary erosion/sedimentation.

  Figure 6with high Mg/Sr ratios (200-600) and high 87 Sr/ 86 Sr values, especially for the Ypresian and Quaternary components. The groundwater of Well A12 near the local Se-rich end-member may typically result from a mixing between T-Q leachates and confined groundwater as found at the confined pole of Well A26. This diagram also highlights the E-H.A recharge pole with lower Mg/Sr and higher 87 Sr/ 86 Sr values, which is clearly distinct from the A26-A12 mixing line. In the study site, however, the confined groundwater has high Sr and Mg concentrations, close to those of the Chalk rock, with a low Mg/Sr ratio. Consequently, the Mg and Sr concentrations in the groundwater of the recharge and

  concentrations recorded from 2000 to 2013. Considering 2005-2013 data, the Se content close to the drinking water standard is clearly linked to the groundwater level rising above the Chalk formation. Here, the Se and Sr contents were well correlated (R 2 =0.92) and the Se contents were negatively correlated with dissolved O 2 (R 2 =-0.90). No clear correlation

Figure 2 :

 2 Figure2: Study site showing the geological formations, the four studied well fields and piezometric levels of the Chalk groundwater (geological scale 1:50,000), and the maxima of Se contents measured in the groundwater in 2011. Quaternary deposits cover the Tertiary formations (Louvil Clay, Ostricourt Sand, Flandres Clay) and the Cretaceous Chalk. The schematic A-B geological cross section is based on geological data (infoterre.brgm.fr) and provides new insights concerning the faulted depressed area (wells A11 to A15) the age of which could be Early Tertiary instead of Quaternary.

Figure 3 :

 3 Figure 3: Stable 2 H and 18 O isotopes of the water molecule measured during the initial sampling and analysis program in June 2011, along with calculated 2012 summer and winter means of the Deûle River. Data from the Chalk aquifer, both regionally[START_REF] Bernard | Acquisition et rassemblement des données géothermiques disponibles et nouvelles dans la zone Franco-Belge de Saint-Ghislain à Saint-Amand-les-Eaux (Nord)[END_REF][START_REF] Simon | Apports de la biogéochimie isotopique à la connaissance de la dénitrification dans les eaux souterraines[END_REF][START_REF] Vachier | Mouvements de l'eau dans la zone non saturée et alimentation de la nappe de la craie de Champagne (France)[END_REF] and locally[START_REF] Gourcy | Apport des outils chimiques et isotopiques dans la connaissance de l'aquifère du calcaire carbonifère[END_REF], plus local Carboniferous data, are also plotted.

Figure 8 :

 8 Figure 8: Se/Sr (mol/mol) vs. strontium isotopes in the sampled groundwater of the study site. This diagram shows the local Se-rich end-member of Well A12 and the chemical relationships with the groundwater of wells H11 and H4. The three poles of the Chalk aquifer are still well defined. Errors are within the points.

Figure 9 :

 9 Figure9: Evolution of the selenium, strontium, dissolved oxygen and SF6 concentrations in the groundwater of a) Les Ansereuilles Well A13, and b) Houplin-Ancoisne Well H4, along with pumping rates and piezometric levels versus time. Chemical data prior to 2011 as well as the pumping rates and piezometric levels are from the water operators.

Figure 10 :

 10 Figure 10: Percentage of young waters vs. Se content in the studied groundwaters. For monitoring in well A13, months are indicated. The end-members identified with the isotopic strontium typology are again well discriminated.

δ 18 O, δ 2 H, and 87 Sr/ 86 Sr variations in the surface water and groundwater

  18 O and δ 2 H values. Apart from the February and March 2012 data, the monitoring in wells A13, A24, S3, and H4 showed only a narrow range of δ 18 O and δ 2 H signatures; the differences during the sampling period were less than 0.6‰ for δ 2 H and 0.07‰ for δ 18 O.The strontium isotopic signatures, only analysed for samples from the initial program, displayed significant spatial disparities. The highest values (0.70840-0.70845) were recorded in Emmerin and Houplin-Ancoisne groundwater (wells E10, E5, E7, H4, H9), whereas the lowest value (0.70789) was recorded from the confined groundwater of Well A26. All the other wells displayed intermediate Sr isotopic values.

	in the E-H.A group. The surface water of the Deûle catchment displayed a low Cl For the groundwater, the δ 18 O and δ 2 H from the initial program fell respectively in the range
	concentration ranging from 1.7 to 2.1 mmol/L. of -7‰ to -5.64‰ and -47.9‰ to -44.2‰ (Fig. 3). The groundwater in Well A26 displayed the
	lowest δ
	The Se content has been measured for control for two decades and high Se content in the
	groundwater was observed since the 2000's. Sixteen of the 27 wells (60%) in Les
	Ansereuilles and 8 of the 11 wells (70%) in Houplin-Ancoisne displayed Se contents >10
	μg/L (0.13 µmol/L). Interestingly, the groundwater of the depressed area in Les Ansereuilles
	displayed a varying Se content ranging from below the quantification limit in Well A14 to 30
	µg/L (0.38 µmol/L) in Well A12 where pumping was subsequently abandoned. Consequently,
	we selected the nearby (100 m) pumped Well A13 for the monitoring program. Any Se in the
	groundwater sampled from wells A26 and A14 was below the detection limit. Well A24, with a
	range of Se content of 10-15 µg/L (0.12-0.19 µmol/L) that was inferior to the drinking water
	standard, was monitored as a reference. In the Salomé well field, 5 out of the 17 wells
	displayed a Se content higher than the drinking water standard, whereas only one of the nine
	wells in the Emmerin well field returned an Se concentration of >10 µg/L.
	Monthly regular sampling of the groundwater displayed evolutionary trends in the major-and
	trace-element concentrations for wells A13, H4, A24 and S3. With wells A13 and H4, this
	basically followed variations in the hydrological conditions(high rainfalls in April (92 mm),
	June (82 mm), July (150 mm), and October (132 mm) 2012). For Well A13, the Sr content
	varied from 11.5 to 18.0 µmol/L respectively in July and April 2012, and the Se content varied
	significantly between 0.16 µmol/L (12.8 µg/L) in January 2013 and 0.4 µmol/L (32 µg/L) in
	July 2012, with an average value of 0.27 µmol/L. Well H4 also displayed varying Se
	concentrations, although with a lower range (0.06 to 0.17 µmol/L) and a slightly lower
	average value (0.25 µmol/L for 2012). Monitored as a reference, Well A24 showed no
	significant variation in Sr and Se concentration (respectively mean of 7.5 and 0.16 µmol/L).
	Similarly, Well S3 in Salomé displayed only a very low variation in its chemistry (respectively
	mean of 15 for Sr and 0.07 µmol/L for Se).
	4.2 For the surface water, the stable isotopes ranged from -7.76‰ to -4.8‰ for δ 18 O and from -
	51.8‰ to -35.6‰ for δ 2 H (see Supporting Information Table 1). The average δ 18 O and δ
	2 Hdata for the Deûle River ranged respectively from -6.84‰ and -46.50‰ in winter to -6.2‰
	and -44.06‰ in summer. The values for the Seclin canal and the Deûle, Naviette and Tortue
	river waters were particularly enriched in δ 2 H and δ 18 O during the 2012-2013 winter.

Acknowledgments

This study was financed by the Lille urban community (LMCU), the French Geological Survey (BRGM), and the water operators Noréades and Eaux du Nord. The Research Direction of the BRGM is thanked for its additional support. The work benefited from the collaboration of C. Guerrot, C. Fléhoc and A.M. Gallas (BRGM Mass-Spectrometry team) who provided the Sr and O-H isotopic analyses. It also benefited from the collaboration of T. Conte (BRGM Chemistry laboratories) who provided the major-and trace-element analyses. Jean-Yves

Caous is thanked for fruitful discussions. The water operators are thanked for their technical help in the field and for all the chemical records at our disposal. We are grateful to Patrick Skipwith for the English review. The authors are grateful to Kimberly Beisner and an anonymous reviewer for thoughtful revisions of the manuscript and constructive comments.

Supporting Information

Table 1: Physical and chemical analyses of the surface water and groundwater from the study area. Deûle data are for the Deûle River, and among groundwater, all samples names of wells with A are for Les Ansereuilles, with E for Emmerin, with H for Houplin-Ancoisne and S for Salomé.