Open Access
How to translate text using browser tools
1 March 2016 Seawater-groundwater Exchange in a Silty Tidal Flat in the South Coast of Laizhou Bay, China
Lijia Hou, Hailong Li, Chunmiao Zheng, Qian Ma, Chaoyue Wang, Xuejing Wang, Wenjing Qu
Author Affiliations +

Hou, L. J.; Li, H. L.; Zheng, C. M.; Ma, Q.; Wang, C. Y.; Wang, X. J., and Qu, W. J., 2016. Seawater-groundwater exchange in a silty tidal flat in the south coast of Laizhou Bay, China.

There were few studies about seawater-groundwater exchange in silty, low-permeability tidal flats with very gentle slopes. This paper reports the monitoring data and preliminary analytical results on a typical transect in a silty tidal flat with large-scale seepage faces at the south coast of Laizhou Bay, China. The “pair-wells method”, which was an improvement of the single-well method used by Ma et al. (2015) was used to estimate the seawater-groundwater exchange rate. We selected 14 locations along a typical transect in the intertidal zone, with a slope of 0.4‰ and cross-shore length of 3514 m, to install pair-wells for monitoring the groundwater head, salinity and temperature at two different depths simultaneously once per hour from 20:00 August 10th to 23:00 September 12th, 2014. The vertical hydraulic conductivity measured in situ ranges from 5.4 ×10−7 to 1.1×10−5 m s−1. The salinity of the underground brine ranges from 25 to 56 g l−1. Based on the observed data and the generalized Darcy's law, the submarine groundwater discharge (SGD) and inflow along the entire transect were estimated to be 164.0 and 6.5 m2 d−1, respectively. It is found that 80.4% of the total SGD occurred between W7 and W10, where the hydraulic conductivity is one or two orders of magnitude greater than that at other wells; and 64.6% of the total inflow occurred between W2 and W4. The single-well method, on the other hand, yields a SGD value of 223.0 m2 d−1 and inflow value of 13.3 m2 d−1. Neglecting the density effect may lead to an overestimation of the SGD by 25.6% and underestimation of the inflow by 27.4%.


With the new technologies and methods developed, the study of seawater-groundwater exchange has been a hot spot of hydrogeology and oceanography. Submarine groundwater discharge (SGD) has significant impacts and implications on ecological environment due to its important role in transport of nutrients, contaminations and other chemicals (Li and Jiao, 2013; Li et al., 1999; Slomp and Cappellen, 2004). An increasing number of scholars have paid attention to SGD estimation.

Many different methods have been developed to estimate SGD. Since 1996, natural radioactive isotopes Ra and Rn were widely used as tracers to estimate SGD in large scale areas (Gu et al., 2012; Lee et al., 2012; Moore, 1996, 2007; Moore and Oliveira, 2008; Wang et al., 2015; Xu et al., 2014). Seepage meters were also used to estimated SGD at local area, such as manual seepage instrument developed by Cable et al. (1997), heat pulse seepage meter developed by Taniguchi and Fukuo (1993) and ultrasonic flow meter developed by Paulsen et al. (2001). Numerical simulations usually focused on limited coastal transect, such as Guo et al. (2010), Robinson et al. (2007), Xia and Li (2012), Xia et al. (2010).

In comparison with many existing studies in gravel or sand beaches, there were few studies of seawater-groundwater exchange in silty beaches. A series of experiments conducted in gravel beaches of Prince William Sound, Alaska, indicated that the two-layered structure of a high-permeability surface layer underlain by a low-permeability layer had significant impacts on the persistence of spilled oil (Guo et al., 2010; Li and Boufadel, 2010; Xia et al., 2010). Xia and Li (2012) made a comparative study between a bald beach and a mangrove beach with similar topography and tidal actions in the intertidal zones of Dongzhaigang in Haikou City of Hainan Province, indicating that the inland freshwater recharge and the mud-sand two-layered structure played key roles in the development of mangrove and their spatial distribution.

Ma et al. (2015) developed a simple method (“single-well method”) for estimating seawater-groundwater exchange rate in a silty beach, taking into account the effects of seepage face and density and based on field measurements of groundwater hydraulic head, temperature and salinity. The exchange rate can be calculated by the generalized Darcy's law, the observed tidal level and the hydraulic head in the observation well, Their case study in a silty beach at Haimiao of Laizhou Bay indicated that the SGD and inflow were 8.8 and 15.3 m3 d−1 m−1, respectively. A freshwater discharge tube was also identified near the low-tide line. Riedel et al. (2010) quantified SGD in a tidal sand flat margin with a seepage face in the German Wadden Sea using the numerical code SUTRA, and found that the magnitude of SGD was controlled by the maximum hydraulic gradient during low tide.

