Abstract

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.

1. 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) [79], 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 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, to ascertain the primary salinity sources of porewater and the mixing between seawater and freshwater, and to determine water-rock reaction processes by using geochemical indicators.

2. Study Area

2.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.

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 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 [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.

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 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. 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.where is the isotope ratio of sample relative to the VSMOW and is the or atomic ratio.

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 () is calculated as follows:where is the measured concentration of Cl in groundwater and and are the concentrations of Cl in the standard seawater and the selected freshwater, respectively.

The seawater contribution is used to calculate the theoretical concentration of component :

The ionic delta (Δ) reflects a comparison between the measured concentration and its theoretical concentration ():

4. Results

4.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.

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.

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.

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+, , 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 , 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.

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, , and = 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.

5. Discussion

5.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 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 ( = 0.60) and −0.66 ( = 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].

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+--) 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., + 2CO2 + 11H2O → 2Na+ + Al2Si2O5 + 4H4SiO4 + 2) 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+--) 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.