Abstract

Although there are many existing analytical studies of tidal groundwater level fluctuations in coastal aquifer systems, few of them focus on an offshore submarine aquifer. Here, we consider tidal groundwater head fluctuations in a submarine leaky confined aquifer overlain by a semipermeable seabed. Both the seabed and the confined aquifer are assumed to extend horizontally infinitely. A one-dimensional mathematical model is established to describe the problem, and the analytical solution is derived. The impacts of the tidal loading efficiency, hydraulic conductivity and elastic storage of the semipermeable layer and aquifer on the groundwater head fluctuations in the aquifer system are analyzed and discussed. Solution analyses indicated that tidal loading effects tend to enhance the amplitude of the tidal groundwater fluctuation in the confined aquifer system and to reduce the phase shift between the groundwater head and the sea tide fluctuations.

1. Introduction

The study of tidal groundwater level fluctuations in coastal aquifer systems is beneficial to solving various coast environmental problems. Many analytical solutions have been derived in this field since the 1950s. Jacob [1] and Ferris [2] firstly derived analytical solutions to describe the tide-induced groundwater fluctuation in a confined aquifer which ends at the coastline. The equations have been widely used in field survey and aquifer parameter estimation in coastal areas (e.g., [39]). Since then, various analytical solutions related to tidal head fluctuations have been derived for coastal confined aquifers. Van der Kamp [10] investigated the tidal groundwater fluctuation in an aquifer with an impermeable roof extending seaward and landward infinitely. Li and Chen [11] considered a coastal confined aquifer with an impermeable roof extending under the sea for a certain length, neglecting the vertical leakage. Jiao and Tang [12] discussed the tide-induced groundwater head fluctuation in an onshore multilayered coastal aquifer system neglecting the elastic storage of the semipermeable layer. Based on Jiao and Tang’s study [12], Li and Jiao [13] took into account the elastic storage of the semipermeable layer, and they [14] considered a leaky aquifer with the roof extending under the sea. Jeng et al. [15] considered the interference of the tidal groundwater level fluctuations in a multilayered aquifer system with the leakage of a semipermeable layer considered. Li and Jiao [16] improved Jeng et al.’s work [15] in the sense that both the leakage and elastic storage of the leaky layer are considered. Li et al. [17] considered a coastal multilayered aquifer system with the confined aquifer and its semipermeable roof extending infinitely under the sea. All the effects related to the semipermeable layer such as the elastic storage, the leakage, and the tidal loading efficiency were taken into account. Guarracino et al. [18] derived an analytical solution to describe the tide-induced head fluctuation in an aquifer system that extends a finite distance under the sea, with which the extending length may be estimated from groundwater head measurement. Wang et al. [19] presented an analytical study on the loading effect of the water table variations in a coastal aquifer system consisting of an unconfined aquifer, a confined aquifer, and an impermeable layer between them.

Due to the tidal fluctuations, the total weight of the pore water (the onshore part) or sea water (the offshore part) in/overlying the aquifers will vary accordingly, which in turn can cause a considerable tidal loading effect on the groundwater flow. Many researchers [11, 14, 1720] considered the tidal loading effect in their analytical studies, whereas almost all of the existing studies of tide-induced groundwater fluctuations focus on the tide-induced groundwater level fluctuation in onshore aquifer. Although there are some analytical solutions of tide-induced groundwater fluctuations in the submarine confined aquifer system in the offshore part (e.g., [14, 17, 18]), the dependency of the tide-induced groundwater head fluctuation in the offshore submarine aquifer on the model parameter was not analyzed or discussed. The offshore submarine aquifers are important both in water resources and in marine environments (e.g., [2126]). Therefore, the study of the groundwater flow in offshore submarine aquifer systems is necessary and important.

This paper considered a submarine leaky confined aquifer system with the confined aquifer and its semipermeable roof extending infinitely under the sea. All effects related to the submarine leaky confined aquifer system such as the elastic storage, the hydraulic conductivity, and the tidal loading efficiency were taken into account. The mathematic model was given, and the analytical solution was derived. The influences of the tidal loading efficiency, the elastic storage, and the hydraulic conductivity on the tidal groundwater head fluctuations in the confined aquifer and in the seabed were analyzed and discussed.