In this paper, “pair-wells method” is used to improve the single-well method used by Ma et al. (2015) to investigate the seawater-groundwater exchange rate in a typical cross-shore silty-beach transect in the intertidal zone of Qingxiang Town (Liutuan Town), Laizhou Bay. Fourteen pair-wells were used to monitor the pore water pressure, salinity and temperature at two different depths simultaneously along the entire transect. Based on these observed data, we quantitatively estimated the seawater-groundwater exchange rate. The impacts of density and groundwater pumping for salt production on SGD and inflow were discussed quantitatively. The difference between the results of this study and those of Ma et al. (2015) were analysed. Comparisons were made between the single-well and pair-wells methods. The advantages of the pair-wells method over the single-well one are discussed.


Field Description

Laizhou Bay, an important fishing and salt production base of China (Figure 1), is located in the south of Bohai Sea, north of Shandong Peninsula. The beach in the southern coast of Laizhou Bay (from Hutouya to Yangjiaogoukou) is very flat, and the width of the intertidal zone in the cross-shore direction is about 4000 m. The tidal flat in the study area consists of silty-muddy sediments, and the aquifers are mainly composed of late Pleistocene to Holocene sediments. Brine in Quaternary aquifers is mainly distributed within 10 km of the shoreline (Han et al., 2014; Zhang and Peng, 1998). There are three layers containing brine in the sediments (Xue et al., 2000), which has resulted from a combination of past and recent processes, including a number of phases of regression, transgression, brine formation and mixing (Han et al., 2014; Li et al., 2000; Xue et al., 2000).

Figure 1.

Location of the study area: (a) Bohai Sea; (b) Laizhou Bay; (c) 14 pair-wells (numbered W1, …,W13, and W14 sequentially seaward) were set along the cross-shore transect. The white patch areas in the lower part of (c) are saltpans


The weather in Laizhou Bay belongs to sub-humid monsoon climate, with an annual average temperature of 11.9°C (Han et al., 2014). The annual mean precipitation is 660 mm, and 60–70% of the precipitation generally concentrates in June–August. The annual average potential evapotranspiration, approximately 1774 mm, is much greater than the precipitation (Bi et al., 2012; Han et al., 2014). Under natural conditions, the unconfined aquifer is recharged by precipitation and groundwater from south (Han, Wu, and Zhang, 2011; Han et al., 2014). The area is exposed to semi-diurnal tides, and the average durations for both rising and falling tides are about 6 hours.

The typical transect that we selected in the intertidal zone is located in Qingxiang Town (Figure 1). Along the entire transect, we dug 14 pits to install pair-wells, which are numbered as W1, W2, …, and W14 sequentially seaward (Figure 2). The beach surface elevations of the pair-wells were surveyed by an electronic total station (TS–800, Johanna). These data and the elevations of each pressure transducer installed at the 14 observation pair-wells were given in Table 1. There were several times of raining during the monitoring period (from 20:00 PM August 11th to 23:00 PM September 12th, 2014, almost 33 days), with the precipitation ranging from 28.3 mm to 51.6 mm.

Figure 2.

Vertical cross-section of the studied transect and locations of transducers. A pair-wells, with two LTC-Divers, was installed at every pit of W1, W2, …, and W14 for monitoring pore water pressure, electrical conductivity and temperature


Table 1.

Elevations of beach surface and each LTC-Diver at the 14 observation wells


Field Methodology

Each of pair-wells consists of two LTC-Divers (Level, Temperature and Conductivity Diver; Levelogger Junior, Solinst, Canada), inner and outer PVC pipes, and a plexiglass rod linking the two LTC-Divers together (Figure 3a). The pipes were slotted uniformly with holes about 5.0 mm in diameter to keep a good hydraulic connection around the pair-wells. Furthermore, to prevent holes from being blocked by silty-muddy sediments generally occurring in the silty beach, fine stainless-steel meshes were screened outside the two pipes, and the room between the two pipes was filled with coarse sands (diameter of 0.5 – 2.0 mm). Then, a LTC-Diver was placed in each inner pipe (Figure 3b). In each pit, a pair of wells were installed to record the pore water pressure, temperature and conductivity hourly at two different depths (Table 1). At the same time, a LTC-Diver was installed in the sea to monitor the tidal fluctuation; and a Baro-Diver (Schlumberger) was placed above the ground surface to obtain the air pressure and temperature once an hour.

Figure 3.

(a) Schema of pair-wells. Not to scale. (b) A LTC-Diver was placed in each inner pipe, and the distance between the two LTC-Divers is a constant of 52.8 cm


Based on the definition of salinity given by UNESCO (1966), the electrical conductivity measured by a LTC-Diver can be transformed into salinity. More details can be seen in Cox, Culkin, and Riley (1967); Ma et al. (2015); and, Wooster, Lee, and Dietrich (1969). Then, the observed pressures of air and ground-water, and elevations of each LTC-Diver, were transformed into freshwater-equivalent hydraulic head (Figure 4).

Figure 4.

