Articles | Volume 26, issue 5
Research article
10 Mar 2022
Research article |  | 10 Mar 2022

Saline groundwater evolution in the Luanhe River delta (China) during the Holocene: hydrochemical, isotopic, and sedimentary evidence

Xianzhang Dang, Maosheng Gao, Zhang Wen, Guohua Hou, Hamza Jakada, Daniel Ayejoto, and Qiming Sun

Since the Quaternary Period, paleo-seawater intrusions have been suggested to explain the observed saline groundwater that extends far inland in coastal zones. The Luanhe River delta (northwest coast of the Bohai Sea, China) is characterized by the distribution of saline, brine, brackish, and fresh groundwater from the coastline inland. The groundwater in this region exhibits a wide range of total dissolved solids (TDS): 0.38–125.9 g L−1. Meanwhile, previous studies have revealed that this area was significantly affected by Holocene marine transgression. This study used hydrochemical, isotopic, and sedimentological methods to investigate groundwater salinization processes in the Luanhe River delta and its links to paleo-environmental settings. Isotopic results (2H, 18O, 14C) allowed old groundwater recharge to be distinguished from new groundwater recharge. Hydrochemical analysis using the PHREEQC code indicated that the salt in saline and brine groundwater originates from a marine source. The 18O–Cl relationship diagram yields three-end-member groundwater mixing, and two mixing scenarios are suggested to explain the freshening and salinization processes in the study area. When this was interpreted along with data from paleo-environmental sediments, we found that groundwater salinization may have occurred since the Holocene marine transgression. The brine is characterized by radiocarbon activities of  50–85 pMC and relatively depleted stable isotopes, which are associated with seawater evaporation in the ancient lagoon during delta progradation and mixing with deeper fresh groundwater, which was probably recharged in the cold Late Pleistocene. The brackish and fresh groundwaters are characterized by river-like stable isotope values, where high radiocarbon activities (74.3–105.9 pMC) were formed after the washing out of the salinized aquifer by surface water in the delta plain. This study presents an approach that utilizes geochemical indicator analysis with paleo-geographic reconstruction to better assess groundwater evolutionary patterns in coastal aquifers.

1 Introduction

It is estimated that around 40 % of the world's population live in coastal areas (UN Atlas, 2010). Groundwater is an important freshwater resource for domestic consumption and agricultural activities in this region (Cary et al., 2015; Jayathunga et al., 2020). However, groundwater salinization poses a significant threat to everyday living and development activities (Tulipano, 2005; de Montety et al., 2008). In recent decades, groundwater salinization in coastal zones has prompted widespread concern and has been widely studied. On the one hand, seawater intrusion due to groundwater pumping is an important salinization process in coastal aquifers (Reilly and Goodman, 1985; Werner, 2010, 2013; Han and Currell, 2018). On the other hand, groundwater salinization caused by paleo-seawater intrusion in response to Quaternary changes in the global sea level has been reported in many coastal zones worldwide (Edmunds, 2001; Akouvi, 2008; Santucci et al., 2016, Larsen et al., 2017).

Coastal aquifers are linked to the ocean and continental hydrological cycles (Ferguson and Gleeson, 2012), both of which are influenced by natural and human-induced change (Jiao and Post, 2019). In the natural state, there is a steady-state seawater–freshwater interface that extends inland from the coastal line (Costall et al., 2020). During the Quaternary Period, however, sea-level fluctuations on geological timescales have caused the interface to change, allowing seawater intrusion during transgression events and freshwater flushing during glacial low-sea-level periods, as is evident from the hydrochemical characteristics of groundwater in coastal aquifers (Kooi et al., 2000; Sanford, 2010; Aquilina et al., 2015; Lee et al., 2016). In addition, the hypersaline groundwater found in coastal zones, particularly brine groundwater with a salinity 2–4 times that of seawater, cannot be explained solely using a seawater intrusion model (Sola et al., 2014; Han et al., 2020), and paleo-environment settings must be taken into consideration (van Engelen et al., 2019). Some studies, for example, attribute the presence of brine in Mediterranean countries to the evaporation of seawater in the lagoon system during the Holocene transgression (Giambastiani et al., 2013; Vallejos et al., 2018).

The Bohai Sea of northern China was affected by Late Pleistocene transgressive–regressive cycles, which caused paleo-saltwater intrusions with various salinities along the coastal aquifers (Du et al., 2015; Li et al., 2017). Several studies have applied geochemical methods to elucidate the origin of saline groundwater and salinization processes with an anthropogenic influence, including the induced mixing of brine water from adjacent aquifers caused by groundwater overexploitation in Laizhou Bay (Han et al., 2011, 2014; Liu et al., 2017; Qi et al., 2019). However, the association between groundwater salinization (especially brine formation) and paleo-environmental implications are still not clear. Thus, this study applied a range of chemical, isotopic, and sedimentary indicators to examine the groundwater salinization processes in the Luanhe River delta (situated along the northwestern coast of the Bohai Sea) in relation to recharge, salt sources, mixing behavior, and paleo-geographic evolution. The overall goal was to understand the evolutionary pattern of groundwater as influenced by transgression/regression events on a geologic timescale. Our findings are relevant to aquifer remediation activities in the region as well as other similar sedimentary environments around the world.

2 Background on the study area

The study area is located in northeastern Hebei Province, China, on the west coast of Bohai (Fig. 1a). It consists of an alluvial fan and coastal delta bounded by a Holocene maximum transgression line (Xue, 2016). The delta area can be further divided into two parts: the old delta between the Douhe River and the Suhe River and the new delta between the Suhe River and the modern Luanhe River (He et al., 2020). In terms of geomorphology, the study area is inclined to the south and southwest with a slope of about 0.04 ‰–2 ‰. The temperate monsoon climate of the area affects its average annual temperature of 12.5 C and annual rainfall of 601 mm (1956–2010), with 80 % of the annual rainfall occurring between July and September.

2.1 Hydrogeology

The thickness of Quaternary sediments in the study area is about 400–500 m. According to their lithology and hydrogeological characteristics, there are four distinct Quaternary aquifers (Fig. 1b). The first Holocene aquifer (Q4) is a phreatic or semiconfined aquifer with a bottom depth of 15–30 m, is primarily composed of fine sand and slit, and involves fresh, brackish, saline, and brine groundwater (Dang et al., 2020). The second Late Pleistocene aquifer (Q3), the third Middle Pleistocene aquifer (Q2), and the fourth Early Pleistocene aquifer (Q1) have bottom depths of 120–170, 250–350, and 350–550 m, respectively. They are confined aquifers primarily made of medium sand and gravel (Niu et al., 2019). The first aquifer is mainly recharged by meteoric precipitation and lateral infiltration of surface water (Li et al., 2013). The groundwater from the first aquifer is widely extracted for irrigation in the alluvial fan area. The largest salt farm in North China, the Daqinghe Salt Farm, uses shallow brine groundwater for salt production in the delta area, where there is only minor agricultural activity. Except in the alluvial fan area, the circulation between phreatic and confined aquifers is weak. The deep groundwater in the second, third, and fourth aquifers is mainly recharged by a surrounding mountain range and mainly discharged by human pumping (Ma et al., 2014).