2. Mathematical Model and Analytical Solution

2.1. Mathematical Model

The configuration of a submarine leaky confined aquifer system is shown in Figure 1. The aquifer is overlain by a semipermeable seabed aquifer/layer. Both the semipermeable seabed and the confined aquifer are assumed to extend horizontally infinitely. This assumption is reasonable when the considered offshore submarine aquifer is far enough from the coastline. There are even many offshore freshwater submarine aquifers extending in all horizontal directions extensively and very far from the coastline [2226]. In addition, several existing analytical studies assume that the coastal aquifer systems extend seaward infinitely (e.g., [10, 14, 17, 27]). A vertically one-dimensional mathematical model is used to describe the problem.

Let the axis be positive upwards (Figure 1). All the layers are homogeneous and with constant thickness. Based on these assumptions and the theories of leaky, elastic aquifers proposed by Hantush and Jacob [28] and Jacob [1], the mathematical model to describe the groundwater fluctuations in Figure 1 is as follows: where , , and are the hydraulic head [L], specific storage [L-1], and vertical hydraulic conductivity [LT-1] of the semipermeable layer or the aquifer, respectively; is the tidal loading efficiency (dimensionless) [‒‒]; is the density ratio (dimensionless) of the seawater to the freshwater with the range from 1.02 to 1.03; and are the thicknesses [L] of the aquifer and the semipermeable layer, respectively; takes the value of in the semipermeable layer and in the confined aquifer; equals in the semipermeable layer and in the confined aquifer; is the hydraulic head of the sea tide [L]; is the amplitude [L] of the tidal fluctuation; is the angular velocity (or frequency) [T-1] of tide and equals with being the tidal period [T]; and is the phase shift (in radian). The bottom of the confined aquifer is set to be the datum of the system.

Equation (1) indicates that the groundwater level fluctuation in the submarine leaky confined aquifer is caused by its elastic compression and expansion due to the tidal loading rate of the sea water above seabed (), in addition to the tidal fluctuation at the seabed boundary as expressed by equation (2). Equations (3) and (4) are the continuity conditions of the hydraulic head and groundwater flux, respectively. Equation (5) states the no-flow boundary condition on the upper surface of the impermeable bottom.

2.2. Analytical Solution

In order to facilitate analysis, four parameters are introduced: the main aquifer’s tidal propagation parameter [L-1] defined as the dimensionless thickness of the seabed defined as where is the dimensionless storativity ratio and defined as

and the dimensionless leakage defined as

The solution of boundary value problem (1), (2), (3), (4), and (5) is given by the following equations (see Appendix A for the derivation):

The analytical solutions (10a) and (10b) of the groundwater heads can be deformed as (see Appendix B for the derivation of the analytical solutions) where is the amplitude ratio and is the comprehensive tidal phase shift in the aquifer defined as

The three parameters , , and are given by the following equations (see Appendix C for the derivation): where is the thickness ratio (dimensionless) and and are two complex constants determined by equations (A.4) and (A.5) and given by with

3. Discussion of the Solution

3.1. Value Ranges of Model Parameters

The value ranges of dimensionless model parameters , , , and for the leaky confined aquifer system reported in the literature are shown in Table 1, according to the summarization of Li and Jiao ([13], Table 1). One can see that ranges from 0.014 to 0.055, from 1.37 to 25.87, from to , and from 1.22 to 3.50. Based on these data and other situations (for example, the thickness of the submarine confined aquifer could be up to several hundred meters), the ranges of the parameters in the following discussion will be chosen: 0.001~10 for , 0~100 for , ~ for , and 1~10 for .

The tidal loading efficiency () in the semipermeable roof of the coastal confined aquifer system is usually very close to 1 [17]. So the value of will be fixed at 0.9 in the following discussion. Previous analytical studies and field surveys [7, 29, 30] show that the range of the tidal loading efficiency () in the confined aquifer changes from 0.05 to 1.

3.2. Comparison with Existing Solutions
3.2.1. Case 1: Seabed without Elastic Storage ()

When the elastic storage and loading efficiency of the seabed are close to zero (; ), a simpler analytical solution of the confined aquifer system can be obtained. In view of equation (7), one has . For the confined aquifer, letting and in equation (13a) yields

