The hydrogeochemical and stable isotope compositions of aquitard porewater samples from three boreholes were investigated to determine the origin, salinization, and hydrochemical evolution of water in the North Jiangsu coastal plain, China. Three porewater groups were identified based on the water-bearing subsystems. The total dissolved solids (TDS) of porewater samples highly varied from 0.03 to 26.1 g/L. Molar Cl/Br ratios and δ18O and δ2H data indicate that the source of Group 1 salinized porewater was the Holocene seawater, whereas Group 3 salinized samples were probably related to the remnant palaeoseawater of the Late Pleistocene. Group 2 samples had low salinity (TDS < 1 g/L) and undetectable Br− concentrations, which were probably recharged during a colder period and without evidence of seawater involvement. The salinized porewater was likely diluted by freshwater, as evidenced by its depleted isotopes and low salinity relative to standard seawater. The ionic ratios and ionic deltas indicate that considerable water-rock interactions (e.g., cation exchange and mineral weathering) also accounted for the hydrochemical constituents of porewaters, and cation exchange seems to be more noticeable at low salinity.
National Natural Science Foundation of China41502231412722581. Introduction
Coastal areas are among the most densely populated areas in the world; it is estimated that half of the world’s population lives in coastal areas [1]. In China, the continental coastline exceeds 18,000 km, and more than 60% of the population lives in 14 coastal provinces [2]. With rising population levels and growing industrial development, vulnerable water resources seem to be overused, especially in these coastal areas [3, 4].
Water management problems can be exacerbated by the presence of saline groundwater. Groundwater salinization is likely associated with modern seawater intrusion due to overpumping [5, 6], contributions of deep saline waters (fossil seawater) [7–9], and water-rock interaction with evaporates [10, 11]. Research suggests that released porewater from the adjacent aquitards is a substantial proportion of the amount of aquifer groundwater that can be exploited [12]. Aquitards act as natural protective barriers for underlying aquifers and also contain large amounts of salt and contaminants that may leak into aquifers [13, 14]. In the context of potential recharge sources for aquifers, the successful description of porewater constituents and hydrochemical processes involved in coastal aquitards will benefit the sustainable exploitation, evaluation, and protection of groundwater resources in aquifers.
Hydrochemical and stable isotopes of groundwater in aquifers have classically been considered tools to deduce the sources of dissolved salts [8, 15, 16]. By contrast, few studies have examined porewater in aquitards even though aquitards commonly have a greater capability to contain water with high salinity relative to aquifers [17, 18]. Diffusion is considered the dominant mechanism of the transport of solutes in the coastal aquitards [19]. Porewater salinity in aquitards cannot simply relate to seawater intrusion and may involve complex geochemical processes (such as water-rock interaction) because of the extremely long residence time of water in aquitards [20, 21]. Further, a combination of past and recent processes, possibly including several phases of regression, transgression, formation of salt water, and mixing, can result in salinization of porewater [8]. In China, the sources of saline porewater in the coastal plain aquitards are still not clear.
In North Jiangsu coastal plain, the hydrogeological conditions have been studied comprehensively [4, 22, 23]. Most of the saline and brackish groundwater samples in the aquifers were found to be derived from relict seawater [4], because three marine invasions and regressions have occurred along the coastal area of eastern China since the Late Pleistocene [24]. As a result of the interest in the hydrochemical processes of porewaters in the clay-rich deposits, major ions and stable isotopes (δ18O, δ2H) are used in the present study to provide insight into the salinity sources, mixing behaviors, and hydrochemical evolution. The objectives of this paper are (1) to analyze the vertical distribution of hydrochemistry and stable isotopic signatures of porewater in three cores (SY1, SY2, and SY3) in North Jiangsu coastal plain, (2) to ascertain the primary salinity sources of porewater and the mixing between seawater and freshwater, and (3) to determine water-rock reaction processes by using geochemical indicators.
2. Study Area2.1. Hydrogeological Setting
This study examined three cores located in North Jiangsu coastal plain, China, adjacent to the south Yellow Sea (SY1: 33.82°N, 120.43°E; SY2: 33.80°N, 120.32°E; SY3: 34.40°N, 119.81°E) (Figure 1(b)). The North Jiangsu coastal plain has a continental and maritime climate, with a mean annual temperature of 13–16°C. The average annual precipitation is approximately 800–1,200 mm, and nearly 30–60% of the annual precipitation falls from June to September. The annual evaporation is 900–1,050 mm.
Location map of the research area. The A-A′ line indicates the location of the cross-section displayed in Figure 2.
A-A′ cross-section of the study area. SY1 and SY2 cores are displayed in the cross-section.
That transgressions and regressions occurring frequently in the Quaternary period have been confirmed in published studies [25, 26]. The depositional facies of aquitards change from alluvial and lacustrine to marine sediments from west to east on the coastal plain, mainly because of the development of the Huai River and Yellow River, climate change, and transgressions and regressions [27, 28]. The major aquitards are widely dispersed in this study area and generally consist of Quaternary sediments with thickness of 15–75 m in the east and up to 100 m in the west [28]. Four porewater types, namely, salt-fresh water inversion, salt-brackish water inversion, salt-fresh-brackish water superimposition, and salt-fresh-salt-fresh superimposition, were confirmed on the basis of the ions and the TDS analysis of porewater from the southern borehole sediments of this region [29]. Data from the northeastern regions were not readily available.
Quaternary aquifers are observable from the single aquifer to multiaquifer system from the west to the east [30]. Three partially connected confined aquifers exist in Quaternary sediments of this study area (Figure 2). These aquifers are composed of sandy gravel, medium-fine sand, and fine sand and are separated by silt- and clay-dominated aquitards. Aquifer groundwater flow in the confined aquifers is generally from west to east under natural conditions with relatively low velocity and with a hydraulic gradient of approximately 2/10,000 [4]. Geochemical analysis indicates that phreatic aquifer can be divided into two zones: the fresher zone (a thin layer about 2.0–2.5 m thick below the water table) and the saline zone. Phreatic aquifer is trapped in marine-terrigenous deposits that correspond to the Holocene (Q4) formation. Phreatic aquifer in the fresher zone generally has TDS concentrations within 1-2 g/L and tends to be less saline than the water in the underlying saline zone as a result of the diluted influence of atmospheric precipitation. Phreatic aquifer in the saline zone and shallow confined aquifers are saline with TDS mostly over 3 g/L. The TDS concentrations of saline water exhibit an increasing trend from the west to east. The saline aquifer groundwater mainly originated from seawater and then was mixed with the recharge of precipitation and irrigation water [4]. Groundwater in the second confined aquifer exhibits the chemical signature of freshwater and was trapped in the fluvial deposits during the Middle Pleistocene (Q2) [30]. The third confined aquifer groundwater has relatively high TDS (>1 g/L), which could be related to residual seawater [22]. Additionally, leakage recharge from the upper shallow aquifer groundwater may play a key role in contributing to the source of deep confined groundwater during groundwater exploitation.
2.2. Local Stratigraphy
Three cores were drilled on-site between 2013 and 2015, to depths of 250 m, 120 m, and 350 m, respectively. The stratigraphy of the study boreholes was composed of Neogene and Quaternary deposits (Figure 3). Clay, silty clay, silt, and sand were noted in the field drilling logs that described the observed core sediments. The dominant sediments generally consisted of fine-grained deposits derived from an alluvial fan.
Stratigraphic columns of the drilled three boreholes in North Jiangsu coastal plain. The observed deposit distributions across the cores are illustrated.
The upper 2-3 m of the deposit was visibly fractured and weathered, and below that layer, the deposits were grey with no visible fractures. Various lithological facts and the aquifers distribution can be observed from Figures 2 and 3. The approximate chronostratigraphic framework of three cores in the column was provided on the basis of Kang [30] and Zheng [31]. The Early Pleistocene sediments were primarily composed of clay, medium-coarse sand, and silt-fine sand interlayers. The Middle and Late Pleistocene layer mainly consisted of clay with silt-fine sand interlayers. The Holocene layer is dominated by clay and silty clay.
The mineralogy of the clay-size fraction of the aquitards from the middle reaches of Huai River and southern Yellow Sea is composed of 52%–82% illite, 6%–16% chlorite, 5%–14% kaolinite, and 1%–15% smectite [32, 33], with only minor amounts of quartz [34]. Illite (17–76%) and chlorite (6%–17%) are the major components of these core samples; a certain amount of quartz, albite, and calcite is also found in the deposits.
3. Methods
Field work was conducted on-site using a rotary drill. Between January and February 2013, continuous core samples were collected at borehole SY1. To improve the sample density for the reliability of hydrochemical data and obtain other high-salinity porewater samples in deep clay-rich sediments, additional continuous core samples were collected at the SY2 and SY3 cores in December 2014. The distances between borehole SY1 and borehole SY2, borehole SY1 and borehole SY3, and borehole SY2 and borehole SY3 were 10 km, 86 km, and 81 km, respectively.
Forty-four core samples (length: 20 cm; diameter: 8 cm) were collected in SY1 core, 26 samples in SY2, and 15 samples in SY3 core to extract porewater for chemical analysis. The collected samples were packaged in a sealable aluminum barrel and sealed with wax. Immediately after retrieval and before packing, the outer 2-3 cm of the samples was removed and discarded, excluding the contamination of the drilling fluid. The collected core samples were squeezed using high pressure to obtain porewater samples for chemical analyses.
From the boreholes SY2 and SY3, 25 and 11 core samples were collected on-site by 10 mL glass bottles to extract porewater for stable isotopic analyses. Subsequently, the selected samples were sealed using raw adhesive tape to avoid moisture loss. All samples were stored in insulated coolers at approximately 4°C prior to analysis to minimize the growth of microorganisms.
3.1. Porewater Extraction Method
Initially, researchers mostly focused on high-porosity clays, because porewater could be easily extracted from those clays by installing piezometers in shallow aquitards (at depths not greater than 100 m) in situ [35, 36]. Recently, for research into the deep, low-permeability clays, concentrations of porewater were mostly based on the ion concentrations of leach solutions from core deposits. Then, ion concentrations of the leach solutions were converted to in situ porewater concentrations [37, 38]. However, many parameters (e.g., dry bulk density and effective porosity) are required in this method and the parameters retained some errors during the test. In this study, the thicker cohesive soil was focused on. To avoid any transforming error, an ex situ mechanical squeezing method was used in this study to collect porewater [21, 39]. Various studies have described such mechanical squeezing methods in detail [40, 41].
The design of the squeezing device was based on the squeezing method of Patterson et al. [41] and N2 pressurizing equipment was added and applied to a piston to produce pressures of up to 8 MPa (Figure 4). The squeezing device consisted of an Inox chamber, two porous plates, a top cover, a bottom cover, several screws, and a N2 gas holder. To prevent any effect from air exposure and ionic reactions on porewater composition, a stainless steel cylinder contained the sediments. The connections between the top and bottom endplates and the Inox chamber were sealed with O-rings. Böttcher et al. [42] suggested that porewater chemistry could be modified when pressures were as high as 60 MPa due to the collapse of soil structures. Thus, this method is considered suitable for fine-grained and organic-matter-sediments [41].
The photograph of the squeezing device.
The outermost 1-2 cm of the prepared core samples was scraped to discard any materials that could have been altered by exposure to the atmosphere (oxidation). Each sample (approximately 1000 g) was trimmed into separate pieces to fit the sample chamber. About 10–20 min was required until the core samples were inserted into the sample chamber. Some tens of millilitres of air might exist in the space of sample chamber; any air present in the chamber would be expelled during the first half hour of each sample squeezing test. Furthermore, the water collection bottle was assembled to the outlet hole after this phase to eliminate air in the bottle. To avoid overconsolidation or destruction of the clay-pore system, the applied stress was increased to 8 MPa gradually rather than in a single step. First, a small stress of approximately 1 MPa was initially applied to expel most of the air in the sample chamber and to ensure that the sample was bedded. Second, the applied stress was increased by 1 MPa every 4-5 hours during the day. Finally, after about 72–96 h, approximately 8–70 mL of porewater was collected in the clean bottle from the center hole of the bottom endplate. The water samples were weighed and filtered through a 0.45 μm membrane and immediately stored in a refrigerator at 4°C before being analyzed.
Because of the high salinity (>1 g/L) of many porewater samples, vacuum distillation was selected to collect porewater for stable isotope analysis. The vacuum distillation technique is widely used to extract porewater from core samples. If the extraction had not been complete, it might have produced a systematic depletion in heavy isotopes due to the influence of isotope fractionation [43]. The extraction temperature, sample grain size, and the extraction time were the major factors that influenced the porewater extraction [44, 45]. In addition, it is suggested that the depletion tended to increase as the water content of the deposits decreased, especially for clay samples [44, 46]. The fractionation effects can be negligible when the sediment samples have water content proportions that exceed 0.15 [45, 47]. In this study, each of the samples (about 15 g) with high water content (0.21–0.27 by weight) was placed in an airtight glass tube and connected to the modified extraction line. All extraction apparatuses were made of glass. Porewater was extracted through vaporization at approximately 120°C and heated for 6-7 h in vacuum. The extraction time was verified to ensure complete extraction.
3.2. Water Sample Analysis
All extracted porewater samples were filtered through a 0.45 μm membrane. High-salinity porewater was diluted using deionized water before chemical analysis. Samples for cation analysis were acidified with 6N HNO3 to pH < 2. Samples for chemical determinations were analyzed at the School of Environmental Studies of the China University of Geosciences. Major cations and anions were determined using a plasma optical emission spectrum analysis instrument (ICP-OES: ICAP-6300) and an ion chromatograph (ICS-1100), respectively, with an analytical reproducibility lower than 0.01 mg/L. HCO3- was measured using the titration method.
Stable isotope analysis was conducted at the Laboratory of the Geological Survey Institute, China University of Geosciences in Wuhan. Stable isotope values of δ18O and δ2H were measured using a liquid water isotope analyzer (LGR, IWA-45 EP). In the course of the experiment, three standard solutions (see Table 1) will be used to check the stability of the system for stable isotope measurements. Every standard solution was tested six times, and the mean measured values of the middle 4 times were used to establish standard curve. For stable isotope analysis, porewater samples were rerun six times, and then the average value of the middle 4 times will be obtained based on the standard curve. The replicates were normalized to Vienna Standard Mean Ocean Water (V-SMOW) [48] on the basis of (1) and presented in δ‰. The standard deviation of the replicates was less than 0.1‰ for δ18O and 0.5‰ for δ2H. Based on the standard deviations for replicates and the standard solutions, the analytical precision was calculated to be 0.18‰ for δ18O and 0.71‰ for δ2H.(1)δsample‰=RsampleRstandard-1×1000‰,where δsample is the isotope ratio of sample relative to the VSMOW and R is the H2/H1 or O18/O16 atomic ratio.
The laboratory (LGR) standards and their assigned values.
LGR3C
LGR4C
LGR5C
δ2H = −97.3 ± 0.5‰ versus
δ2H = −51.6 ± 0.5‰ versus
δ2H = −9.2 ± 0.5‰ versus
VSMOW
VSMOW
VSMOW
δ18O = −13.39 ± 0.15‰ versus
δ18O = −7.94 ± 0.15‰ versus
δ18O = −2.69 ± 0.15‰ versus
VSMOW
VSMOW
VSMOW
3.3. Seawater Mixing Content
Groundwater in the North Jiangsu coastal plain was often affected by seawater during the transgressions. Determining the proportion of seawater in the groundwater could be used to understand the hydrological geochemical process. The extent of seawater mixing was calculated between two end-member samples by using a mass-balance model [8].
In this study, Cl− is used as the tracer to estimate the seawater mixing proportion. The seawater contribution (fsea) is calculated as follows:(2)fsea=CCl,sample-CCl,fCCl,sea-CCl,f,where CCl,sample is the measured concentration of Cl− in groundwater and CCl,sea and CCl,f are the concentrations of Cl− in the standard seawater and the selected freshwater, respectively.
The seawater contribution fsea is used to calculate the theoretical concentration of component j:(3)Cj,mix=fsea∗Cj,sea+1-fsea∗Cj,f.
The ionic delta (ΔCj) reflects a comparison between the measured concentration and its theoretical concentration (Cj,mix):(4)ΔCj=Cj,sample-Cj,mix.
4. Results4.1. Salinity Distribution of Water Samples
The chemical compositions of the porewater samples are displayed in Figure 5. The TDS of the porewater in the aquitards is lower than current seawater (35.0 g/L), ranging from 0.03 to 26.1 g/L. On the basis of the various tested TDS concentrations, porewater in the aquitards is divided into three water types: fresh (TDS of <1 g/L), brackish (TDS of 1–3 g/L), and saline (TDS of 3–30 g/L) water. Fresh porewater is completely absent from the SY3 profile, but presents in cores SY1 and SY2. Three water-bearing subsystems termed the “shallow, middle, and deep systems” had been defined in this study area [29]. Accordingly, porewater samples were classified into three groups.
Solute concentrations and stable isotopes (δ18O, δ2H) of porewater versus depth.
(1) Group 1 corresponds to porewater samples that were collected from clay-rich sediments in the shallow water system with depths of 2.1–62.6 m (SY1), 10.6–83.8 m (SY2), and 6.4–80.0 m (SY3). These porewater samples were highly saline, and their TDS concentrations ranged from 0.85 g/L to a maximum of 26.1 g/L.
(2) Group 2 corresponds to porewater samples from clay-rich sediments in the middle water system at depths of 62.6–180.0 m (SY1) and 83.8–120.0 m (SY2). Most of these porewater samples indicated the chemical signature of the group with TDS values below 1.0 g/L.
(3) The samples of Group 3 were classified as mixtures of brackish and saline porewater samples taken from deep water system, at depths of 180–250 m (SY1) and 80–253.6 m (SY3). The TDS concentrations were 1.2–7.0 g/L, and some small variations of ion concentrations with depths were evident in each core.
4.2. Ion Trends
Cl−, Na+, Ca2+ Mg2+, SO42-, and Br− concentrations in the porewater chemical compositions along the depth were presented in Figure 5. The dominant anion in most porewater samples was Cl− (up to 98% of the total anion concentration), followed by SO42-, especially in saline and brackish porewater samples. The dominant cation in most porewater samples was Na+ (up to 96% of the total anion concentration), followed by Ca2+ and Mg2+. In Group 1, their concentrations surged at a depth of approximately 20–50 m and then decreased rapidly. The maximum Cl− value reached 455.2 meq/L, which is higher than that of the near-shore seawater values (420–448.7 meq/L). In Group 2, the porewater samples with the lowest salinity had the characteristics of freshwater. Group 3 had brackish and saline porewater characteristics, and the total ion concentrations are higher than Group 2 samples.
Both Cl− and Br− concentrations are valuable indices for investigating the origin and evolution of groundwater [49]. In this study, Group 2 samples had relatively low Br− values, whereas Group 1 had higher values, some (0.83–1.10 meq/L) of which even exceeded the Br− concentrations of standard seawater (0.81 meq/L [50]) (Figure 5). The Group 1 samples had a wide range of Cl/Br molar ratios ranging from 181.6 to 904.8. Except for some samples from the SY1 core, most of the Cl/Br molar ratios were close to or slightly below the near-shore seawater values, as depicted in Figure 6(c). The molar Cl/Br ratios in most Group 3 samples were relatively variable and lower than the standard seawater value and ranged from 308.9 to 726.5.
Evolution of hydrochemical ratios with increasing salinity of porewaters. The near-shore seawater samples were given by Fu [22]. Brine samples were collected from Han et al. [10]. R indicates ratios in rainwater, based on data from Zhang et al. [51]. The seawater dilution line between standard seawater and fresh aquifer groundwater is illustrated in red line, and the numbers along with the mixed line stand for the percentage (%) of seawater mixing. The seawater evaporation trajectory (red dashed line) was given by Fontes and Matray [52]; SW, G, H, E, and S represent points of standard seawater, precipitation of gypsum, halite, epsomite, and sylvite, respectively.
As illustrated in Figure 6(a), Group 1 and Group 3 saline porewaters were characterized by low Na/Cl ratios compared with Group 2 samples. For brackish porewater samples in Group 3, Na/Cl ratios were higher than normal seawater. By contrast, Group 2 samples presented a wide range and exceeded the seawater and rainwater value. The Ca/Cl ratios in most Group 1 and 3 porewater samples were higher than seawater. For Group 2 samples, Ca/Cl ratios were below that of the rainwater (Figure 6(b)).
4.3. Stable Isotope Compositions (δ18O and δ2H)
The δ18O and δ2H values in Group 1 varied between −8.5‰ and −2.1‰, and δ2H values were between −54.6‰ and −17.3‰. The observed stable isotope values generally decreased with depth (Figure 5). Group 2 samples had narrow ranges, between −8.5‰ and −7.6‰ for δ18O and between −54.6‰ and −53.2‰ for δ2H. The δ18O and δ2H values in Group 3 samples ranged from −7.1‰ to −9.4‰ and from −56.8‰ to −66.9‰, respectively.
A local meteoric water line (LMWL, δ2H = 5.8δ18O − 6.79, n=44, and R2 = 0.87, Figure 7) was estimated on the basis of the average monthly rainfall monitored from January to September between 1986 and 1990 at Yantai station, which is the nearest coastal station with similar climatic conditions [53]. The stable isotope compositions of aquifer groundwater and surface water were taken from Fu [22] and Ha et al. [28] (Figure 7) to identify the sources of porewaters. In the plot of δ2H against δ18O, three groups of points, surface water, and groundwater in the aquifers were close to the LMWL. The δ18O and δ2H values in the Group 3 samples and the deep confined groundwater of the Huai River basin were characterized by more depleted δ-values compared with Group 1 and Group 2 porewater.
Plots of δ2H versus δ18O for porewater, aquifer groundwater, and river water samples. Samples with (α) represent the water samples in the study area which were collected from Fu [22] based on aqueous analyses. Samples with (β) represent aquifer groundwater samples in the Huai River basin that were obtained from Ha et al. [28]. The LMWL (δ2H = 5.8δ18O-6.79) was plotted based on the δ18O and δ2H values of the average monthly rainfall at the nearest coastal station. The mixing line between the standard seawater and fresh aquifer groundwater (δ2H: −60.1‰; δ18O: −8.4‰) is depicted in red.
5. Discussion5.1. Origin of Porewater
Group 1 samples had δ18O and δ2H values between −8.5‰ and −2.1‰ and between −54.6‰ and −17.3‰, which were more enriched than those of the groundwater in the first confined aquifer beneath the aquitard, as depicted in Figure 7. The δ18O and δ2H values of Group 1 samples were mostly scattered around the LMWL, suggesting that atmospheric precipitation was their predominant origin. The maximal δ18O and δ2H values in Group 1 were −2.1‰ and −18.0‰ consistent with the values of near-shore seawater. Based on the mass-balance mixing model (see (2)), the seawater mixing proportions were estimated to be between 2% and 82% in Group 1. In addition, the plots of δ18O and δ2H fell along the mixing line between standard seawater and fresh aquifer groundwater. This indicates that the Group 1 porewater originated chiefly from seawater and was subsequently diluted by the infiltration of meteoric water during a regression period. The major ion concentrations and δ18O and δ2H values tended to decrease with depth between 20 m and 80 m, proving that the sooner the saline pore water develops, the stronger the dilution effect is [54]. Furthermore, no obvious discontinuous vertical distribution occurred in the stable isotopic profiles. Presumably, the Group 1 samples were derived from seawater formed in the Holocene period around 10 ka BP–1,128 AD, because seawater from the Holocene period had generally retreated from the study area by 1,128 AD (Figure 1(b)) [55].
By contrast, Group 2 samples with depleted stable isotope values deviated from the mixing trend line, indicating no evidence for the involvement of seawater. On the other hand, these samples had Br− concentrations below detection, further corroborating freshwater recharge for Group 2. Their δ18O values ranging from −8.5‰ to −7.6‰ were close to the LMWL, suggesting Group 2 porewater was mainly derived from atmospheric water. Some overlaps in Group 2 porewater and the second confined groundwater samples were observed (Figure 7), indicating that the vertical diffusive mixing between the porewater and aquifer groundwater possibly existed. Based on an analysis of carbon isotopic data (14C), groundwater ages in the second and third aquifers turned out to be around 30 ka BP [22, 30]. Thus, recharge under cold climate conditions in the Late Pleistocene (e.g., during the glacial period) is supported by the depleted and identical stable isotopic signatures.
The δ18O and δ2H values of Group 3 samples and deep aquifer groundwater of the Huai River basin were similar; however, the TDS concentrations for most of the deep aquifer groundwater in the Huai River basin were less than 1.0 g/L [28]. Two transgressions (MIS 5, about 110–70 ka BP; MIS 1, around 10 ka BP) had been reported and accepted as having occurred in this study area [56, 57]. Given that the estimated age of groundwater in the overlying third confined aquifer was around 30 ka BP [22, 30], Group 3 samples with high salinity (TDS: 1.2–7.0 g/L) were mainly attributed to seawater origin in transgression with an age of 70–110 ka BP during the Late Pleistocene. As shown in Figure 6(d), the proportion of seawater mixing in porewater was less than 20%, and most of the samples contained less than 5%. The fact that these samples were diluted by palaeofreshwater with depleted stable isotopes similar to those of the deep aquifer groundwater in the Huai River basin is supported by the fact that a small contribution of seawater has a great and measureable influence on water chemistry but a limited effect on δ18O and δ2H values [58].
5.2. Salinization of Porewater
Old brine has been identified in Laizhou Bay, approximately 200 km north of the study area [8, 59] (Figure 1(a)). Old brine was confirmed to originate from ancient seawater, and evapotranspiration was the dominant mechanism of brine formation [10]. Based on the 14C data, the old brine was deduced to form in the Holocene over a period of approximately 2.3–7.0 ka BP. Around the same time, some palaeolagoons formed in the Jiangsu coastal plain [60, 61]. After the sea level dropped, the high lagoon surface and abundant surface water (e.g., many branches of the Huai River and palaeo-yellow river) resulted in the brine being flushed out [60]. In addition, based on the archaeological evidence, a large number of marine salterns were present along the paleocoastline between 1,300 and 1,900 AD [62]. Old brine likely existed in the study area and later was eventually displaced, just like Qingdao (another area in the coastal plain) [63].
A well-correlated relationship between Cl− and Br− was observed, and most Cl/Br molar ratios in Group 1 samples were scattered around the dilution line of standard seawater (Figures 6(c) and 6(d)). It is deemed that the salinity and Br− of Group 1 samples were mainly derived from seawater. Some samples had Cl/Br molar ratios less than the normal seawater. The SY1 core samples were likely influenced by the old brine, leading to the low Cl/Br ratios with respect to normal seawater. Considering that Group 1 samples were mainly collected from marine sediments, the release of Br− from diagenesis of marine organic material, the adsorption of Br− on clay minerals, and selective concentration of Br− through the clay membrane effects may also be responsible for the observed low Cl/Br ratios [50, 64]. Furthermore, two samples at depths of 15.0 and 17.1 m had notably rich Br− concentrations, consistent with the elevated NO3- concentration (up to 0.99 meq/L), which seems to probably be also attributed to artificial factor (such as pesticides) [44].
The Cl/Br ratios in most of Group 3 samples were low with respect to normal seawater, as illustrated in Figure 6(c). In the diagram of Br− and Cl−, most of the samples tended to scatter around the mixing lines. A linear relationship between Br− and Cl− was observed, with some Br− enriched relative to the expected path for standard seawater and freshwater. Groundwater in the adjacent aquifers and deep aquifer groundwater in the Huai River basin were proved to mainly contain brackish water or freshwater from the Late Pleistocene [22, 28, 30]. In addition, the upper thick low-permeability clay layers may have the ability to impede the leakage of brine formed in the Holocene. For the porewater in the deep aquitards, a contribution of brine to porewater composition seems less likely in comparison with the groundwater in the adjacent aquifers. The salinity possibly was derived from the relict palaeoseawater. The low Cl/Br ratios relative to seawater might be related to diagenesis of marine organic matter with Br− release, the adsorption of Br− on clay minerals, and selective concentration of Br− through the clay membrane effects [50, 64].
5.3. Water-Rock Interaction
Cl− was deemed to behave relatively conservatively and to be a useful tracer to distinguish the origins of salinity [65]. To highlight the different diagenetic reactions contributing to porewater compositions, some plots of relations among major elements and ion/Cl molar ratios are displayed in Figure 6.
As indicated in Figure 6, Na/Cl ratios in the Group 1 and Group 3 porewater samples ranging wildly from 0.31 to 2.90 relative to seawater (0.86) indicated the occurrence of both Na+ enrichment and depletion processes. The Ca/Cl ratios were mostly higher than seawater. To identify the processes that modified the original contents of the Group 1 and Group 3 samples, the ionic deltas of Na+ and Ca2+ were calculated based on (2), (3), and (4).
As illustrated in Figure 8(a), the values of ΔCa had inverse linear correlations with ΔNa in the Group 1 porewater, with slopes of −0.41 (R2 = 0.60) and −0.66 (R2 = 0.97) for the Group 3 porewater. These results indicate that cation exchange plays a crucial role in Na+ depletion. Na+ adsorbed on the clay surface and albitization may also be responsible for the high Na+ deficiency [66].
(a) ΔCa versus ΔNa; (b) (Ca2++Mg2+-HCO3--SO42-) versus (Na++K+-Cl−) of aquitard porewater. Symbols are the same as Figure 6.
In contrast to the Group 1 and Group 3 samples, the Group 2 samples had high Na/Cl molar ratios ranging from 1.00 to 6.93, whereas lower Ca/Cl ratios relative to rainfall were observed (Figures 6(a) and 6(b)). To examine the cation exchange process of porewater, the plot of (Ca2+ + Mg2+-HCO3--SO42-) versus (Na+ + K+-Cl−) was used [67]. As illustrated in Figure 8(b), most of the samples were plotted around the −1 : 1 line, implying that cation exchange was mainly responsible for the elevated Na+. The portion of the samples plot below the −1 : 1 line implies that the minor excess (Na+ + K+) was contributed by hydrolysis of silicate minerals (e.g., 2NaAlSi3O8(Albite) + 2CO2 + 11H2O → 2Na+ + Al2Si2O5OH4(Kaolinite) + 4H4SiO4 + 2HCO3-) because the aquitards of the study area contain abundant silicate weathering products such as kaolinite, illite montmorillonite, and mica minerals [8, 34] and the release of Na+ from exchange sites of clay minerals [10]. Some other points are above the −1 : 1 line, indicating that the additional (Ca2+ + Mg2+-HCO3--SO42-) inputs may have been derived from the weathering of Ca-bearing plagioclase and calcite, while yet again cation exchange was most noticeable for Na+ enrichment and Ca2+ depletion in the Group 2 porewater.
6. Conclusions
This study provides insights into the origin and evolution of the porewater in aquitards using hydrochemical and stable isotopic evidence based on samples from three deep boreholes drilled into the North Jiangsu coastal plain, China. Three groups were classified by the water-bearing subsystems of “shallow, middle, and deep system” in this study area.
Regarding porewater samples, the δ2H and δ18O values were plotted around the LMWL, indicative of signatures of atmospheric precipitation. Chemical data indicated significant variations in salinity and hydrochemical constituents of porewaters, which were interpreted in terms of different origins. The trapped seawater during Holocene transgressions was the major source producing salinity of the Group 1 samples. The Group 2 samples had freshwater characteristics and resulted from atmospheric precipitation in a cold climate. According to the analysis of the plots of Br− versus Cl− and δ2H versus δ18O, it is deduced that the Group 3 samples evolved from the remaining palaeoseawater from the Late Pleistocene and were subsequently diluted by palaeofreshwater.
Water-rock interaction also accounted for the differences between the initial and current porewater solutions as a result of the long residence times of the water samples. Na+ adsorbed on clay surfaces and albitization were likely associated with the Group 1 and Group 3 constituents, as well as cation exchange between Ca2+ and Na+. By contrast, cation exchange was more noticeable for the Group 2 samples with low salinity. The dissolved Ca2+ is exchanged mainly for Na+, and weathering of silicate minerals and Ca-bearing plagioclase may provide only a minor contribution.
The present study attempted to ascertain the distributions of porewaters, the origins of salt, and the hydrochemical evolution in porewaters with the long residence times. The results have crucial implications for assessing the salinization of porewater in coastal aquitards. The observations of this study could serve as a reference for porewater quality in the coastal aquitards and verify that thick aquitards have a function for preserving ancient water.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
Financial support was provided by the National Natural Science Foundation of China (Grant nos. 41502231 and 41272258). The authors thank Yalei Liu and Yunquan Wang for their useful comments and advice, as well as all members of the Geological Survey of Jiangsu Province, China, for their assistance and support in the process of sampling and data collection.
PostV. E. A.Fresh and saline groundwater interaction in coastal aquifers: is our technology ready for the problems ahead?200513112012310.1007/s10040-004-0417-22-s2.0-17044410545ShiM. Q.2012Master thesis of Shanghai Normal UniversityShiL.JiaoJ. J.Seawater intrusion and coastal aquifer management in China: a review2014728281128192-s2.0-8490770424210.1007/s12665-014-3186-9ShiX.JiangF.FengZ.YaoB.XuH.WuJ.Characterization of the regional groundwater quality evolution in the North Plain of Jiangsu Province, China2015747558756042-s2.0-8494118384510.1007/s12665-015-4575-4Pulido-LeboeufP.Seawater intrusion and associated processes in a small coastal complex aquifer (Castell de Ferro, Spain)20041910151715272-s2.0-294271542310.1016/j.apgeochem.2004.02.004SchiavoM. A.HauserS.PovinecP. P.Stable isotopes of water as a tool to study groundwater-seawater interactions in coastal south-eastern Sicily20093641-240492-s2.0-5764917165110.1016/j.jhydrol.2008.10.005LarsenF.TranL. V.Van HoangH.TranL. T.ChristiansenA. V.PhamN. Q.Groundwater salinity influenced by Holocene seawater trapped in incised valleys in the Red River delta plain201710537638110.1038/ngeo2938HanD. M.SongX. F.CurrellM. J.YangJ. L.XiaoG. Q.Chemical and isotopic constraints on evolution of groundwater salinization in the coastal plain aquifer of Laizhou Bay, China201450812272-s2.0-8488764262510.1016/j.jhydrol.2013.10.040AquilinaL.LadoucheB.DoerfligerN.SeidelJ. L.BakalowiczM.DupuyC.Le StratP.Origin, evolution and residence time of saline thermal fluids (Balaruc springs, southern France): Implications for fluid transfer across the continental shelf20021921-21212-s2.0-003712166510.1016/S0009-2541(02)00160-2HanD.PostV. E. A.SongX.Groundwater salinization processes and reversibility of seawater intrusion in coastal carbonate aquifers2015531106710802-s2.0-8494810644810.1016/j.jhydrol.2015.11.013MongelliG.MonniS.OggianoG.PaternosterM.SinisiR.Tracing groundwater salinization processes in coastal aquifers: a hydrogeochemical and isotopic approach in the Na-Cl brackish waters of northwestern Sardinia, Italy2013177291729282-s2.0-8488090391710.5194/hess-17-2917-2013CaoG.ZhengC.ScanlonB. R.LiuJ.LiW.Use of flow modeling to assess sustainability of groundwater resources in the North China Plain20134911591752-s2.0-8487653744710.1029/2012WR011899HendryM. J.WassenaarL. I.Controls on the distribution of major ions in pore waters of a thick surficial aquitard200036250351410.1029/1999WR9003102-s2.0-0033959666ThorstensonD. C.FisherD. W.CroftM. G.The geochemistry of the fox hills‐basal hell creek aquifer in southwestern North Dakota and northwestern South Dakota19791561479149810.1029/WR015i006p014792-s2.0-0018704835KhaskaM.Le Gal La SalleC.LancelotJ.teamA.MohamadA.VerdouxP.NoretA.SimlerR.Origin of groundwater salinity (current seawater vs. saline deep water) in a coastal karst aquifer based on Sr and Cl isotopes. Case study of the La Clape massif (southern France)2013371021222710.1016/j.apgeochem.2013.07.0062-s2.0-84883176362NajibS.FadiliA.MehdiK.RissJ.MakanA.GuessirH.Salinization process and coastal groundwater quality in Chaouia, Morocco201611517312-s2.0-8495026833110.1016/j.jafrearsci.2015.12.010CyrM.DaidiéA.Optimization of a high-pressure pore water extraction device2007782, article 02390610.1063/1.2428273WassenaarL. I.LeonardI.HendryM. J.Improved piezometer construction and sampling techniques to determine pore water chemistry in aquitards20103745645712-s2.0-0032780880HendryM. J.WassenaarL. I.KotzerT.Chloride and chlorine isotopes (36Cl and δ37Cl) as tracers of solute migration in a thick, clay-rich aquitard system20003612852962-s2.0-003398813410.1029/1999WR900278TimmsW. A.HendryM. J.Quantifying the impact of cation exchange on long-term solute transport in a clay-rich aquitard20073321-21101222-s2.0-3384557147710.1016/j.jhydrol.2006.06.025LiJ.LiangX.JinM.XiaoG.HeJ.PeiY.Geochemistry of clayey aquitard pore water as archive of paleo-environment, western Bohai Bay20152634454522-s2.0-8493081795510.1007/s12583-014-0491-xFuC. C.2015Jilin UniversityFuY. L.Evaluation of groundwater resource forecast in Yancheng2005265659YiL.Reconciling late Quaternary transgressions in the Bohai Sea, China to the global sea level changes, and new linkage of sedimentary records to three astronomical rhythms2013YangH. R.ChenX. Q.Quaternary transgressions, eustatic changes and shifting of shoreline in east China198555980WangY.ZhangZ. K.ZhuD. K.YangJ. H.MaoJ. L.LiS. H.Riversea interaction and the North Jiangsu Plain formation200626301320301320XueC. T.ZhouY. Q.WangG. L.Reviews of the Yellow river delta super lobes since 700 BC2003232330HaC. Y.ZhuJ. Q.YeN. J.HuangJ. J.GongJ. S.The forgotten delta the formation and evolution of Huaihe Delta River20052410941106ZhouH. C.Origin, distribution, and evolution of salt-fresh groundwater in North Jiangsu Basin, China1988121152210.1007/BF025748222-s2.0-0024055753KangB.2017Jilin UniversityZhengG. Y.Quaternary stratigraphic division of the hole QC2 in the southern Huanghai Sea19884513ZhangJ. Q.LiuJ.KongX. H.XueC. T.LiuX. B.SunL. S.YueB. J.WenC.Clay mineral assemblange of sediments from middle reach of Huaihe River2011312129ChenZ. X.Clay minerals from drill hole QC1 in South Yellow Sea19903ShuQ.ChenY.ZhangM. H.ZhaoZ. J.Environmental significance and clay mineral characteristics of YZ1 core sediments in Zhouzhuang, Northern Jiangsu Plain20122238243DesaulniersD. E.DonaldE.CherryJ. A.Origin and movement of groundwater and major ions in a thick deposit of Champlain Sea clay near Montreal19892618089HendryM. J.WassenaarL. I.Implications of the distribution of σD in pore waters for groundwater flow and the timing of geologic events in a thick aquitard system19993561751176010.1029/1999WR9000462-s2.0-0033045024HendryM. J.HarringtonG. A.Comparing vertical profiles of natural tracers in the Williston Basin to estimate the onset of deep aquifer activation2015508649665062-s2.0-84905636339WaberH. N.SmellieJ. A. T.Characterisation of pore water in crystalline rocks2008237183418612-s2.0-4504908423910.1016/j.apgeochem.2008.02.007MazurekM.Alt-EppingP.BathA.GimmiT.Niklaus WaberH.BuschaertS.CannièreP. D.CraenM. D.GautschiA.SavoyeS.VinsotA.WemaereI.WoutersL.Natural tracer profiles across argillaceous formations2011267103510642-s2.0-7995627458910.1016/j.apgeochem.2011.03.124RidoutP. S.PagettR. M.A large volume sediment squeezing system for the analysis of rare earth elements in deep ocean pore waters19841531932012-s2.0-4854910835510.1016/0304-4203(84)90017-3PattersonR. J.FrapeS. K.DykesL. S.McLeodR. A.A coring and squeezing technique for the detailed study of subsurface water chemistry197815116216910.1139/e78-0142-s2.0-0017909348BöttcherG.BrumsackH.-J.HeinrichsH.PohlmannM.A new high-pressure squeezing technique for pore fluid extraction from terrestrial soils1997943-42892962-s2.0-0030613913SacchiE.MichelotJ.-L.PitschH.LalieuxP.AranyossyJ.-F.Extraction of water and solutes from argillaceous rocks for geochemical characterisation: methods, processes, and current understanding200191173310.1007/s1004000001132-s2.0-0345776846Araguás-AraguásL.RozanskiK.GonfiantiniR.LouvatD.Isotope effects accompanying vacuum extraction of soil water for stable isotope analyses19951681-41591712-s2.0-002941306610.1016/0022-1694(94)02636-PAltinierM. V.SavoyeS.MichelotJ.-L.BeaucaireC.MassaultM.TessierD.WaberH. N.The isotopic composition of pore-water from tournemire argillite (France): An inter-comparison study2007321-720921810.1016/j.pce.2006.02.0472-s2.0-33847743917IngrahamN. L.ShadelC.A comparison of the toluene distillation and vacuum/heat methods for extracting soil water for stable isotopic analysis19921401-43713872-s2.0-002699770810.1016/0022-1694(92)90249-UMoreau-Le GolvanY.MichelotJ.-L.BoissonJ.-Y.Stable isotope contents of porewater in a claystone formation (tournemire, France): assessment of the extraction technique and preliminary results199712673974510.1016/S0883-2927(97)00044-92-s2.0-0031435347GonfiantiniR.Standards for stable isotope measurements in natural compounds197827156455345362-s2.0-001787383110.1038/271534a0EdmundsW. M.MaJ.Aeschbach-HertigW.KipferR.DarbyshireD. P. F.Groundwater recharge history and hydrogeochemical evolution in the Minqin Basin, North West China20062112214821702-s2.0-3375094315910.1016/j.apgeochem.2006.07.016DavisS. N.WhittemoreD. O.JuneF.Uses of chloride/bromide ratios in studies of potable water199836233835010.1111/j.1745-6584.1998.tb01099.x2-s2.0-0032033019ZhangZ. F.ShenL. J.ZhuB.YeJ.WangH. L.ZhangQ. C.Chemical chanracteristics and potential sources of precipitation in Nanjing2015384473482FontesJ. C.MatrayJ. M.Geochemistry and origin of formation brines from the Paris Basin, France: 1. Brines associated with Triassic salts19931091-414917510.1016/0009-2541(93)90068-T2-s2.0-0027837922IAEA/WMOGlobal network of isotopes in precipitation the gnip database viennahttp://www-naweb.iaea.org/napc/ih/IHS_resources_gnip.html, 2006ZhaoJ. C.LiangJ.CaiH. S.Relation between formation of saltfresh water and waterbearing medium in the plain of North Jiangsu1993202527XueC. T.LiuJ.KongX. H.Preliminary study of holocene Huai river delta on west coastal plain of yellow sea2010305892901YiL.YuH.-J.OrtizJ. D.XuX.-Y.ChenS.-L.GeJ.-Y.HaoQ.-Z.YaoJ.ShiX.-F.PengS.-Z.Late Quaternary linkage of sedimentary records to three astronomical rhythms and the Asian monsoon, inferred from a coastal borehole in the south Bohai Sea, China2012329-33010111710.1016/j.palaeo.2012.02.0202-s2.0-84862820260LingS.Changes of the palaeogeographic environment in north Jiangsu Plains since the Holocene1991842028WarnerN. R.JacksonR. B.DarrahT. H.OsbornS. G.DownA.ZhaoK.WhiteA.VengoshA.Geochemical evidence for possible natural migration of Marcellus formation brine to shallow aquifers in Pennsylvania201210930119611196610.1073/pnas.11211811092-s2.0-84864370882XueY.WuJ.YeS.ZhangY.Hydrogeological and hydrogeochemical studies for salt water intrusion on the South Coast of Laizhou Bay, China2000381384510.1111/j.1745-6584.2000.tb00200.x2-s2.0-0033987983LingS.Study on models of the changes about Sheyang lake in historical period20053207379ZhangZ. K.XieL.ZhangY. F.XuJ.LiS. H.WangY.Sedimentary records of the MIS 3 transgression event in the north Jiangsu plain, China2010305883891BaoJ. L.2014Shanghai, ChinaFudan UniversityYou-SongH.Shao-QingW.Guang-IanM.The forming-storing environment and abundance conditions of quaternary underground brine in the embayed coast of Qingdao199715433234110.1007/BF028505672-s2.0-77951500750BensenouciF.MichelotJ. L.MatrayJ. M.SavoyeS.MassaultM.VinsotA.Coupled study of water-stable isotopes and anions in porewater for characterizing aqueous transport through the Mesozoic sedimentary series in the eastern Paris Basin201453881012-s2.0-8489737961710.1016/j.marpetgeo.2013.12.012LiJ.LiangX.ZhangY.LiuY.ChenN.AbubakariA.JinM.Salinization of porewater in a multiple aquitard-aquifer system in Jiangsu coastal plain, China20172582377239010.1007/s10040-017-1622-0KarroE.MarandiA.VaikmäeR.The origin of increased salinity in the Cambrian-Vendian aquifer system on the Kopli Peninsula, northern Estonia20041244244352-s2.0-464426608310.1007/s10040-004-0339-zFisherR. S.MullicanW. F.IIIHydrochemical evolution of sodium-sulfate and sodium-chloride groundwater beneath the Northern Chihuahuan Desert, Trans-Pecos, Texas, USA1997524162-s2.0-000177737410.1007/s100400050102