Abstract: Abstract Sloping coastal aquifers in reality are ubiquitous and well documented. Steady state sharp-interface analytical solutions for describing seawater intrusion in sloping confined and unconfined coastal aquifers are developed based on the Dupuit-Forchheimer approximation. Specifically, analytical solutions based on the constant-flux inland boundary condition are derived by solving the discharge equation for the interface zone with the continuity conditions of the head and flux applied at the interface between the freshwater zone and the interface zone. Analytical solutions for the constant-head inland boundary are then obtained by developing the relationship between the inland freshwater flux and hydraulic head and combining this relationship with the solutions of the constant-flux inland boundary. It is found that for the constant-flux inland boundary, the shape of the saltwater interface is independent of the geometry of the bottom confining layer for both aquifer types, despite that the geometry of the bottom confining layer determines the location of the interface tip. This is attributed to that the hydraulic head at the interface is identical to that of the coastal boundary, so the shape of the bed below the interface is irrelevant to the interface position. Moreover, developed analytical solutions with an empirical factor on the density factor are in good agreement with the results of variable-density flow numerical modeling. Analytical solutions developed in this study provide a powerful tool for assessment of seawater intrusion in sloping coastal aquifers as well as in coastal aquifers with a known freshwater flux but an arbitrary geometry of the bottom confining layer.

Abstract: The Roussillon sedimentary Basin (South France) is a complex multi-layered aquifer, close to the Mediterranean Sea facing seasonally increases of water abstraction and salinization issues. We report geochemical and isotopic vertical variability in this aquifer using groundwater sampled with a Westbay System® at two coastal monitoring sites: Barcarès and Canet. The Westbay sampling allows pointing out and explaining the variation of water quality along vertical profiles, both in productive layers and in the less permeable ones where most of the chemical processes are susceptible to take place. The aquifer layers are not equally impacted by salinization, with electrical conductivity ranging from 460 to 43,000μS·cm−1. The δ2H–δ18O signatures show mixing between seawater and freshwater components with long water residence time as evidenced by the lack of contribution from modern water using 3H, 14C and CFCs/SF6. S(SO4) isotopes also evidence seawater contribution but some signatures can be related to oxidation of pyrite and/or organically bounded S. In the upper layers 87Sr/86Sr ratios are close to that of seawater and then increase with depth, reflecting water–rock interaction with argillaceous formations while punctual low values reflect interaction with carbonate. Boron isotopes highlight secondary processes such as adsorption/desorption onto clays in addition to mixings. At the Barcarès site (120m deep), the high salinity in some layers appear to be related neither to present day seawater intrusion, nor to Salses-Leucate lagoonwater intrusion. Groundwater chemical composition thus highlights binary mixing between fresh groundwater and inherited salty water together with cation exchange processes, water–rock interactions and, locally, sedimentary organic matter mineralisation probably enhanced by pyrite oxidation. Finally, combining the results of this study and those of Caballero and Ladouche (2015), we discuss the possible future evolution of this aquifer system under global change, as well as the potential management strategies needed to preserve quantitatively and qualitatively this water resource.