The tidal head (50%gray line), observed hydraulic heads at the up LTC-Diver (black dots) and low LTC–Diver (40%gray dots) are shown in the primary axis, salinity (g l−1) at the up LTC-Diver (black line) and low LTC-Diver (40%gray line) are shown in the secondary axis, of W1, …, W13, W14 in the intertidal zone. The elevations of each LTC-Diver and beach surface are shown as well


Figure 4.



Figure 4.



Figure 4.



The in-situ multiple diameter falling-head method (Wang et al., 2014) was used to measure the vertical hydraulic conductivity Kv(Table 2) of intertidal sediments at the locations of some pair-wells (W2, W4, W6, W8 and W12). More details about the apparatus and method can be found in Li et al. (2010), Ma et al. (2015), and Wang et al. (2014).

Table 2.

Kv–values obtained by the method in Wang et al. (2014)


Seepage Face Development Condition

The presence/absence of seepage face development is defined by the characteristics of tidal actions and tidal flat face in the study area (Turner, 1995).

where, Σ is the seepage face parameter: if Σ>1, there will be seepage face development during low tide; if Σ⩽1, no seepage face development. The parameter RT denotes the tidal range [L]; T refers to the tidal period [T]; n denotes the effective porosity of sediment [−]; KH refers to the horizontal hydraulic conductivity [LT−1 ]; β is the beach slope [−].

In our study area, RT = 0.25–1.6 m, T = 12 h, KH = 1.1 × 10−5 m s−1(Table 2), β = 0.4‰, n = 0.03 – 0.35. Considering that the horizontal hydraulic conductivity is greater than or equal to the vertical one (KHKV) generally, a maximum value of the vertical hydraulic conductivity (1.1×10−5 m s−1) was used for conservative estimation. No measurements were conducted for the effective porosity of the silty beach sediments. According to Fetter (1994), the porosity of silt, silt sands, fine sands and sandy clay ranges from 0.03 to 0.35, which is used here for conservative estimation. According to these parameters and Equation (1), the seepage face parameter ranged from 108 to 1010, much greater than 1. Therefore, seepage faces will develop during low tides in the gentle silty tidal flat.

Direct Estimation of Water Exchange Rate

The observed data and preliminary analysis (Figure 4) indicated that the beach surface between W2 and W14 was saturated during the entire monitoring period. The generalized Darcy's law considering the density effect can be used to estimate the vertical flow rate by the following equation (Ma et al., 2015):

where qn is the vertical flow rate [L T−1 ]; KV is the vertical hydraulic conductivity [L T−1 ]; hup and hlow denote the freshwater-equivalent hydraulic head [L] at the up and low LTC–Divers, respectively; ΔL is the distance between the up and low LTC–Divers (52.8 cm); cup and clow denote the salinity [g l−1 ] of the pore water at the up and low LTC-Divers, respectively; δ is the ratio of the dynamic viscosity of freshwater to saltwater [−]:
where μ0 and μ denote the viscosity [ M L−1 T−1] of freshwater and saltwater, respectively; ξ; is a constant, equaling 1.566×10−3 m3 kg−1; c is the salinity of pore water [ M L−3 ], which is the average of salinities at up and low LTC-Divers, considering that the salinity at the same location did not change much.

The parameter ϵ is a constant, equaling 7.143×10−4 m3 kg−1, which is used to describe the linear relationship between density and salinity, which is defined by

where ρ and ρ0 denote the density [M L−3 ] of saltwater and freshwater, respectively.

Therefore, the generalized hydraulic gradient J is defined as:


The generalized hydraulic gradients J of the 12 pair-wells (W2–W8, W10–W14) were shown in Figure 5.

Figure 5.

The generalized hydraulic gradients of (a) W2–W5; (b) W6–W10 (except W9); (c) W11–W14. Observed tidal head (25%–gray line) are shown in the secondary axis


If J >0 (i. e., qn >0), the groundwater flows upward (SGD); otherwise, it flows downward (inflow). Normally, the density effect leads to downward flow due to the positive parameter ϵ.


As shown in Figure 2, the beach surface of W1 is located in the high intertidal zone, and the low LTC-Diver was inundated only during 7.92–9.67 d, 17.25–19.54 d, and 21.92–23.67 d (Figure 4a). In other words, the LTC – Divers installed at W1 were above the water table and did not able to record the ground-water data during most of the observation time. Besides, the salinities at up and low LTC-Divers were 0 since no groundwater entered the inner pipe when the LTC-Diver were above the water table. Thus, the seawater-groundwater exchange estimation later did not include W1.