Letting in equations (12a) and (12c) yields

Assuming in equations (14b) and (14c) yields

and then yields

Substituting equations (14f) and (14g) back into equation (11), the groundwater head in the confined aquifer can be written as which is essentially the same as the equation given by Li and Jiao [14] when the roof length and . Therefore, our analytical solution is a generalization of Li and Jiao’s solution [14] when the roof length and .

For the semipermeable seabed, letting and in equations (12b) and (12d) yields

Substituting equations (16a) and (16b) back into equation (11), the groundwater head in the semipermeable seabed can be obtained as follows:

3.2.2. Case 2: Impermeable Seabed ()

When the semipermeable seabed becomes impermeable (i.e., ), in view of equation (9), one has . Letting in equation (13a) yields

Letting in equations (12a) and (12c) yields

Substituting equations (18b) and (18c) back into equation (11), the groundwater head in the confined aquifer can be written as

The solution given by equation (19) is essentially the same as the solution of Van der Kamp [10] and Li et al. [17] when and except some notation differences.

3.3. The Amplitude Ratio and Phase Shift

The amplitude ratio and phase shift of the groundwater head fluctuation in the confined aquifer and seabed are given by equations (12a), (12b), (12c), and (12d). Here, we will concentrate on the dependency of and on the model parameters such as the tidal loading effect of the confined aquifer () and of the seabed () and the dimensionless thickness of the semipermeable seabed and of the confined aquifer. Since the amplitude ratio and the phase shift depend on , we discuss the changes of and when is fixed at the middle of the confined aquifer () and at the middle of the semipermeable seabed (). The representative values of and in the literature are chosen for the following discussion: i.e., for and 3 for (see Table 1).

3.3.1. Changes of and in the Confined Aquifer

Figure 2 shows how the amplitude ratio at the middle of the confined aquifer changes with for different values of and when . From Figure 2, one can see that the amplitude ratio decreases from a constant close to 1 to (which equals 0.8, 0.5, and 0.1, respectively) as increases from −3 to 1. When is small (), with the θ increasing, decreases from 1 firstly and then increases to 0.9, which is the value of the tidal loading efficiency () of the semipermeable seabed. When , approaches the fixed values of and becomes independent of . For intermediate values of (approximately in the range ), decreases significantly with . The curves for different values may cross each other near , indicating that is not always monotonic with respect to . In addition, also increases with for fixed values of and . This can be observed when one compares the curves for the same values of and .

Figure 3 shows how the phase shift at the middle of the confined aquifer changes with for different values of and when . From Figure 3, one can see that the value range of phase shift decreases significantly with the tidal efficiency . The value range of is about (0, 0.96) for , (0, 0.34) for , and (0, 0.11) for . When ab is very small (), with the increasing of , increases from 0 firstly and then decreases to 0. When is very large (), the phase shift tends to be zero for any values of and . From Figure 3, one can also see that is nonmonotonic with respect to and , which is due to the superposition of the multieffects such as the tidal loading, the downward propagation of the tidal signal, and the elastic storages of the semipermeable seabed and the aquifer. According to the curve when and , the maximum phase shift is close to 1.0 (rad) which can lead to a time lag about 2 h for a semidiurnal tide and 4 h for a diurnal tide.

According to the discussions about the amplitude ratio and the phase shift in the confined aquifer, one can conclude that tidal loading effects tend to enhance the amplitude and reduce the phase shift of the tidal groundwater head fluctuations in the confined aquifer. Figures 2 and 3 show that the increase of the dimensionless thickness of the upper semipermeable seabed may significantly enhance the tidal head fluctuation in the confined aquifer. Furthermore, Figures 2 and 3 demonstrate that an active range of () exists in which both the amplitude and phase shift change with significantly and both become insensitive to outside the range (-2.5, 1).

3.3.2. Changes of and in the Seabed