Figure 1(a) Location map of the study area. Also shown are the locations of the sampling site and published cores in the Luanhe River delta. Cores LT05, HCZ, BXZK01, BXZK02 and BXZK03 are from He et al. (2020); cores NP05, NP03, DY01, DQ03, DQ04, DY02, MT04, BG07, FB01, A02, and TH04 are from Xu et al. (2020); core LQZ04 is from Cheng et al. (2020); core FG01 is from Xu et al. (2011); core Bai03 is from Li and Wang (1983); and core HCZ is from Peng et al. (1981). (b) Hydrogeological cross section (A–A' in Fig. 1a) of the study area; modified from Ma et al. (2014).

2.2 Sedimentary evolution since the Late Pleistocene

Previous studies have shown that in the study region, the salt–fresh groundwater interface gradually deepens from land (depth of  5 m) to sea (depth of  100 m), as shown in Fig. 1b, with salt groundwater primarily occurring in the first aquifer of the delta area (Li et al., 2013; Ma et al., 2014). According to a stratigraphic transect along the present coastline (Fig. 2), the series stratigraphic architecture of the first aquifer consists of Late Pleistocene continental facies – Holocene marine facies – Holocene delta facies – modern continental facies or artificial fill, indicating that the sediments of the first aquifer have been deposited through lowstand continental accumulation, marine transgression, and highstand progradation since the Late Pleistocene.

The seawater did not reach the modern coastline from the Last Glacial Maximum to the early Holocene (about 30–9 ka BP). The Luanhe alluvial fan was active during this period (He et al., 2020). Since about 9000 a BP, the Holocene marine transgression has approached the present coastline (Xu et al., 2020), and Holocene marine sediments developed due to the sea-level rise from 9–7 ka BP. The Holocene marine transgression had reached its maximum inland area 20 km from the modern coastline by about 7 ka BP (Gao et al., 1981; Peng et al., 1981; Xue, 2014, 2016) (see the Holocene transgression line in Fig. 1). The modern coastal plain was formed through the accumulation of the highstand prograding delta on top of Holocene marine strata, together with artificial fill. In addition, lagoons are important components of the Luanhe River delta (Feng and Zhang, 1998). Based on the records of lagoon facies in the published cores for this region, the approximate distribution range of the buried lagoon is shown as a purple dashed line in Fig. 1a.

Figure 2Stratigraphic transect along the present coastline of the Luanhe River delta; modified from He et al. (2020).

3 Methods

In total, 45 water samples were collected from the Luanhe River delta, including 38 groundwater samples, five surface water samples, one local rain water sample, and one Bohai seawater sample, during four sampling campaigns from October 2016 to June 2020. Groundwater samples were divided into shallow groundwater samples and deep groundwater samples, which were pumped from unconfined and confined aquifers, respectively. The surface water included two Suhe River water samples and two Luanhe River water samples. Due to the artificial fill that has modified the coastal landscape, it was difficult to locate the modern lagoon environment. However, during the investigation, it was found that the Daqinghe Salt Farm in this area extracts seawater into the evaporation pond. The mixture of seawater and meteoric water is subject to evaporation, leading to the formation of concentrated saline water (CSW) in the pond, similar to the formation of CSW in a coastal lagoon (Stumpp et al., 2014). Thus, one CSW sample (sample P18) from the evaporation pond was collected.

Water types were classified according to Zhou (2013) into freshwater (TDS < 1 g L−1), brackish water (TDS = 1–3 g L−1), saline water (TDS = 3–50 g L−1), and brine (TDS > 50 g L−1). Groundwater sampling depths and pH values were measured on site using Eureka Manta+(water-quality multiprobe sonde). The concentrations of K+, Na+, Ca2+, Mg2+, and Br ions were measured using inductively coupled plasma analysis (ICAP-7400), while SO42- and Cl ions were determined using ion chromatography (ICS-600). The HCO3- concentrations of samples were measured using titration. The resulting hydrochemical data are listed in Table S1 (see the Supplement). The stable isotope concentrations (δD, δ18O) in the water samples G02-10, G06-10, G03-05, G04-40, G05-10, G05-46, G07-27, P07-20, P08-30, P09-30, P10-30, P11-20, P12-40 P14-15, P07-100, P13-200, P14-300, P15-150, P16-100, P17-200, P18, LH01, LH02, SH01, SH02, SW01, and R1 were tested at the Experimental & Testing Center of Marine Geology, Ministry of Natural Resources, China, using high-temperature pyrolysis–isotope ratio mass spectrometry. The values of δ18O and δD were calculated with respect to Vienna Standard Mean Ocean Water (VSMOW), and the uncertainties in δD and δ18O were ±1.0 ‰ and ±0.2 ‰, respectively. The radioisotope measurements (AMS 14C) of the groundwater samples P14-300, P15-150, and P16-100 were performed at the Pilot National Laboratory for Marine Science and Technology. Analysis of the stable isotopes (δD, δ18O, 13C) and the radioisotope in groundwater samples G10-10, G03-20, G04-15, G05-30, G06-15, G07-15, G08-15, G08-40, G09-15, G09-40, G10-10, and G10-30 was performed at the Beta Analytic testing laboratory, where the δ18O and δD values were also calculated with respect to VSMOW, and the uncertainties in δD and δ18O are listed in Table S1. The 14C age of groundwater was calculated using the following equation (Clark and Fritz, 1997): bt=-8267lnat14C/qa014C, where t is the radiocarbon age in years Before Present (a BP); at14C is the measured 14C activity in % of modern carbon (pMC); a014C is the derived modern 14C activity of soil; and q is a corrective factor accounting for the dissolution of calcite, which is assumed to be free of 14C and therefore dilutes the initial 14C activity of aqueous dissolved inorganic carbon (DIC) in recharged water. The results for 13C and 14C and the uncorrected residence times are listed in Table S2.

4 Results

4.1 Hydrochemistry

Except for P13-200 (TDS = 1.617 g L−1, indicating brackish water), all the deep groundwater samples in the study region are freshwater. Deep groundwater hydrochemical forms shift from Ca-HCO3 to Na-HCO3 upon moving from the land to the sea (Fig. 3). For the shallow aquifer, the horizontal salt–fresh groundwater interface corresponds more closely to the maximum Holocene transgression line (see Fig. 1a). The Ca-HCO3 type of shallow fresh groundwater is primarily distributed in the alluvial fan region. The brackish and low-TDS saline groundwaters, which vary from Ca-HCO3 to Na-HCO3 and Na-Cl types, occur mainly in the upper aquifer (depths of 0–15 m) of the delta area, while the lower part (depths of 20–40 m) contains Na-Cl-type saline and brine groundwaters with high TDS. Moreover, in terms of the horizontal distribution of salinity, the groundwater TDS tends to decrease from west to east, with the TDS values of saline and brine groundwater generally ranging from 16.57 to 125.97 g L−1 in the old delta (western delta) and from 3.26 to 52.48 g L−1 in the new delta (eastern delta).

Figure 3Piper diagram of the various water samples.


4.2 The stable isotopes 2H and 18O

Figure 4 shows the relationship between deuterium and oxygen-18. The global meteoric water line (GMWL, δ2H = 8 δ18O + 10) is from Craig (1961), while the local meteoric water line (LMWL, δ2H = 6.6 δ18O + 0.3) is based on δ2H and δ18O isotope data (1985–2003, mean monthly rainfall values) from the Tianjin station about 100 km southwest of the study area (IAEA/WMO, 2006). The deep groundwater samples exhibit depleted values of stable isotopes, with values of δ2H ranging from 75.52 ‰ to 57.06 ‰ and those of δ18O from 9.82 ‰ to 7.61 ‰. Shallow groundwater samples have higher hydrogen and oxygen isotope levels, ranging from 64.6 % to 22.46 % for δ2H and 8.74 % to 2.07 % for δ18O. While the relatively small overall values of fresh and brackish groundwater samples are similar to those of the river samples, saline and brine groundwater samples generally plotted below the LMWL or GMWL, meaning that the water was subjected to evaporation prior to becoming groundwater recharge (Gibson et al., 1993), or that multiple-end-member mixing processes were involved (Han et al., 2011).

Figure 4Stable isotope compositions of different water samples. Seawater mixing line: mixing between deep fresh groundwater and seawater; CSW mixing line: mixing between deep fresh groundwater and CSW.


4.3 Groundwater residence times

The measured 14C activities of groundwater samples range from 0.774 to 105.9 pMC (Table S2). A plot of sampling depth versus 14C is shown in Fig. 5; this elucidates a negative correlation, showing that the variation in 14C activity can be attributed to radioactive decay in the aquifer. There are multiple processes that can impact the 14C properties, including groundwater mixing and dispersion, long-term variation in atmospheric 14C, and free 14C dilution (e.g., carbonate dissolution) (Cartwright et al., 2020). Because the relative impacts of these processes are not well established in the study area, the uncertainty regarding the correction of radiocarbon ages to real groundwater ages is very high. Consequently, we estimate the groundwater age as a residence time range. The uncorrected age is considered the maximum age, while the corrected age is considered the minimum age. Corrected ages are determined based on two hypothetical models of carbonate dissolution, which is the main influence on the 14C contents of water samples (Lee et al., 2016).

Figure 5 shows that the 14C activities in the shallow groundwater samples are within 30.6–105.9 pMC. These values indicate a relatively modern recharge before the atmospheric nuclear testing period of the 1950s and 1960s. The radiocarbon activities in the deep fresh groundwater samples are less than 12 pMC, which is consistent with the paleo-water recharge. This indicates that there are weak connections between the shallow and deep aquifers. Therefore, we assume that the shallow aquifer is an open system and the deep aquifer is a closed system. The δ13C mixing and chemical mass balance (CMB) models are used to estimate the corrective factor q (Clark and Fritz, 1997).

For the δ13C mixing model (Pearson and Hanshaw, 1970), q=(δ13CDIC-δ13CCARB)/(δ13CRECH-δ13CCARB), where δ13CDIC is the measured δ13C of DIC in groundwater; δ13CCARB is the δ13C of DIC from dissolved soil minerals (δ13CCARB= 1.5 ‰ from Chen et al., 2003); and δ13CRECH is the δ13C in water when it reaches the saturation zone. In this study, we use a δ13CRECH of 15 ‰, which has been suggested as appropriate for soils dominated by C4 plants in northern China (Currell et al., 2010). This model yields some relatively low q values (0.59 for G06-15 and 0.65 for G08-15), possibly because several unaccounted-for factors contribute to variable δ13CRECH values, e.g., local methanogenesis and pH or temperatures in the soil zones.

For the CMB model, q=mDICrech/mDICfinal, where mDICrech is the DIC molar concentration in the recharge water and mDICfinal is the DIC molar concentration in the final groundwater. mDICfinal was calculated using mDICfinal=mDICrech+[mCa+Mg-SO4+0.5(Na+K-Cl)] (Fontes and Garnier, 1979). DICrech was mainly HCO3 in the recharge water when pH values were between 6.4 and 10.3, and the carbonate equilibrium constant varied with temperature (Clark and Fritz, 1997). mDICrech was calculated from the estimated pH and temperature conditions in the recharge environment; e.g., at pH = 6 and T= 15 C, mDICrech= 10 mmol L−1 (Currell et al., 2010).

The corrected radiocarbon ages are shown in Table S2. The residence time of deep groundwater ranges from 15 959 to 39 050 a BP, which is significantly longer than that of the groundwater in the shallow aquifer (9510 a BP to modern). Moreover, most brackish and fresh groundwater ages are modern, while brine has a longer residence period (5590–1245 a BP) and a broad range of residence times.

Figure 514C activity with sampling depth in groundwater.


5 Discussion

5.1 Isotopic analysis of groundwater origin and recharge

Deuterium and oxygen-18 are good tracers of groundwater origin and climatic conditions during recharge periods (Clark and Fritz, 1997). When combined with groundwater residence times, they can also identify modern recharge and paleo-recharge (Han et al., 2014).

The depletion of 18O and 2H values in the deep fresh groundwater (Fig. 4) can be attributed to a cold climate (Kreuzer et al., 2009) and the residence times (33 951–39 050 a BP) of the P15-150 and P14-300 samples, which may suggest that there was a recharge during the last glacial maximum. The stable isotopes in P16-100 are heavier, reflecting a recharge history linked to the warm climate in the previous deglaciation (Hendry and Wassenaar, 2000). The stable isotope values of river samples are similar to those of the shallow brackish and fresh groundwater samples, which are approximately modern, indicating lateral recharge of surface water locally. Meanwhile, in Fig. 4, G03-5 is close to the rainfall sample, indicating that modern precipitation is a new recharge source. The trend towards δ2H and δ18O enrichment in brine and saline groundwater can be attributed to seawater infiltration during the Holocene transgression period, which has been confirmed by other studies of the Bohai Sea coast (Li et al., 2017; Du et al., 2016). Additionally, due to the mixing of meteoric water and the subsequent nonequilibrium fractionation of hydrogen isotopes during evaporation (Clark and Fritz, 1997), the CSW sample is characterized by 18O enrichment but 2H depletion compared to seawater.

5.2 Hydrochemical analysis of sources of salinity

To distinguish the sources of groundwater salinity, the PHREEQC code (Parkhurst and Appelo, 2013) was used to measure and plot the theoretical seawater–freshwater mixing line (“mixing line”) and seawater evaporation line (“evaporation line”) using hydrogeochemical modeling. The two simulation effects serve as references for groundwater hydrochemical characteristics (Figs. 6 and 7). For the Na-Cl (Fig. 6a), Mg-Cl (Fig. 6b), and Br-Cl (Fig. 7a) diagrams, the measured brackish, saline, and brine groundwater samples fit quite well to modeling mixing lines, and evaporation lines follow a linear trend from the least to the most saline. This strongly demonstrates that the salt in these water samples is mainly of marine origin. The major ion concentrations in some samples (such as brine) are higher than those in seawater, suggesting that the enriched ions are associated with evaporation processes rather than seawater intrusion (Colombani et al, 2017).

Moreover, the samples deviate from the modeling lines (Fig. 6c and d), indicating that there may be other hydrogeochemical processes that are responsible for the modified ionic compositions (Giambastiani et al., 2013). (1) Ca2+ depletion is seen for the P18 and P12 samples in Fig. 6d. This phenomenon is likely explained by gypsum (CaSO4) precipitation. The evaporation line reveals that the Ca2+ composition of evaporating seawater follows a hooked trajectory (Fig. 6d). During evaporation to the point of gypsum saturation, the Ca2+ concentration of the residual CSW progressively decreases. (2) Ca2+ and SO42- are present in excess in most of the fresh and brackish samples (Fig. 6c and d), which can be attributed to mineral dissolution along with stream water recharging (such as gypsum dissolution), highlighting the occurrence of some degree of dilution with continental runoff since the Holocene regression. (3) Decomposition of the organic matter abundant in marine or lagoon facies sediments can result in the release of bromide ions, thus pushing the Br / Cl ratios of saline groundwater samples higher than the mixing line (Fig. 7b).

Figure 6Hydrochemical relationships between Cl and major ions of the measured samples and simulated results in groundwater. Seawater–freshwater mixing line: theoretical mixing between seawater and deep fresh groundwater. Blue numbers: mixing ratios of seawater. Seawater evaporation line: theoretical evaporation of Bohai seawater. Red numbers: evaporation rates. G: precipitation of gypsum. H: precipitation of halite. Dissolution: possible gypsum dissolution along with river water recharging. The symbols used for the samples are the same as in Fig. 6.


Figure 7Relationship between the chloride and bromide contents in water samples. Symbols used are the same as in Fig. 6.


5.3 Mixing processes

Figure 8 depicts the relationship between δ18O and Cl in different water samples. There are higher Cl concentrations and lower δ18O values in brine samples than in seawater, meaning that simple two-end-member mixing cannot adequately explain the groundwater salinization. Stable isotopes of high-TDS saline and brine samples fall between the seawater and CWS mixing lines, suggesting a potential for three-end-member mixing processes (Douglas et al., 2000). Therefore, we considered SW01 (seawater) and P18 (mostly saline, but with relatively depleted stable isotopes) as two saline end-members. P16-100, which was most likely recharged during the last deglaciation, was chosen to represent fresh end-members that could have been impacted by overlying seawater or CSW during the Holocene transgression. In Fig. 8, an inferred salinization zone is established that includes almost all saline and brine groundwater samples, demonstrating the salinization processes in which fresh groundwater mixed with either seawater, CSW, or a mixture of both.

The fresh and brackish groundwater samples, on the other hand, have low Cl concentrations and less 18O, so they deviate from the assumed salinization zone and approach the river samples in Fig. 8, implying a river water–groundwater mixing trend. LH02 (with low δ18O) and SH02 (with high δ18O) were selected to represent river water end-members for different river recharge of groundwater in the study area, while G09-15 (saline but with river-like stable isotope levels) was considered as a groundwater end-member. It is presumed that a freshening zone could form between the two river water-groundwater mixing lines, indicating the occurrence of freshening processes, which would be in agreement with the continental runoff dilution discussed in Sect. 5.2.

Figure 8Plot showing the relationship between Cl and δ18O in different water samples as a means to probe various mixing processes in the Luanhe River delta. The symbols are the same as used in Fig. 6. The green area is the assumed freshening zone, and the purple area is the assumed salinization zone.


6 Interpretation of paleo-environmental development

Based on an analysis of a range of evidence related to Quaternary geographic evolution, it is possible to understand the changes in hydrogeological conditions in the past (van Engelen et al., 2018). Other authors (Wang et al., 1981; Peng et al., 1981; Li et al., 1982, 1983; Xu et al., 2018) found that the Pleistocene transgression events – related to marine isotope stage (MIS) 3 and 5 – reached the study area once, which would have resulted in groundwater salinization. Since the last deglaciation (about 15 ka BP), the paleo-coastline has dropped to a depth of approximately 100 m below the present sea level along the shelf edge (Li et al., 2014). Stronger river downcutting and flushing in the study region would have promoted a large fresh recharge of groundwater. For example, P16-100 (freshwater) was sampled from a relatively deep position (100 m below the surface) and has an estimated groundwater age of between 15 959 to 17 490 a BP, which likely provides evidence that the salinization groundwater related to MIS 5 and/or 3 marine transgression could have been flushed out until the latest Late Pleistocene. Accordingly, we believe that the observed saline groundwater in the Luanhe River delta is probably related to the subsequent Holocene marine transgression. This research illuminates the evolutionary pattern of saline groundwater, as shown in Table 1 and Fig. 9, where three phases are synthesized and reconstructed.

Table 1Saline groundwater evolution processes in the study area.

Download Print Version

The global sea level was affected by deglaciation of the ice sheet (Fairbanks, 1989), causing the sea level to rise rapidly during the deglaciation period (15.4–7 ka BP) (Li et al., 2014). The Holocene transgression stage, which occurred between 9 and 7 ka BP, resulted in the study area being inundated by seawater (Xu et al., 2015; Xue, 2009, 2014) (Fig. 9a). On the one hand, there would have been a tendency for the denser seawater to infiltrate through the aeration zone (Santucci et al., 2016); on the other hand, the sea-level rise would have caused the seawater–freshwater interface to move landward (Ferguson and Gleeson, 2012), both of which contributed to paleo-seawater intrusion. G08-40 contains 27.173 g L−1 of TDS, which is similar to that of SW01. Simultaneously, the residence time (9810–6884 a BP) indicates the existence of paleo-seawater trapped by low-permeability aquitard sediments, which may be another critical salinity source for neighboring aquifers in the coastal zone (Post and Kooi, 2003; Lee et al., 2016).

Figure 9Diagram of paleo-environmental development during the Holocene and evolutionary pattern of saline groundwater.


The presence of paleo-seawater intrusion during the Quaternary has been recorded in other coastal regions worldwide (Groen et al., 2000; Bouchaou et al., 2009; Wang and Jiao, 2012; Delsman et al., 2014; Tran et al., 2020; Han et al., 2020). According to the works described above, the salinity of the groundwater after salinization could not exceed that of seawater due to paleo-seawater intrusion.

Other salinization processes that occurred during paleo-environmental growth are likely to be correlated with this brine groundwater.

6.1 Phase 2: highstand system tract–old Luanhe River delta development (7–3.5 ka BP)

The good fit between the measured hydrochemistry and simulated evaporation lines (Figs. 6 and 7) is an indicator that the brine samples were associated with the seawater that was exposed to evaporation during geological history. Previous research has revealed that the lagoon was active during the progradation of the old Luanhe River delta between 7 and 3.5 ka BP (He et al., 2020; Xu et al., 2020). Meanwhile, the region has had a relatively arid climate since 5500 a BP, which may have led to increased evaporation (Jin, 1984). The ancient lagoon would have been an ideal location for evaporating seawater that had been trapped due to storms or tides (Fig. 9b). As a result, concentrated saline water (CSW) with a salinity greater than that of seawater would have been created, and the CSW would have undergone two processes: (1) it would have infiltrated and descended to the lower part of the aquifer due to its higher density, before combining with the salinized groundwater from phase 1, resulting in a three-end-member mixing scenario in the relationship diagram (Fig. 8); (2) after reaching saturation during the later stages of evaporation, precipitation of minerals such as gypsum, calcite, and halite would have occurred, which would have been redissolved by meteoric waters or seawater, resulting in high-salinity water which would then have been subjected to the above process. The Br / Cl ratios in certain fresh or brine groundwater samples deviate from the evaporation line (Fig. 7b), which may be related to halite precipitation and redissolution. These two processes caused the groundwater salinity to rise even further, resulting in the formation of brine groundwater with 3 times the TDS of seawater, such as G03-20, with a residence time range of 4323–5590 a BP.

6.2 Phase 3: new Luanhe River delta development (3.5 ka BP to present)

Since about 3500 a BP, a nearly 90 diversion of the Luanhe River channel in the study area has resulted in new delta development (Wang et al., 2007; Xue, 2016). There are some signs of a lagoon environment in the new Luanhe River delta (Cheng et al., 2020), and, as previously discussed, the brine groundwater sample G10-30 can be attributed to evaporation in a lagoon setting (Fig. 9c). However, some factors are likely to limit CSW formation in the study area: (1) the relatively low evaporation capacity due to the semi-humid climate since about 2.5 ka BP (Jin, 1984), (2) diluvial deposits or artificial reclamation would have filled the coastal lowlands such as lagoons, and (3) offshore levees prevent the seawater from flooding inland during storms or tides. Unlike the old Luanhe River delta, these factors may also explain why the current Luanhe River delta does not have high-TDS brine groundwater.

In addition, the brackish and low-TDS saline groundwater with a relatively modern age (e.g., G09-15) and the river-like stable isotope levels (Figs. 4 and 8) are compelling evidence that freshening processes have occurred in the delta plain. With the semi-humid paleoclimate, some abandoned channels have developed into small rivers such as the Suhe River and Shahe River following the diversion of the ancient Luanhe River (Gao, 1981). Firstly, the lateral recharge from the surface stream plays a role in washing out the salty groundwater. Secondly, due to the unsuitability of saline groundwater throughout human history, river irrigation has been commonly used for agricultural activities in the study region, which has freshened the upper saline aquifer (Fig. 9c). Some groundwater samples found above the seawater mixing line in the Ca-Cl and SO4-Cl relationship diagrams (Fig. 6c, d) may be related to mineral dissolution during river water or irrigation recharge. However, saline groundwater can be washed out over time in coastal zones with low-permeability marine layers and a low hydraulic gradient (van Engelen et al., 2019; Han et al., 2020).

In summary, the evolution of saline groundwater in the study area results from paleo-environment development such as sea-level changes, paleogeography, and paleoclimate, and is significantly affected by human activities. The coastal brine groundwater found on the Bohai Sea coast, such as Bohai Bay (Li et al., 2017) and Laizhou Bay (Han et al., 2014), is a special product of geological evolution. The change in sea level during the Late Pleistocene would have favored marine intrusion and a similar sedimentary environment on the Bohai coast, allowing this study to infer the following conditions for brine formation in this region: (1) stable evaporative environments (e.g., a lagoon), (2) suitable climatic conditions (e.g., aridity), (3) entry of seawater into evaporative environments (e.g., via storms or tides), and (4) a long time period for salinity accumulation.

7 Conclusions

In this study, we used a range of isotopic–geochemical methods to analyze the groundwater recharge and salinity sources in the Luanhe River delta. The isotopic results (2H, 18O, 14C) show that deep confined groundwater was recharged during the Late Pleistocene cold period, shallow saline and brine groundwater was recharged during the warm Holocene period, and shallow brackish and fresh groundwater was mainly recharged by surface water. The hydrogeochemical modeling (PHREEQC) results show that seawater or evaporated seawater is the primary source of the salt in salinized groundwater. Variation among multiple water samples in the 18O–Cl relationship further indicates multiple-end-member mixing, which is useful for assessing the salinization and/or freshening processes in the aquifers. Our study shows that multiple water types are particularly associated with complex geographic evolution in coastal areas. Rising sea levels caused lowland coastal areas to be inundated by seawater, inducing paleo-seawater intrusion. Coastal deltas developed after a significant drop in sea levels. The concentration of saline water in the lagoon environment at the front of the delta continuously added salinity to the groundwater. Thus, brine groundwater was formed under the effects of evaporation, mixing, and dissolution. In contrast, the lateral recharge by surface water and irrigation return would have caused the slow washing out of salinized groundwater in the delta plain.

Given that most coastal zones around the world have experienced transgression/regression events in the Quaternary Period, this work's findings will promote a better understanding of the origin of salinization in coastal aquifers. In addition, it is important to recognize the potential for connate saline groundwater previously preserved in adjunct aquifers to leak due to the overextraction of deep groundwater. To effectively prevent pollution from saline groundwater movement, this study recommends extensive characterization of the groundwater interface dynamics, such as the fresh/saline, fresh/brine, and brine/seawater interfaces, and that continuous monitoring of water quality and levels across the aquifers should be maintained.

Data availability

The data used in this paper are available in the Supplement.


The supplement related to this article is available online at:

Author contributions

XD and MG carried out the study and wrote the manuscript with support from ZW. MG and GH helped supervise the project. QS and GH helped collect the data. ZW, DA, and HJ helped improve the English language of the manuscript. All co-authors contributed to reviewing and editing the manuscript.

Competing interests

The contact author has declared that neither they nor their co-authors have any competing interests.


Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.


This study was financially supported by the National Natural Science Foundation of China (U2106203, 41977173), National Geological Survey Project of China Geology Survey (no. DD20211401). The authors would like to thank Sen Liu, Chenxin Feng, Chen Sheng, Xueyong Huang, and Haihai Zhuang for their help and support in collecting the field data and conducting the geological survey.

Financial support

This research has been supported by the National Natural Science Foundation of China (grant nos. 41977173 and U2106203) and the National Geological Survey Project of China Geology Survey (no. DD20211401).

Review statement

This paper was edited by Nadia Ursino and reviewed by three anonymous referees.


Akouvi, A., Dray, M., Violette, S., de Marsily, G., and Zuppi, G. M.: The sedimentary coastal basin of Togo: example of a multilayered aquifer still influenced by a palaeo-seawater intrusion, Hydrogeol. J., 16, 419–436,, 2008. 

Aquilina, L., Vergnaud-Ayraud, V., Les Landes, A. A., Pauwels, H., Davy, P., Petelet-Giraud, E., Labasque, T., Roques, C., Chatton, E., and Bour, O.: Impact of climate changes during the last 5 million years on groundwater in basement aquifers, Sci. Rep.-UK, 5, 14132,, 2015. 

Bouchaou, L., Michelot, J. L., Qurtobi, M., Zine, N., Gaye, C. B., Aggarwal, P. K., Marah, H., Zerouali, A., Taleb, H., and Vengosh, A.: Origin and residence time of groundwater in the Tadla basin (Morocco) using multiple isotopic and geochemical tools, J. Hydrol., 379, 323–338,, 2009. 

Cartwright, I., Currell, M., Cendon, D., and Meredith, K.: A review of the use of radiocarbon to estimate groundwater residence times in semi-arid and arid areas, J. Hydrol., 580, 124247,, 2020. 

Cary, L., Petelet-Giraud E., Bertrand, G., Kloppmann, W., Aquilina, L., Martins, V., Hirata, R., Montenegro, S., Pauwels, H., Chatton, E., Franzen, M., and Aurouet, A.: Origins and processes of groundwater salinization in the urban coastal aquifers of Recife (Pernambuco, Brazil): A multi-isotope approach, Sci. Total Environ., 530–531, 411–429,, 2015. 

Chen, Z. Y., Qi, J. X., Xu, J. M., Xu, J. M., Ye, H., and Nan, Y. J.: Paleoclimatic interpretation of the past 30 ka from isotopic studies of the deep confined aquifer of the North China Plain, Appl. Geochem., 18, 997–1009,, 2003. 

Cheng, L. Y., Xu, Q. M., Guo, H., Li, M., Yang, N., Liu, J. B., Zhao, J. J., and Guo, J. J.: The Late Holocene Stratum and evolution in the Luanhe River Delta, Quaternary Sciences, 40, 751-763,, 2020 (in Chinese with English abstract). 

Clark, I. D. and Fritz, P. (Eds): Environmental Isotopes in Hydrogeology, Lewis Publishers, New York,, 1997. 

Colombani, N., Cuoco, E., and Mastrocicco, M.: Origin and pattern of salinization in the Holocene aquifer of the southern Po Delta (NE Italy), J. Geochem. Explor., 175, 130–137,, 2017. 

Costall, A. R., Harris, B. D., Teo, B., Schaa, R., Wagner, F. M., and Pigois, J. P.: Groundwater Throughflow and Seawater intrusion in High Quality coastal Aquifers, Sci. Rep.-UK, 10, 9866,, 2020. 

Craig, H.: Standard for reporting concentration of deuterium and oxygen-18 in natural water, Science, 133, 1833–1834,, 1961. 

Currell, M. J., Cartwright, I., Bradley, D. C., and Han, D. M.: Recharge history and controls on groundwater quality in the Yuncheng Basin, north China, J. Hydrol., 385, 216–229,, 2010. 

Dang, X. Z., Gao, M. S., Wen, Z., Jakada, H., Hou, G. H., and Liu, S.: Evolutionary process of saline groundwater influenced by palaeo-seawater trapped in coastal deltas: A case study in Luanhe River Delta, China, Estuar. Coast. Shelf S., 244, 106894,, 2020. 

de Montety, V., Radakovitch, O., Vallet-Coulomb, C., Blavoux, B., Hermitte, D., and Valles, V.: Origin of groundwater salinity and hydrogeochemical processes in a confined coastal aquifer: case of the Rhone delta (Southern France), Appl. Geochem., 23, 2337–2349,, 2008. 

Delsman, J. R., Hu-a-ng, K. R. M., Vos, P. C., de Louw, P. G. B., Oude Essink, G. H. P., Stuyfzand, P. J., and Bierkens, M. F. P.: Paleo-modeling of coastal saltwater intrusion during the Holocene: an application to the Netherlands, Hydrol. Earth Syst. Sci., 18, 3891–3905,, 2014. 

Douglas, M., Clark, I. D., Raven, K., and Bottomley, D.: Groundwater mixing dynamics at a Canadian Shield mine, J. Hydrol., 235, 88–103,, 2000. 

Du, Y., Ma, T., Chen, L. Z., Shan, H. M., Xiao, C., Lu, Y., Liu, C. F., and Cai, H. S.: Genesis of salinized groundwater in Quaternary aquifer system of coastal plain, Laizhou Bay, China: Geochemical evidences, especially from bromine stable isotope, Appl. Geochem., 59, 155–165,, 2015. 

Du, Y., Ma, T., Chen, L., Xiao, C., and Liu, C. F.: Chlorine isotopic constraint on contrastive genesis of representative coastal and inland shallow brine in China, J. Geochem. Explor., 170, 21–29,, 2016. 

Edmunds, W. M.: Palaeowaters in European coastal aquifers-the goals and main conclusions of the PALAEAUX project, Geological Society London Special Publications, 189, 1–16,, 2001. 

Fairbanks, R. G.: A 17,000-year glacio-eustatic sea level record: influence of glacial melting rates on the Younger Dryas event and deep ocean circulation, Nature, 342, 637–647,, 1989. 

Feng, J. and Zhang, W.: The evolution of the modern Luanhe River delta, north China, Geomorphology, 25, 269–278,, 1998. 

Ferguson, G. and Gleeson, T.: Vulnerability of coastal aquifers to groundwater use and climate change, Nat. Clim. Change, 2, 342–345,, 2012. 

Fontes, J. C. and Garnier, J. M.: Determination of the initial 14C activity of the total dissolved carbon: a review of the existing models and a new approach, Water Resour. Res., 15, 399–413,, 1979. 

Gao, S. M.: Facies and sedimentary model of the Luan River delta, Acta Geographica Sinica, 48, 303–314, 1981 (in Chinese with English abstract). 

Giambastiani, B. M. S., Colombani, N., Mastrocicco, M., and Fidelibus, M. D.: Characterization of the lowland coastal aquifer of Comacchio (Ferrara, Italy): Hydrology, hydrochemistry and evolution of the system, J. Hydrol., 501, 35–44,, 2013. 

Gibson, J. J., Edwards, T. W., Bursey, G. G., and Prowse, T. D.: Estimating evaporation using stable isotopes: quantitative results and sensitivity analysis for two catchments in Northern Canada, Nord. Hydrol., 24, 79–94,, 1993. 

Groen, J., Velstra, J., and Meesters, A.: Salinization processes in paleowaters in coastal sediments of Suriname: evidence from δ37Cl analysis and diffusion modelling, J. Hydrol., 234, 1–20,, 2000. 

Han, D. M., Kohfahl, C., Song, X. F., Xiao, G. Q., and Yang, J. L.: Geochemical and isotopic evidence for Palaeo-Seawater intrusion into the south coast aquifer of Laizhou Bay, China, Appl. Geochem., 26, 863–883,, 2011. 

Han, D. M., Song, X. F., Currell, M. J., Yang, J. L., and Xiao, G. Q.: Chemical and isotopic constraints on the evolution of groundwater salinization in the coastal plain aquifer of Laizhou Bay,China, J. Hydrol., 508, 12–27,, 2014. 

Han, D. and Currell, M. J.: Delineating multiple salinization processes in a coastal plain aquifer, northern China: hydrochemical and isotopic evidence, Hydrol. Earth Syst. Sci., 22, 3473–3491,, 2018. 

Han, D. M., Cao, G. L., Currell, M. J., Priestley, S. C., and Love, A. J.: Groundwater salinization and flushing during glacial-interglacial cycles: insights from aquitard porewater tracer profiles in the North China Plain, China, Water Resour. Res., 56, e2020WR027879,, 2020. 

He, L., Amorosi, A., Ye, S. Y., Xue, C. T., Yang, S. X., and Laws, E. A.: River avulsions and sedimentary evolution of the Luanhe fan-delta system (North China) since the late Pleistocene, Mar. Geol., 425, 106194,, 2020. 

Hendry, M. J. and Wassenaar, L. I.: Controls on the distribution of major ions in pore waters of thick surficial aquitard, Water Resour. Res., 36, 503–513,, 2000. 

IAEA/WMO: Global Network of Isotopes in Precipitation, The GNIP Database, Vienna, (last access: 18 November 2013), 2006. 

Jayathunga, K., Diyabalanage, S., Frank, A. H., Chandrajith, R., and Barth, J. A. C.: Influences of seawater intrusion and anthropogenic activities on shallow coastal aquifers in Sri Lanka: evidence from hydrogeochemical and stable isotope data, Environ. Sci. Pollut. R., 27, 23002–23014,, 2020. 

Jiao, J. J. and Post, V. (Eds): Coastal Hydrology, Cambridge University Press, New York,, 2019. 

Jin, X. F.: The spore-pollen assemblages and the stratigraphy and palaeogeography in western Bohai Sea since late Pleistocene, Marine Science Bullition 3, 16–24, 1984 (in Chinese with English abstract). 

Kooi, H., Groen, J., and Leijnse, A.: Modes of seawater intrusion during transgressions, Water Resour. Res., 36, 3581–3589,, 2000. 

Kreuzer, A. M., Rohden, C. V., Friedrich, R., Chen, Z. Y., Shi, J. S., Hajdas, I., Kipfer, R., and Aeschbach-Hertig, W.: A record of temperature and monsoon intensity over the past 40 kyr from groundwater in the North China Plain, Chem. Geol., 259, 168–180,, 2009. 

Larsen, F., Tran, L. V., Van Hoang, H., Tran, L. T., Christiansen, A. V., and Phan, N. Q.: Groundwater salinity influenced by Holocene seawater trapped in incised valleys in the Red River delta, Nat. Geosci., 10, 376–381,, 2017. 

Lee, S., Currell, M., and Cendon, D. I.: Marine water from mid-Holocene sea level highstand trapped in a coastal aquifer: Evidence from groundwater isotopes, and environmental significance, Sci. Total Environ., 544, 995–1007,, 2016. 

Li, G. X., Li, P., Liu, Y., Qiao, L. L., Ma, Y. Y., Xu, J. S., and Yang, Z. G.: Sedimentary system response to the global sea level change in the East China Seas since the last glacial maximum, Earth-Sci. Rev., 139, 390–405,, 2014. 

Li, H. M. and Wang, J. D.: Palaeomagnetic study on drill core from northern Bohai coastal plain, Geochimica, 2, 196–204, 1983 (in Chinese with English abstract). 

Li, J., Liang, X., Jin, M. G., and Mao, X. M.: Geochemical signature of aquitard pore water and its paleo-environment implications in Caofeidian Harbor, China, Geochem. J., 47, 37–50,, 2013. 

Li, J., Liang, X., Jin, M. G., Yang, J. L., Ma, B., Ge, Q.: Origin and Evolution of Aquitard Porewater in the Western Coastal Plain of Bohai Bay, China, Groundwater, 55, 917—925,, 2017. 

Li, Y. F., Gao, S. M., and An, F. T.: A preliminary study of the Quaternary marine strata and its paleogeographic significance in the Luanhe delta region, Oceanologia et Limnologia Sinica, 13, 433–439, 1982 (in Chinese with English abstract). 

Liu, S., Tang, Z. H., Gao, M. S., and Hou, G.H.: Evolutionary process of saline-water intrusion in Holocene and Late Pleistocene groundwater in southern Laizhou Bay, Sci. Total Environ., 607–608, 586–599,, 2017. 

Ma, F. S., Wei, A. H., Deng, Q. H., and Zhao H. J.: Hydrochemical Characteristics and the Suitability of Groundwater in the Coastal Region of Tangshan, China, J. Earth Sci.-China, 26, 1067–1075,, 2014. 

Niu, Z. X., Jiang, X. W., and Hu, Y. Z.: Characteristics and causes of hydrochemical evolution of deep groundwater in the Luanhe delta, Hydrogeology and Engineering Geology, 46, 27–34, 2019 (in Chinese with English abstract). 

Parkhurst, D. L. and Appelo, C. A. J.: Description of Input and Examples for PHREEQC Version 3 – A Computer Program for Speciation, Batch-Reaction, One-Dimensional Transport, and Inverse Geochemical Calculations, in: Techniques and Methods, book 6, chap. A43, US Geological Survey, Denver, p. 497, (last access: October 2019), 2013. 

Peng, G., Jiao, W. Q., Li, D. M., and Li, G. Y.: Division and correlation of the late Quaternary stratigraphy and discussion on the recent tectonic movement in the region of the Luanhe River Delta, Seismoloqy and Geology, 3, 31–36, 1981 (in Chinese with English abstract). 

Pearson, F. J. and Hanshaw, B. B.: Sources of dissolved carbonate species in groundwater and their effects on carbon-14 dating, in: IAEA, Isotope Hydrology, IAEA, Vienna, 271–285, 1970. 

Post, V. E. and Kooi, H.: Rates of salinization by free convection in high-permeability sediments: insight from numerical modeling and application to Dutch coastal area, Hydrogeol. J., 11, 549–559,, 2003. 

Qi, H. H., Ma, C. M., He, Z. K., Hu, X. J., and Gao, L.: Lithium and its isotopes as tracers of groundwater salinization: A study in the southern coastal plain of Laizhou Bay, China, Sci. Total Environ., 650, 878–890,, 2019. 

Reilly, T. E. and Goodman, A. S.: Quantitative analysis of saltwater-freshwater relationships in groundwater systems-a historical perspective, J. Hydrol., 80, 125–160,, 1985. 

Sanford, W. E.: Groundwater hydrology: Coastal flow, Nat. Geosci., 3, 671–672,, 2010. 

Santucci, L., Carol, E., and Kruse E.: Identification of palaeo-seawater intrusion in groundwater using minor ions in a semi-confined aquifer of the Río de la Plata littoral (Argentina), Sci. Total Environ., 566–567, 1640–1648,, 2016. 

Sola, F., Vallejos, A., Daniele, L., and Pulido-Bosch, A.: Identification of a Holocene aquifer-lagoon system using hydrogeochemical data, Quaternary Res., 82, 121–131,, 2014. 

Stumpp, C., Ekdal, A., Gönenc, I. E., and Maloszewski, P.: Hydrological dynamics of water sources in a Mediterranean lagoon, Hydrol. Earth Syst. Sci., 18, 4825–4837,, 2014. 

Tulipano, L., Fidelibus, M. D., and Panagopoulos, A.: COST ACTION 621 Final Report, Groundwater Management of Coastal Karstic Aquifers, Office for the Official Publications of European Communitied, Luxembourg, Vol. II, 366 pp., ISBN 92-894-0015-1, 2005. 

Tran, D. A., Tsujimura, M., Vo L. P., Nguyen, V. T., Kambuku, D., and Dang, T. D.: Hydrogeochemical characteristics of a multi-layered coastal aquifer system in the Mekong Delta, Vietnam, Environ. Geochem. Hlth., 42, 661–680,, 2020. 

UN Atlas: 44 Percent of us Live in Coastal Areas, (last access: 12 June 2013), 2010. 

Vallejos, A., Sola, F., Yechieli, Y., and Pulido-Bosch, A.: Influence of the paleogeographic evolution on the groundwater salinity in a coastal aquifer. Cabo de Gata aquifer, SE Spain, J. Hydrol., 557, 55–66,, 2018. 

van Engelen, J., Oude Essink, G. H. P., Kooi, H., and Bierkens, M. F. P.: On the origins of hypersaline groundwater in the Nile Delta aquifer, J. Hydrol., 560, 301–317,, 2018. 

van Engelen, J., Verkaik, J., King, J., Nofal, E. R., Bierkens, M. F. P., and Oude Essink, G. H. P.: A three-dimensional palaeohydrogeological reconstruction of the groundwater salinity distribution in the Nile Delta Aquifer, Hydrol. Earth Syst. Sci., 23, 5175–5198,, 2019. 

Wang, P. X., Min, Q. B., Bian, Y. H., and Cheng, X. R.: Strata of Quaternary transgressions in east China: A preliminary study, Acta Geol. Sin., 1981, 1–13, 1981 (in Chinese with English abstract). 

Wang, Y. and Jiao, J. J.: Origin of groundwater salinity and hydrogeochemical processes in the confined Quaternary aquifer of the Pearl River Delta, China, J. Hydrol., 438–439, 112–124,, 2012. 

Wang, Y., Fu, G., and Zhang, Y.: River-sea interactive sedimentation and plain morphological evolution, Quaternary Science, 27, 674–689,, 2007 (in Chinese with English abstract). 

Werner, A. D.: A review of seawater intrusion and its managementin Australia, Hydrogeol. J., 18, 281–285,, 2010. 

Werner, A. D., Bakker, M., Post, V. E. A., Vandenbohede, A., Lu, C. H., Ataie-Ashtiani, B., Simmons, C. T., and Barry, D. A.: Seawater intrusion processes, investigation and management: Recent advances and future challenges, Adv. Water Resour., 51, 3–26,, 2013. 

Xu, Q. M., Yuan, G. B., Zhang, J. Q., and Qin, Y. F.: Stratigraphic division of the Late Quaternary strata along the coast of Bohai bay and its geology significance, Acta Geol. Sin., 85, 1352–1367, 2011 (in Chinese with English abstract). 

Xu, Q. M., Yang, J. L., Yuan, G. B., Chu, Z. X., and Zhang, Z. K.: Stratigraphic sequence and episodes of the ancient Huanghe Delta along the southwestern Bohai Bay since the LGM, Mar. Geol., 367, 69–82,, 2015. 

Xu, Q. M., Yang, J. L., Hu, Y. Z., Yuan, G. B., and Deng, C. L.: Magnetostratigraphy of two deep boreholes in the southwestern BohaiBay: its tectonic implications and constraints on ages of volcanic layers, Quat. Geochronol., 43, 102–114,, 2018. 

Xu, Q. M., Meng, L. S., Yuan, G. B., Teng, F., Xin, H. T., and Sun, X. M.: Transgressive wave-and tide-dominated barrier-lagoon system and sea-level rise since 8.2 ka recorded in sediments in northern Bohai Bay, China, Geomorphology, 352, 106978,, 2020. 

Xue, C. T.: Historical changes of coastlines on west and south coasts of Bohai Sea since 7000 a B.P., Scientia Geographic Sinica, 29, 217–222,, 2009 (in Chinese with English abstract). 

Xue, C. T.: Missing evidence for stepwise postglacial sea level rise and an approach to more precise determination of former sea levels on East China Sea Shelf. Mar. Geol., 348, 52–62,, 2014.  

Xue, C. T.: Extents, type and evolution of Luanhe River fan-delta system, China, Marine Geology & Quaternary Geology, 36, 13–22, 2016 (in Chinese with English abstract). 

Zhou, X.: Basic characteristics and resource classification of subsurface brines in deep-seated aquifers, Hydrogeology & Engineering Geology, 40, 4–10, 2013 (in Chinese with English abstract). 

Short summary
Due to the widespread salinity in coastal zones, freshwater supplies are now under pressure. This study analyzes the recharge and salinity sources of coastal groundwater (brine, saline, brackish, fresh), as well as salinization and freshening processes, using hydrochemical and isotopic methods. The evolution of saline groundwater is studied using sedimentary characteristics as multiple lines of evidence. This study can be used to better understand groundwater evolution in the coastal zone.