As shown in Figure 4, the observed hydraulic heads at the up and low LTC-Divers of observation wells W2–W14 decreased with the falling tide when the tidal level was higher than the beach surface, and became almost constant when the tidal level fell near the beach surface until the next rising tide coming. Note that during the exposure of beach surface, the hydraulic heads of the two divers were higher than the beach surface. Furthermore, the hydraulic head of the low LTC-Diver was always higher than that of up LTC-Diver, especially at W2, W4 – W8, and W13. The generalized hydraulic gradient between the up and low LTC-Divers was positive (except W11) in most time which is about 0.2 at W7 and W8 (Figure 5b), indicating considerable groundwater discharge.

The salinity difference between the up and low LTC-Divers was very small for W2-W14, and was less than 2 g l−1 at W10 and W11 (Figures 4j and 4k). Note that as shown in Figure 4, the salinities at low LTC-Diver at W2, W4, W7 and W8 were 0, probably because the pipe was blocked by fine sediments, which exist extensively in the silty tidal flat, or there were bubbles around the LTC-Diver preventing water from entering. The sudden decrease of hydraulic heads at W3 – W5 (Figures 4c, 4d, and 4e) after 18 d may be caused by the pumping of groundwater for salt production nearby (Feng, Tan, and Liu, 2013; Yu, 2014).

At W9, the low LTC-Diver did not worked (Figure 4i), thus the seawater-groundwater exchange estimation at W9 was not included. As shown in Figure 4i, the hydraulic head at up LTC-Diver was about 0.15 m higher than the beach surface, indicating that the groundwater discharge from the seepage face developed here and might be significant.

Evaluation of SGD and Inflow

The exchange rate qn given by Equation (2) was averaged over the entire monitoring period for the outflow qout and inflow qin at each well using the following equations:

where t0 is the calculated initial time (t0 = 0 h was set at 0:00 on August 17th, 2014, 5–6 days after the installation of the pair-wells to eliminate the disturbance to groundwater table); te is the calculated end time (te = 640 h, at 15:00 on September 12th, 2014). Let W1 be the origin of the x–axis (X1 = 0), and the coordinates of W2, …, W13, W14 be X2, …, X13, X14, correspondingly (Table 3), then the total SGD and inflow could be evaluated by the integral along the x–axis, namely:

Table 3.

Exchange rate qin, qout of every well (except W1 and W9), and corresponding inflow and SGD between the adjacent wells estimated by the pair-wells method along the entire transect; and those (Inflow-s and SGD-s) obtained by the single-well method


The values of qout (or qin) between Xi and Xi+1 are determined by linear interpolations of their values at Xi and Xi+1(i = 2, …, 8, 10,…, 12, 13).

The estimated results of the seawater-groundwater exchange rate at each well are shown in Table 3. The calculated SGD along the entire transect was 164.00 m2 d−1 and inflow was 6.5 m2 d−1. As shown in Figure 6, the inflow concentrating between W2 and W4, accounted for about 64.6% of the total inflow. The SGD concentrated between W7 and W10 (Figure 6), accounting for about 80.4% of the total SGD, which is probably due to the high-permeability (10−5 m s−1 there, one or two orders of magnitudes greater than those of other pair-wells locations) nearby (Table 2). Under the rainfall condition, the SGD always increases with more inland freshwater recharge. The SGD and inflow are season-dependent. It would be different in the wet and dry season with varying inland freshwater recharge, and it still needs further research.

Figure 6.

Outflow rates (squares, volume of groundwater discharging in unit area per day (qout obtained by the pair-wells method (solid line) and qouts obtained by the single-well (gray dotted blue))) are shown in the primary axis and inflow rates (rhombus, volume of seawater infiltrating in unit area per day (qin obtained by the pair-wells method (solid line) and qins obtained by the single-well method (gray dotted blue))) are shown in the secondary axis averaged over the observation period between W2 to W14 (except W9)


When removing the observation periods of W3 to W5 influenced by groundwater pumping (i. e., when the water pressure equaled essentially zero), the SGD and inflow along the entire transect were 170.7 and 5.6 m2 d−1, i. e., the total SGD increased by 6.7 m2 d−1 (about 4.1%), and the total inflow decreased by 0.9 m2 d−1 (about 14.1%), respectively. Though groundwater pumping had little effect on the water exchange of the entire transect, it did have great impact on the zone between W3 and W5 (details will be given later).

The Impacts of Groundwater Pumping for Salt Production

As shown in Figures 4c–4e, the hydraulic heads and salinities of W3–W5 showed a sudden drop between 18.5 d ~ 20 d, even close to the elevations of LTC-Divers sometimes. This was most likely caused by groundwater pumping for salt production nearby. With pumping effect, the generalized hydraulic gradient changed from positive to negative suddenly, and gradually recovered to normal state after t = 20 d (Figure 5a). Meanwhile, the hydraulic heads and salinities recovered gradually within two days. This phenomenon was not observed at W2 and W6, indicating that pumping only had an impact on local area.