Figure 4 shows how the amplitude ratio at the middle of the seabed semipermeable layer changes with for different values of and when . It is obviously that no matter how the and change, the value of becomes constant when . However, the changes of become complex when . When is very small (), as is increasing from 0 to 10, decreases from 1 firstly and then increases to 0.9, which is the value of the tidal loading efficiency of the semipermeable seabed. As is increasing, the trends of are different with values. When is small (), decreases firstly and then becomes a constant independent of . However, when is in the range , increases firstly and then becomes a constant independent of . When is around 3.0, decreases firstly and then increases with and finally becomes a constant. The curves for different values cross each other near , indicating that the amplitude ratio of the seabed is not always monotonic with respect to . Comparing Figures 4(a)4(c), one can find that when is larger, tends to be a constant value of 0.9 in Figure 4(a); however, tends to have different values shown in Figures 4(b) and 4(c). It may be caused due to the tidal loading effect in the confined aquifer.

Figure 5 shows how the phase shift at the middle of the seabed semipermeable layer changes with for different values of and when . The value range of changes significantly with , which is about (0.00, 0.05) for , (–0.11, 0.14) for , and (–0.28, 0.29) for . Thus, the smaller the tidal efficiency , the wider the value range of the phase shift. In general, the tidal loading effect tends to reduce the phase shift between the groundwater head and the tide fluctuations. No matter how and change, the value of becomes constant when . When is equal to 0 or 1, the phase shift increases firstly and then decreases to a constant with the increasing of . For , the phase shift decreases with monotonically to a constant. When is equal to or greater than 7, the phase shift increases with monotonically to a constant. Comparing Figures 5(a)5(c), one can also find that when is larger, tends to be a constant value of 0.9 in shown Figure 5(a); however, tends to have different values shown in Figures 5(b) and 5(c). It may be also caused due to the tidal loading effect in the confined aquifer.

3.4. Effect of the Thickness Ratio and Hydraulic Conductivity Ratio

Figure 6 shows how the amplitude ratio and phase shift φ at the middle of the confined aquifer change with for different values of the thickness ratio . The value range of used is from 1 to 10, which is based on Table 2. Figure 6(a) demonstrates that when increases from -3 to 1, the amplitude ratio decreases from to for any values of . Given any values of in the range of (-2, 0.4), the amplitude ratio increases with . Outside the range of (-2, 0.4), the amplitude ratio becomes insensitive to . In Figure 6(b), the phase shift is a monopeak function of for any values of . The increasing of the thickness ratio can decrease the phase shift when and increase the phase shift when . The impacts are complex in other ranges. Overall, the phase shift ranges from 0 to around 0.12, a relatively narrow value range.

Figure 7 shows how the amplitude ratio and phase shift φ at the middle of the semipermeable seabed layer change with for different values of . Figure 7(a) shows that the amplitude ratio decreases with first and then tends to be a constant with a slight oscillations. The increasing of the thickness ratio can lead to the increase of the amplitude for any value of in the range . Figure 7(b) shows that the phase shift decreases with monotonically from about 0.035 to about -0.11, a very narrow value range. The increasing of the thickness ratio can lead to the increase of the phase shift for any value of in the range .

Figure 8 shows how the changes of the amplitude ratio and phase shift at the middle of the confined aquifer with for different values of the hydraulic conductivity ratio , which ranges from 10-5 to 10-2 according to Table 1. Figure 8(a) shows that the amplitude ratio decreases with monotonically from 0.9 () to 0.5 () for any fixed value of . The increasing of can lead to the increase of the amplitude ratio for any value of in the range . Figure 8(b) shows that the phase shift is a monopeak function of for any given value of . The values of where the curve peaks increase with . The value range of the phase shift is from 0 to less than 0.15.

Figure 9 shows how the amplitude ratio and phase shift at the middle of the semipermeable seabed layer change with for different values of . Figure 9(a) shows that the amplitude ratio decreases with first and then tends to be a constant with slight oscillations. The increasing of hydraulic conductivity ratio can lead to the increase of the amplitude for any value of in the range . Figure 9(b) shows that the phase shift decreases with monotonically from about 0.035 to about -0.11, a very narrow value range. The increasing of hydraulic conductivity ratio can lead to the increase of the phase shift for any value of in the range .

3.5. Groundwater Head at Different Depths of the Confined Aquifer