The SGD and inflow between W3 and W5 were 6.9 m2 d−1 and 46.7 m2 d−1 during the pumping period, respectively. The SGD and inflow were 12.3 and 0.6 m2 d−1 when ignoring the pumping effects. The SGD decreased by 5.4 m2 d−1, and the inflow increased by 46.1 m2 d−1 with the pumping effects. During the pumping period, the maximum inflow between W3 and W5 increases to 217.3 m2 d−1. These results indicated that short-term high-rate pumping had great influence on the interactions of seawater-groundwater between W3 and W5.


Comparison with the Single-well Method

We also estimated the SGD and inflow using the data from low LTC-Divers and tidal level. The difference between the two estimation methods is how to calculate the vertical groundwater flow rate qn. The vertical groundwater flow rate is calculated by the single-well method as following:

where, hw and cw represent the freshwater-equivalent hydraulic head [ L] and salinity [ g l−1 ] of the pore water at low LTC-Diver; csea denotes the seawater salinity [g l−1 ]; hsurface depends on the surface elevation of the location of LTC-Diver and the tidal head when tide coming; δ, KV and ϵ have the same meaning as those in Equation (2); ΔL is the depth of the low LTC-Diver from the beach surface [ L] (Figure 2 and Table 1). The parameter hsurface is defined as following:
where, Esurface is the beach surface elevation [L], Etide denotes the observed tidal level [L], htide refers to the freshwater-equivalent hydraulic head of the tidal level [L], and htide = (1 + ϵcsea) × Etide. When the beach surface is submerged by seawater, hsurface equals the freshwater-equivalent hydraulic head of the tidal level; otherwise, it equals the beach surface elevation during the beach surface exposure. More details about the single-well method are shown in Ma et al. (2015).

The calculated seawater-groundwater exchange rates using the single-well method are shown in Figure 6. As shown in Table 3, the estimated SGD–s and inflow–s (SGD–s and inflow–s, the SGD and inflow calculated using the single-well method) are 223.0 and 13.3 m2 d−1, much greater than those estimated by the pair-wells method. This might be related to the following aspects:

(1) The salinity difference between seawater and underground brine. As shown in Figure7, the underground brine salinities are between 30 and 40 g l−1, up to about 56 g l−1, which was much higher than the averaged seawater salinity about 23 g l−1. When the averaged salinity between seawater and groundwater was used to estimate the vertical groundwater flow rate qn, which will be less than the averaged salinity of pair-wells.

(2) The different tidal levels at different locations. The tidal level will decrease with the spread of tide landward (Lian and Zhao, 1995; Wang et al., 1995), due to the influence of friction of beach surface (Kamphuis, 1978), and the coefficient of friction increases with decreasing depth (Yu and Li, 2002; Zhang, Lv, and Sun, 2012). But the differences of the tidal level at different locations are neglected in the single-well method. What's more, the SGD and inflow estimation in the pair-wells method used the hydraulic head at the up LTC-Diver, which is more close to the real hydraulic head at the observation well locations than the tidal level or beach surface here used in the single-well method.

Based on these reasons, the estimated results and errors of the single-well method would be greater than those by the pair-wells method.

Density Effect and Salinity Distribution

The generalized Darcy's law used in both the pair-wells method and the single-well method above considers the density effect. The salinity along the entire transect is great, which will have significant impacts on the seawater-groundwater exchange.

Ignoring the density effect, the vertical groundwater flow rate qn can be written as:

where, KV is the vertical hydraulic conductivity of the sediment at observation well locations [ LT−1 ]; the parameters hup and hlow denote the freshwater-equivalent hydraulic head at up and low LTC–Divers [L]; ΔL is the distance between the up and low LTC–Divers, equaling 52.8 cm.

The estimation of SGD by Equation (12) increased by 25.6% to 206.0 m2 d−1, and inflow decreased by 27.4% to 4.7 m2 d−1, compared to those estimated by pair-wells method. Due to the high salinity of underground brine, the influence of density effect can change the exchange rate by nearly 25%, indicating that much downward groundwater water flow driven by high density brine.

Considering that the salinity at every LTC–Diver was relatively stable (Figure 4), the salinity observed at each LTC–Diver was averaged and the results were shown in Figure 7. One can see that along the entire transect the salinity decreased from 56.2 to 25.7 g l−1 seaward gradually, but all of them were higher than the seawater salinity, averaged 22.6 g l−1(Figure 7). The reasons for this might include the following aspects: (1) evaporation (the exposure time of the tidal flat increases landward, leading to a landward increasing-evaporation and in turn a landward increasing of salinity), (2) human activities such as exploitation of the underground brine and construction of many salt pans, as shown in Figure 1c (such activities often discharge waste salt brines into the coastal areas, leading to the high salinity there), and (3) the inland fresh groundwater discharge near the low tide zone (Heiss and Micheal, 2014; Li, Boufadel, and Weaver, 2008; Robinson et al., 2007), which may lead to the decrease of salinity of the groundwater in the lower intertidal zone.

Figure 7.

Observed salinity at the up LTC–Diver (solid dots) and low LTC–Diver (hollow dots) of W2, …, W13, W14 and the seawater (hollow dots)


Comparison with Ma et al. (2015)

The estimated SGD by the pair-wells method in this study was 164.0 m2 d−1, much greater than the SGD about 8.8 m2 d−1 calculated in Ma et al. (2015); while the inflow along this transect was 6.5 m2 d−1, less than that nearly 15.3 m2 d−1 in Ma et al. (2015) (Table 4). The differences in SGD and inflow between the two results are significant since there are many differences between the two transects.

Table 4.

The comparison of the two typical transects in this paper and Ma et al. (2015)


The length of the cross-shore transect in this paper is 3514 m, which is about nine times long of the transect (392 m) in Ma et al. (2015) (Table 4). There were six single-well along the entire transect in Ma et al. (2015), while 14 pair-wells were set along our transect in this paper (Figures 1 and 2). The selected cross-shore transect in Ma et al. (2015) has two distinct slopes, with the landward part of a relatively steep slope about 2%, and the seaward part of a gentle slope about 0.1%. The slope of the selected transect in this paper is about 0.4 ‰, which is much smaller than the slopes in Ma et al. (2015) (Table 4).

Field measurements were conducted from 11:00 September 9th to 15:00 September 14th, 2012, about 125 hours in Ma et al. (2015) (Table 4). The observation time from 20:00 August 10th to 23:00 September 12th, 2014, almost 800 hours in this paper, is very long compared with most previous studies (Table 4).

The selected transect in this paper is heterogeneous along the cross-shore direction. The permeability ranges from 10−7–10−5 m s−1, which is one or two orders of magnitude less than that about 10−5 m s−1 in Ma et al. (2015) (Table 4). There are little water flowing downward due to the low-permeability surface sediments and the higher hydraulic head of the groundwater when surface was submerged.

The salinity of groundwater along the transect in Ma et al. (2015) ranged from 16 to 27 g l−1, which is less than the salinity of seawater. The transect in this study area is located in brine zones, and the salinity of groundwater at observation locations ranged from 25 to 56 g l−1, which is greater than the seawater salinity averaged about 23 g l−1(Figure 7). The details of the salinity distribution have been discussed in the “Density effect and salinity distribution” section.


This paper investigated the seawater-groundwater exchange in a silty tidal flat with a large-scale seepage face in the south coastal of Laizhou Bay. We selected 14 locations along a typical transect in the intertidal zone with the slope of 0.4‰ and the cross-shore width of 3514 m, to install observation pair-wells monitoring the groundwater head, salinity and temperature at two different depths simultaneously. The vertical hydraulic conductivity ranges from 5.4×10−7 to 1.1×10−5 m s−1, and the salinity of the underground brine is from 25 to 56 g l−1 along the selected transect. Basing on the observed data from 20:00 August 10th to 23:00 September 12th, 2014 and generalized Darcy's law, the submarine groundwater discharge (SGD) and inflow along the studied transect were estimated to be 164.0 and 6.5 m2 d−1, respectively; within 80.4% of the total SGD concentrated between W7 and W10 due to the high-permeability here, and 64.6% of the total inflow between W2 and W4.

The comparative study between the pair-wells method and single-well method indicated that the pair-wells method is more suitable for the long transect with a very gentle slope than the single-well method, due to the pair-wells can monitor the groundwater head, salinity and electrical conductivity at two different elevations simultaneously. The fixed distance of the pair-wells can avoid the uncertain errors of the single-well method.

The observed data indicated that pumping of groundwater has significant impacts on the exchange between seawater and ground-water. During the groundwater pumping period, the outflow decreased by 5.4 m2 d−1 and the inflow increased by 46.1 m2 d−1, even up to 217.3 m2 d−1 within the tidal zone from W3 to W5.

The density effect had great impacts on the seawater-groundwater exchange in the brine zones due to its high salinity ranging from 25 to 56g l−1 along the entire transect. Leaving out the density effect, the exchange rate can change 25%, indicating that much downward groundwater flow driven by the density effect.

The salinity decreased seaward gradually along the entire transect, probably because the evaporation, human activities such as exploitation of the underground brine and construction of many salt pans, and the inland fresh groundwater discharge near the low tide zone.


This work was supported by the Key Program of National Natural Science Foundation of China (Grant No. 41430641), the National Basic Research Program of China (“973” Program, Grant No. 2015CB452902). We would like to thank Kai Xiao, Zongzhong Song, Shuaishuai Liu, Yan Zhang and Shaohong Li for their field work.



Y. F. Bi H. J. Yu X. Y. Xu Q. Su 2012. Study on the groundwater hydrochemical characteristics in the southern Laizhou Bay. Marie Science Bulletin, 31(3), 241–248. Google Scholar