In order to investigate the differences of groundwater head fluctuations at different depths of the confined aquifer, we considered the changes of and at the top and the bottom of the confined aquifer. The parameters and are fixed in the following discussion: for . Figure 10 shows how the amplitude ratio at the top and bottom of the confined aquifer changes with for different values of and when . The variation patterns of at the top and bottom of the aquifer are similar to that at the middle of the confined aquifer (Figure 2). When , Figure 10(a) shows that the curves of at the top and bottom of the aquifer coincide completely with each other, and so do the curves for the phase shift in Figure 10(b). Figures 10(c) and 10(e) show that when , the curves of at the top and bottom of the aquifer are also entirely coincident; however, at the top of the aquifer becomes larger than that at the bottom of the aquifer with when . The difference in the value of between the top and bottom of the aquifer decreases gradually as increases. Figures 10(d) and 10(f) show that the phase shift at the top of the aquifer becomes smaller than that at the bottom of the aquifer for large values of (10 and 100). The difference in the value of between the bottom and top of the aquifer decreases gradually as increases. It can be concluded that the discrepancy between the tidal groundwater fluctuations at the top and bottom of the confined aquifer becomes larger with the increasing of and the decreasing of .

4. Conclusion

This paper investigated the tidal groundwater head fluctuation in a submarine confined leaky aquifer system. A vertical one-dimensional mathematical model is given with the analytical solution derived. Both the leakage and elastic storage of the semipermeable seabed are considered in this model. The derived analytical solution is different from previous analytical solutions by Li and Jiao [14] and Guarracino et al. [18] which focus on the analysis in the inland portion. Although Li et al. [17] give the analytical solutions in the offshore part, they did not do sensitivity analysis of parameters of the aquifers in the offshore part to head fluctuations of the aquifers there. Our solution quantifies the tidal groundwater head fluctuations in an offshore submarine aquifer system, which contains a confined aquifer and a semipermeable seabed.

The effects of several dimensionless model parameters, including the tidal loading efficiencies ( and ) of the confined aquifer and the semipermeable seabed, the dimensionless thickness of the confined aquifer, and the dimensionless thickness θ of the seabed, on the tidal groundwater fluctuations are quantitatively discussed and analyzed. The following conclusions can be drawn based on the solution analyses and discussion: (1) The tidal loading effect tends to enhance the amplitude of the tidal groundwater fluctuation in the confined aquifer system, and it can also tend to reduce the phase difference between the fluctuation of the groundwater head in the aquifer (or an aquitard) and the tidal fluctuations. (2) The increase of the thickness ratio or the hydraulic conductivity ratio can not only result in a larger amplitude in the confined aquifer but also lead to a larger phase shift. (3) The discrepancy between the tidal groundwater fluctuations at the top and bottom of the confined aquifer increases when increases and/or decreases.

Appendix

A. Derivation of Solutions (10a) and (10b)

Let be a complex function of the real variables and that satisfies equations (1), (2), (3), (4), and (5) after in equations (1) and (2) is replaced by , where .

Let be the solution to (1), (2), (3), (4), and (5); it follows that where denotes the real part of the followed complex expression. Now suppose where is an unknown function of . Substituting (A.2) into the five equations which satisfies and extending the five resultant real equations into complex ones with respect to the unknown complex function , then dividing the results by , yield

The solution to (A.3), (A.4), (A.5), (A.6), (A.7), and (A.8) is where and are the dimensionless model parameters given by equations (6) and (7), respectively, and , , and are the complex constants given by equations (13a), (13b), (13c), (13d), (13e), and (13f).

B. Derivation of Solution (11)

The following is the derivation deformation of the analytical solutions (10a) and (10b).

C. Derivation of Solutions (13a), (13c), and (13d)

Putting equation (A.9) into equation (A.6) yields and then, the following equation can be obtained from (C.1):

Getting the derivative of equation (A.9) and then putting them into equation (A.6) yield

Putting equation (A.9) into equation (A.5) yields

The solutions to (C.2), (C.3), and (C.4) are with

Data Availability

The data used to support the findings of this study are included within the article.

Disclosure

This article was modified based on the master thesis of Mr. Zongzhong Song.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the Key Program of the National Natural Science Foundation of China (Grant No. 41430641) and the National Basic Research Program of China (“973” Program, Grant No. 2015CB452902). This work was also sponsored by Guangdong Provincial Key Laboratory of Soil and Groundwater Pollution Control (No. 2017B030301012) and State Environmental Protection Key Laboratory of Integrated Surface Water-Groundwater Pollution Control.