J. E. Cable W. C. Burnett J. P. Chanton D. R. Corbett P. H. Cable 1997. Field Evaluation of Seepage Meters in the Coastal Marine Environment. Estuarine, Coastal and Shelf Science, 45, 367–375. Google Scholar


R. A. Cox F. Culkin J. P. Riley 1967. The electrical conductivity/chlorinity relationship in natural sea water. Deep Sea Research and Oceanographic Abstracts, 14 (2), 203–220. Google Scholar


S. T. Feng X. F. Tan G. Liu 2013. Potentiality Analysis on Underground Brine Resources in Shandong Province. Shandong Land and Resources, 29, 69–73. Google Scholar


C. W. Fetter 1994. Applied Hydrogeology. New Jersey: Prentice–Hall. Google Scholar


H. Q. Gu W. S. Moore L. Zhang J. Z. Du J. Zhang 2012. Using radium isotopes to estimate the residence time and the contribution of submarine groundwater discharge (SGD) in the Changjiang effluent plume, East China Sea. Continental Shelf Research, 35, 95–107. Google Scholar


Q. N. Guo H. L. Li M. C. Boufadel Y. Sharifi 2010. Hydrodynamics in a gravel beach and its impact on the Exxon Valdez oil. Journal of Geophysical Research, 115, C12077, doi: 10.1029/2010JC006169. Google Scholar


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


F. Han J. C. Wu Y. X. Zhang 2011. Numerical simulation for the salt/brine water intrusion in a leaky aquifer: A case study of the southern beach of Laizhou Bay, eastern China. Journal of Nanjing University, 47(3), 281–290. Google Scholar


J. W. Heiss H. A. Michael 2014. Saltwater-freshwater mixing dynamics in a sandy beach aquifer over tidal, spring-neap, and seasonal cycles. Water Resources Research, 50(8), 6747–6766. Google Scholar


J. W. Kamphuis 1978. Attenuation of gravity waves by bottom friction. Coastal Engineering, 2, 111–118. Google Scholar


C. M. Lee J. J. Jiao X. Luo W. S. Moore 2012. Estimation of submarine groundwater discharge and associated nutrient fluxes in Tolo Harbour, Hong Kong. Science of the Total Environment, 433, 427–433. Google Scholar


D. G. Li M. H. Zhao M. Han A. X. Jiang Z. L. Zhang 2000. A study of the shallowly-buried paleochannel zones in the south coast plain of the Laizhou Bay. Marine Geology & Quaternary Geology, 20 (1), 23–29. Google Scholar


H. L. Li M. C. Boufadel 2010. Long-term persistence of oil from the Exxon Valdez spill in two-layer beaches. Nature Geoscience, 3(2), 96–99. Google Scholar


H. L. Li M. C. Boufadel J. W. Weaver 2008. Tide-induced seawater-groundwater circulation in shallow beach aquifers. Journal of Hydrology, 352, 211–244. Google Scholar


H. L. Li J. J. Jiao 2013. Quantifying tidal contribution to submarine groundwater discharges: a review. Chinese Science Bulletin, 58 (25), 3053–3059. Google Scholar


H. L. Li P. P. Sun S. Chen Y. Q. Xia S. Liu 2010. A falling-head method for measuring intertidal sediment hydraulic conductivity. Ground Water, 48(2), 206–211. Google Scholar


L. Li D. A. Barry F. Stagnitti J. Y. Parlange 1999. Submarine groundwater discharge and associated chemical input to a coastal sea. Water Resource Research, 35 (11), 3253–3259. Google Scholar


J. J. Lian Z. D. Zhao 1995. Interactions between wave-current and muddy sea-bed. Journal of Sediment Research, 3, 40–50. Google Scholar


Q. Ma H. L. Li X. J. Wang C. Y. Wang L. Wan X. S. Wang X. W. Jiang 2015. Estimation of seawater-groundwater exchange rate: case study in a tidal flat with a large-scale seepage face (Laizhou Bay, China). Hydrogeology Journal, 23(2), 265–275. Google Scholar


W. S. Moore 1996. Large groundwater inputs to coastal waters revealed by 226Ra enrichments. Nature, 380, 612–614. Google Scholar


W. S. Moore 2007. Estimating submarine groundwater discharge in the range 101 to 108 km2 using Ra isotopes. Geochimica Et Cosmochimica Acta, 71(15), A684–A684. Google Scholar


W. S. Moore J. D. Oliveira 2008. Determination of residence time and mixing processes of the Ubatuba, Brazil, inner shelf waters using natural Ra isotopes. Estuarine Coastal and Science, 76(3), 512–521. Google Scholar


R. J. Paulsen C. F. Smith D. O'Rourke T. F. Wong 2001. Development and Evaluation of an Ultrasonic Ground Water Seepage Meter. Groundwater, 39(6), 904–911. Google Scholar


T. Riedel K. Lettmann M. Beck H. J. Brumsack 2010. Tidal variations in groundwater storage and associated discharge from an intertidal coastal aquifer. Journal of Geophysical Research, 115, C04013, doi: 10.1029/2009JC005544. Google Scholar


C. Robinson B. Gibbes H. Carey L. Li 2007. Salt-freshwater dynamics in a subterranean estuary over a spring-neap tidal cycle. Journal of Geophysical Research-oceans, C9 (112), C09007, doi: 10.1029/2006JC003888. Google Scholar


C. Robinson L. Li D. A. Barry 2007. Effect of tidal forcing on a subterranean estuary. Advances in Water Resources, 30(4), 851–865. Google Scholar


C. P. Slomp P. V. Cappellen 2004. Nutrient inputs to the coastal ocean through submarine groundwater discharge: controls and potential impact. Journal of Hydrology, 295 (1–4), 64–86. Google Scholar


M. Taniguchi Y. Fukuo 1993. Continuous measurements of ground-water seepage using an automatic seepage meter. Ground Water, 31(4), 675–679. Google Scholar


I. L. Turner 1995. Simulating the influence of groundwater seepage on sediment transported by the sweep of the swash zone across macro-tidal beaches. Marine Geology, 125(1–2), 153–174. Google Scholar


L. Wang Y. Jin Y. Li J. Ding 1995. Approact to Characteristics of Wave Breaking. The Ocean Engineering, 13 (3), 81–89. Google Scholar


X. J. Wang H. L. Li J. J. Jiao D. A. Barry L. Li X. Luo C. Y. Wang L. Wan X. S. Wang X. W. Jiang Q. Ma W. J. Qu 2015. Submarine fresh groundwater discharge into Laizhou Bay comparable to the Yellow River flux. Science Reports. doi:10.1038/srep08814. Google Scholar


X. J. Wang H. L. Li J. Z. Yang L. Wan X. S. Wang X. W. Jiang H. M. Guo 2014. Measuring in-situ vertical hydraulic conductivity in tidal environments. Advances in Water Resources, 70, 118–130. Google Scholar


W. S. Wooster A. J. Lee G. Dietrich 1969. Redefinition of salinity. Limnology and Oceanography, 14 (3), 437–438. Google Scholar


Y. Q. Xia H. L. Li 2012. A combined field and modeling study of groundwater flow in a tidal marsh. Hydrology and Earth System Science, 16(3), 741–759. Google Scholar


Y. Q. Xia H. L. Li M. C. Boufadel Y. Sharifi 2010. Hydrodynamic factors affecting the persistence of the Exxon Valdez oil in a shallow bedrock beach. Water Resources Research, 46(10), W10528, doi: 10.1029/2010WR009179. Google Scholar


B. C. Xu D. Xia W. C. Burnett N. T. Dimova H. J. Wang L. J. Zhang M. S. Gao X. Y. Jiang Z. G. Yu 2014. Natural 222Rn and 220Rn indicate the impact of the Water-Sediment Regulation Scheme (WSRS) on submarine groundwater discharge in the Yellow River estuary, China. Applied Geochemistry, 51, 79–85. Google Scholar


Y. Q. Xue J. C. Wu S. J. Ye Y. X. Zhang 2000. Hydrogeological and Hydrogeochemical Studies for Salt Water Intrusion on the South Coast of Laizhou Bay, China. Ground Water, 38 (1): 38–45. Google Scholar


Y. Yu Y. C. Li 2002. Dissipation of wave energy on very mild slope. Journal of Hydrodynamics, 17 (5), 553–564. Google Scholar


Y. H. Yu 2014. The development of salt industry on south of Laizhou Bay in northern Shandong. Salt Industry History Research, 2, 8–16. Google Scholar


J. C. Zhang X. Q. Lv L. Y. Sun 2012. Study on bottom friction parameterizations for two-dimensional tidal models using data assimilation. Chinese Journal of Computational Mechanics, 29(1), 49–54. Google Scholar


W. H. Zhang Y. Q. Cai C. C. Wu 2002. A research into some environmental aspects of beach erosion. Journal of Safety and Environment, 2 (3), 8–12. Google Scholar


Z. L. Zhang L. M. Peng 1998. The groundwater hydrochemical characteristics on seawater intruded in eastern and southern coasts of Laizhou Bay. China Environmental Science, 18 (2), 121–125. Google Scholar
Lijia Hou, Hailong Li, Chunmiao Zheng, Qian Ma, Chaoyue Wang, Xuejing Wang, and Wenjing Qu "Seawater-groundwater Exchange in a Silty Tidal Flat in the South Coast of Laizhou Bay, China," Journal of Coastal Research 74(sp1), 136-148, (1 March 2016).
Received: 23 February 2015; Accepted: 12 July 2015; Published: 1 March 2016
Laizhou Bay.
seawater-groundwater exchange
seepage face
silty beach
Back to Top