Exact simulation of the groundwater flow in anisotropic confined aquifer near time-varying streams by comprehensive new analytical solutions

This paper presents some new analytical solutions to an accurate prediction of the behavior of the groundwater flow in aquifers in response to changes in surface water. The new analytical solutions are obtained using integral transforms. An anisotropic rectangular confined aquifer bounded with four time-varying streams is undertaken. The effects of anisotropy on groundwater head and flow rate near time-varying streams are investigated. Depending on the change rates of the streams level, an anisotropic aquifer may render either lower or higher hydraulic head than an isotropic aquifer. In addition, an anisotropic aquifer has provided less water exchange at the interfaces than an isotropic one. The sensitivity of the hydraulic head to change rate of the streams level in both isotropic and anisotropic aquifers is evaluated. It is shown that the aquifer response is more sensitive to change rate of the streams parallel to y-direction and less sensitive to change rate of the streams parallel to x-direction in an anisotropic aquifer and vice versa in an isotropic aquifer. The results of the present new analytical solutions are compared with numerical model of MODFLOW. The results obtained from the presented solutions showed good agreement with the results of MODFLOW. The results show that the presented new analytical solutions are accurate, robust and efficient. Therefore, the results indicate that the presented new analytical solutions are very effective in the simulation of the groundwater flow in river–aquifer systems. Furthermore, one of the advantages of the new analytical solutions is to investigate the sensitivity analysis of aquifer parameters, which has been carried out in this paper. Also, some other new analytical solutions for steady-state conditions and sudden fall in streams level are provided as well. Feasibility of the proposed new analytical solutions is presented via calculating and simulating the hydraulics of groundwater flow in river–aquifer systems by means of integral transforms.


Introduction
Numerical and analytical models play an important role in scientific research in engineering. Therefore, researchers have been doing research on new techniques to solve differential equations. In the recent years, many computational methods have been proposed and developed for analyses of engineering problems (Avazzadeh et al. 2020;Nikan and Avazzadeh 2021;Nikan et al. 2022;Rasoulizadeh et al. 2021).
Analytical models have proposed in fluid mechanics subjects for analyses of engineering problems. For example, an accurate prediction of the behavior of the groundwater flow in aquifers in response to changes in surface water is of considerable importance in obtaining solutions for groundwater flow problems. Therefore, prediction of groundwater head in porous media is among the most important topics in the study of groundwater-surface water systems. Groundwater and surface water should be treated as an integrated system. Numerical and analytical models play an important role in assessing the future behavior of water table fluctuations in the groundwater-surface water systems. Although the numerical methods can easily deal with complex geometries as well as heterogeneity and anisotropy, the analytical methods are preferred because they consume rather less time to compute the problem than numerical methods. Moreover, analytical solutions are useful tools for analyzing the sensitivity of aquifer parameters.
Many analytical models have been developed to estimate groundwater head variations in response to recharge or pumping from wells (Zlotnik and Tartakovsky 2008;Lu et al. 2015), surface recharge (Rai and Manglik 1999;Teloglou et al. 2008), tidal fluctuations (Sun 1997;Tang and Jiao 2001;Huang et al. 2015) and stage changes in adjacent water bodies (Singh 2004a, b;Jiang and Tang 2015). The aforementioned studies assumed that the aquifer was isotropic. In fact, many soils exhibit a certain degree of anisotropy due to stratification associated with soil forming processes such as sedimentation, illuviation, compaction and particle orientation (Assouline and Or 2006). There are also several analytical attempts that take into account the anisotropy of the aquifer (Park and Zhan 2002;Intaraprasong and Zhan 2009;Fen and Yeh 2012;Wang et al. 2015).
The effects of anisotropy on the nature of groundwater variation have been investigated by several researchers. For example, Chang and Yeh (2007) presented an analytical solution to describe the hydraulic head distribution and flow system in an anisotropic unconfined aquifer with a sloping bed and arbitrarily located multiwells under transient recharge. They demonstrated that the water table is steeper in y-direction than in x-direction as the ratio of hydraulic conductivity in x-direction to y-direction increases. Manglik et al. (2013) presented an analytical solution of groundwater flow equation for unconfined, anisotropic, two-dimensional rectangular aquifer to predict water table variations in the aquifer in response to general time-varying intermittent recharge from multiple basins. They assumed that the aquifer is bounded with four constant head boundaries. The solution is obtained by using extended finite Fourier sine transform. The results showed a significant effect of anisotropy in hydraulic conductivity on the pattern and magnitude of the water table variations. It was found that the growth of the water table for isotropic aquifer always maintains higher elevation than the level of water table for anisotropic aquifer. Singh (2010) presented some generalized analytical solutions for groundwater head in a horizontal aquifer in the presence of subsurface drains. The aquifer is homogeneous and anisotropic and interacts with four surrounding streams of constant head. It was found that the isotropic case is overall characterized with higher water levels as compared to the anisotropic one. Manglik and Rai (2015) developed an analytical model to predict water table variations in an anisotropic aquifer in response to intermittently applied time-varying recharge from multiple heterogeneous basins and pumping from multiple wells. They considered no-flow boundary conditions and solved the equations by using the finite Fourier transform. In the process of investigation of the effects of anisotropy on water table, it was found that the growth of water table below the recharge basin and its surrounding region is more for the isotropic aquifer than that for the anisotropic aquifer. However, in the region away from the boundaries of the recharge basin, growth of the water table for anisotropic aquifer is more than that for the isotropic aquifer. Huang et al. (2011) presented an analytical solution for describing the head distribution in an anisotropic unconfined aquifer with a single pumping horizontal well parallel to a fully penetrating stream. They assumed the anisotropic aquifer has smaller vertical hydraulic conductivity than horizontal one. Given that the stream has a constant stage during the pumping period, they indicated that the anisotropic aquifer has larger stream depletion rate than the isotropic one.
In the present work, an effort is being made to investigate the effects of anisotropy on hydraulic head and flow rate in a rectangular confined aquifer adjoining to time-varying streams. In this paper, an anisotropic aquifer generally refers to an aquifer in which the hydraulic conductivity of the aquifer along x-direction is more than that along y-direction. Hence, the novelty of this paper is evaluating the anisotropy effects on groundwater hydraulic head as well as flow rate in a confined aquifer in contact with varying level boundaries. Therefore, a set of new analytical expressions are obtained by means of the Laplace and Fourier transforms and the solutions applicability is shown by the help of hypothetical examples. The results show that the presented new analytical solutions are accurate, robust and efficient. Therefore, the results indicate that the presented new analytical solutions are very effective in the simulation of the groundwater flow in river-aquifer systems. Furthermore, one of the advantages of the new analytical solutions is to investigate the sensitivity analysis of aquifer parameters, which has been carried out in this paper. Also, some other new analytical solutions for steady-state conditions and sudden fall in streams level are provided as well. The main applicability of the new analytical solutions is to investigate interactions between stream and aquifer. These new analytical solutions can also be used to evaluate aquifer response to gradual and sudden drop in stream stage. Also, the derived new analytical solutions could be used inversely to find the aquifer parameters. It can be mentioned that the presented new analytical solutions could further be used in many practical problems in stream-aquifer systems. Furthermore, it could be utilized for the validation of experimental and numerical models. Also, the results of the present new analytical solutions obtained will enable a better understanding regarding the modeling of the interaction between the river and the aquifer. Therefore, this research is a contribution to a better understanding of the fluxes between the river and the aquifer. Finally, the current study contributes to overcome common weaknesses of model applications, fulfills gaps in the existing literature and highlights the importance of the modeling process in planning sustainable management of groundwater resources.
The remainder of this paper is structured as follows: At first, new analytical solutions based on integral transforms to calculate and simulate the hydraulics of groundwater flow in river-aquifer systems are presented. Then, results and discussion are provided. Finally, in the last section conclusions are drawn.

Methodology
The geometry of the stream-aquifer system considered for study is shown in Fig. 1. A finite confined, anisotropic, incompressible and homogeneous aquifer is assumed to be surrounded with four streams of varying levels. The partial differential equation governing hydraulic head in a two-dimensional confined aquifer is taken to be: The initial and boundary conditions of the problem are: where h is the hydraulic head,S s is the specific storage and K x and K y are the hydraulic conductivity of the aquifer along x-axis and y-axis, respectively.x L and y L are the length and width of the aquifer, respectively. 1 , 2 , 3 and 4 are the positive constants signifying change rate of the streams at the north, east, south and west of the aquifer, respectively (thereafter called the rates of streams), h L + h 0 is the initial head of the system and h L is the final level of the streams. Defining (2)-(6) may be rewritten as follows: Dimensionless variables can be introduced as follows: where is the anisotropy ratio. Using Eq. (12), the governing equation with boundary conditions can be rewritten as follows: (2) h(x, y, 0) = h L + h 0 , The problem is decomposed into four parts, each of which has a three constant level boundaries and one varying level boundary, as shown in Fig. 2. Accordingly, the function H (X, Y, ) can be divided into four parts: Details of derivations of the solution for H 1 (X, Y, ) are as follows:

Part 1
The governing equation with boundary conditions for H 1 (X, Y, ) is: The Laplace transform can be defined as: where Λ denotes the Laplace transform of H and p is the Laplace variable. Taking the Laplace transform of Eq. (20) results in: and the associated Laplace-transformed initial and boundary conditions are as follows: The finite Fourier sine transform with respect to X can be defined as: where is the finite Fourier sine transform of Λ and n is the transform parameter. Applying the finite Fourier sine transform with respect to X in Eqs. (27)-(32) results in: Taking inverse finite Fourier sine transform of Eq. (42) results in: The following equation can be used to invert the Laplacedomain solutions into time-domain solutions: Dimensionless flow rate through the two-dimensional aquifer with a unit cross-sectional area can be stated as: Hence, the dimensionless flow rate along x-axis and y-axis for part 1 can be stated as follows: Similarly, the procedures to get the solutions of transient hydraulic head and flow rate for parts 2-4 can be summarized as follows:

Part 2
The governing equation with boundary conditions for H 2 (X, Y, ) is: In this part, the Laplace transform is applied with respect to and the finite Fourier sine transform is applied with respect to Y. Finally, the solution for H 2 (X, Y, ) is obtained as follows: And the components of dimensionless flow rate are as follows: Equations (57), (58) and (59) provide the transient hydraulic head, flow rate along x-axis and flow rate along y-axis, respectively, for a confined aquifer with a varying level boundary at the east end and three constant level boundaries at the other ends.

Part 3
The governing equation with boundary conditions for H 3 (X, Y, ) is: In this part, the Laplace transform is applied with respect to and the finite Fourier sine transform is applied with respect X. Finally, the solution for H 3 (X, Y, ) is obtained as follows: . and the components of dimensionless flow rate are as follows: Equations (66), (67) and (68) provide the transient hydraulic head, flow rate along x-axis and flow rate along y-axis, respectively, for a confined aquifer with a varying level boundary at the south end and three constant level boundaries at the other ends.

Part 4
The governing equation with boundary conditions for H 4 (X, Y, ) is: In this part, the Laplace transform is applied with respect to and the finite Fourier sine transform is applied with respect Y. Finally, the solution for H 4 (X, Y, ) is obtained as follows: and the components of dimensionless flow rate are as follows: Equations (75), (76) and (77) provide the transient hydraulic head, flow rate along x-axis and flow rate along y-axis, respectively, for a confined aquifer with a varying level boundary at the west end and three constant level boundaries at the other ends.

Comparison of the new analytical solutions with numerical model of MODFLOW
In order to demonstrate the capability of the proposed new analytical solutions in the estimation of groundwater hydraulic head in an anisotropic aquifer, another hypothetical aquifer is simulated in numerical model of MOD-FLOW. The rates of streams are: 1 = 0.05 day −1 , 2 = 0 , 3 = 0.1 day −1 and 4 = 0.2 day −1 . The anisotropy ratio is = 3 ( K x = 30 m∕day , K y = 10 m∕day ). The length and width of the aquifer are selected as: x L = 100 m and y L = 50 m . The other parameters are given in Table 1. For the numerical model of MODFLOW, the domain is discretized into 500 × 250 rectangular cells with Δx = Δy = 0.2 m . The computational time interval for numerical model of MODFLOW is 12 h. The type of boundary conditions is considered as specified hydraulic head. The solver and flow packages are PCG2 and LPF, respectively. LPF has the ability to enter horizontal anisotropy values on a cell-bycell basis. The values of hydraulic head at t = 10 days and at y = 25 m are shown in Fig. 3. Figure 3 shows that the results of the new analytical solutions are in good agreement with those results obtained from numerical model of MODFLOW.

Results and discussion
The configuration of the personal computer used to perform the simulation results is that the CPU is Intel(R)/ Core(TM)4/i7-8550U, CPU@4.000 GHz. As mentioned earlier, Eq. (19) presents the new analytical solution of an aquifer bounded with four time-varying level boundaries, which is sum of Eqs. (45), (57), (66) and (75). Notice that setting i = 0 (i = 1, 2, 3, 4) gives constant level boundaries. Therefore, the new analytical expressions for other types of boundary configurations with both constant and time-varying levels can readily be obtained. Table 2 presents the new analytical expressions for 15 different boundary configurations. Thereafter, BC stands for boundary configuration and the following number indicates the type of boundary configuration given in Table 2. For example, BC5 refers to the boundary configuration number 5 in Table 2.
Steady state of Eq. (45) is: And the associated flow rates can be expressed as follows: Steady state of Eq. (57) is: And the associated flow rates can be expressed as follows: Steady state of Eq. (66) is: And the associated flow rates can be expressed as follows: Steady state of Eq. (75) is: And the associated flow rates can be expressed as follows:  (84) and (87) are very useful new exact expressions describing the groundwater hydraulic head distribution in an anisotropic 2D aquifer in steadystate conditions. In a one-dimensional aquifer connected to two streams with different stages, the initial groundwater head can be assumed to vary linearly between the streams. However, estimation of the initial groundwater hydraulic head in a two-dimensional aquifer connected to four streams with different stages is more complicated. Equations (78), (81), (84) and (87) can be used for such these conditions.

Sudden fall in streams level
For the case of sudden fall in streams level i (i = 1, 2, 3, 4) goes to infinity in Eqs. (45), (57), (66) and (75). Sudden fall in level of the north stream is: Sudden fall in level of the east stream is: Sudden fall in level of the south stream is: Sudden fall in level of the west stream is: Groundwater hydraulic head distribution for three scenarios of sudden fall in streams level at t = 1 s is illustrated in Fig. 4. The figures are plotted for three different BCs. Scenario 1: Fig. 4a shows groundwater hydraulic head distribution for sudden changes in streams level of BC1. In this scenario, the water level falls suddenly at the north stream and is constant at the other streams. Thus, the results can be obtained by using Eq. (90). Figure 4b shows groundwater hydraulic head distribution for sudden changes in streams level of BC5. Scenario 2: In this scenario, the water level falls suddenly at the north and east streams and is constant at the other streams. Here, the results can be obtained by sum of Eqs. (90) and (91). Scenario 3: Fig. 4c shows groundwater hydraulic head distribution for sudden changes in streams level of BC14. In this scenario, the water level falls suddenly at the east, south and west streams and is constant at the north stream. Here, the results can be obtained by sum of Eqs. (91), (92) and (93). The anisotropy ratio is = 3 ( K x = 30 m∕day ,K y = 10 m∕day ), and the other parameters are as before. Figure 5a, b and c shows the groundwater hydraulic head variations with time at different points of the aquifer for scenarios 1, 2 and 3, respectively. Five different points are considered, namely P1:x = 20 , y = 30 m , P2:x = 80 , y = 20 m , P3:x = 50 , y = 50 m , P4:x = 30 , y = 80 m and P5:x = 70 , y = 60 m . It can be seen that the groundwater hydraulic head stabilizes about less than 0.01 s.

Effects of anisotropy on aquifer response
In this section, several examples are conducted to investigate the effects of anisotropy on hydraulic head and flow rate. Here, four cases are defined. The rates of streams for these cases are given in Table 3, the anisotropy ratio is = 3 and the other parameters are kept as before. Figure 6 shows the values of hydraulic head over distance at 10th day for both isotropic and anisotropic aquifer and for cases 1-4. It can be observed from Fig. 6 that with increases in K y the hydraulic head increases in cases 1 and 3 but decreases in cases 2 and 4. This phenomenon is valid for any arbitrary time, which is also found in Fig. 7. Figure 7 shows the hydraulic head variations with time at a point located at x = 50 m, y = 25 m. This phenomenon is because that a higher hydraulic conductivity along y-direction deviates the groundwater flow direction along this direction. On the other hand, in cases 2 and 4 the streams level at the south and north boundaries are falling rather faster than the other streams. This causes the groundwater direction to also deviate along y-axis. Therefore, an isotropic  = 1 s a scenario 1, b scenario 2 and c scenario 3 aquifer K y = 30 (m∕day) allows the water to flow faster in y-direction, and as a result, the groundwater hydraulic head decreases. This is exactly vice versa in cases 1 and 3 in which the groundwater hydraulic head is deviated along x-direction due to the available gradient as well as a smaller K y . Therefore, in the presence of anisotropy the groundwater hydraulic head decreases in cases 1 and 3. Figures 8, 9, 10 and 11 show the sensitivity of hydraulic head to changes in rates of the streams in the isotropic ( = 1 ∶ K y = K x = 30 m∕day ) and anisotropic ( = 3 ∶ K y = 10 m∕day, K x = 30 m∕day ) aquifers for t = 10 days and at a point located at x = 70 m, y = 40 m . These figures are plotted for different rates 1 , 2 , 3 and 4 , keeping the other parameters as before. Figure 7 shows that the magnitude of difference between corresponding hydraulic head for 1 = 0.01 and 1 = 1 day −1 is about 0.52 m for the isotropic aquifer and it is about 0.14 m for the anisotropic aquifer. Also, Fig. 8 shows that the magnitude of difference between corresponding hydraulic head for 2 = 0.01 and 2 = 1 day −1 is about 1.64 m for the isotropic aquifer and it is about 2.24 m for the anisotropic aquifer. Moreover, Fig. 9 shows that the magnitude of difference between corresponding hydraulic head for 3 = 0.01 and 3 = 1 day −1 is about 1.05 m for the isotropic aquifer and it is about 0.42 m for the anisotropic aquifer. In addition, Fig. 10 shows that the magnitude of difference between corresponding hydraulic head for 4 = 0.01 and 4 = 1 day −1 is about 0.41 m for the isotropic aquifer and it is about 0.82 m for the anisotropic aquifer.

Effects of anisotropy on hydraulic head
In summary, it can be concluded from Figs. 7, 8, 9 and 10 that the hydraulic head in an anisotropic aquifer is more sensitive to changes in 2 and 4 and less sensitive to changes in 1 and 3 than that in an isotropic aquifer. This is because a large value of hydraulic conductivity along y-direction raises the groundwater flow velocity along this direction. Hence, with increases in K y the groundwater flow tends to move toward either the north or south streams and consequently the hydraulic head gets more sensitive to variations of 1 and 3 in isotropic aquifer. Table 4 gives the values of flow rate and the direction of the flow ( ) at different points of the aquifer. is the angle  between vector q and positive axis. To investigate the effects of anisotropy on flow rate and its direction, two values of hydraulic conductivity are considered, namely K y = 30 and K y = 10 (m∕day) . Here, 1 = 2 = 3 = 4 = 0.02 (day −1 ) and the other parameters are given in Table 1. As expected, Table 4 shows that as K y increases, the groundwater flow tends to deviate along y-axis and vice versa. This is because the higher hydraulic conductivity allows the water to flow easier in aquifer. Tables 5, 6, 7 and 8 provide the dimensionless flow rate at the left, right, north and south interfaces for different values of K y , respectively. Table 5 shows that at the left interface, the average difference between corresponding dimensionless flow rate ( Q ) of K y = 30 and K y = 10 (m∕day) is 0.73, 0.1, 0.13 and 0.1 for cases 1 to 4, respectively. Table 6 shows that at the right interface, the average difference between corresponding dimensionless flow rate ( Q ) of K y = 30 and K y = 10 (m∕day) is 0.13, 0.1, 0.73 and 0.1 for cases 1-4, respectively. Table 7 shows that at the north interface, the average difference between corresponding dimensionless flow rate ( Q ) of K y = 30 and K y = 10 (m∕day) is 0.42, 0.12, 0.42 and 0.04 for cases 1 to 4, respectively. Table 8 shows that at the south interface, the average difference between corresponding dimensionless flow rate ( Q ) of K y = 30 and K y = 10 (m∕day) is 0.42, 0.04, 0.42 and 0.12 for cases 1-4, respectively.

Effects of anisotropy on flow rate
Tables 5, 6, 7 and 8 show that with rises in K y , values of flow rate at different points of interfaces may increase or decrease. However, it can be seen that the average dimensionless flow rate at interfaces increases with increasing K y . The average of the dimensionless flow rate of interfaces in cases 1 to 4 is increased about 0.425 , 0.09 , 0.425 and 0.09 , respectively, in the isotropic aquifer. It can be highlighted that in case 1 (when the rate of the left stream is more than that of the other streams), the dimensionless flow rate by a change in K y is more influenced at the left boundary than that at the other boundaries. Similarly, in case 3 (when the rate of the right stream is more than that of the other streams) the dimensionless flow rate by a change in K y is more influenced at the right boundary than that at the other boundaries. However, this phenomenon is vice versa in cases 2 and 4. In case 2 (when the rate of the south stream is more than that of the other streams), the dimensionless flow rate by a change in K y is less influenced at the south boundary than that at the other boundaries. Moreover, in case 4 (when the rate of the north stream is more than that of the other streams) the dimensionless flow rate by a change in K y is less influenced at the north boundary than that at the other boundaries.

Conclusions
A set of dimensionless new analytical solutions are proposed to describe transient hydraulic head and flow rate in a finite confined, anisotropic, incompressible and homogeneous aquifer with four time-varying streams located at its ends. Using superposition principle, the problem is decomposed into four parts and the final solution is obtained by the sum of the solutions. In part one, a new analytical solution describing the hydraulic head for an aquifer with a time-varying stream at the north boundary and three constant streams at the other boundaries is derived. The associated new analytical expression for this case is Eq. (45). Part 2 deals with an aquifer with a time-varying stream at the east boundary and three constant streams at the other   are very effective in the simulation of the groundwater flow in river-aquifer systems. Furthermore, one of the advantages of the new analytical solutions is to investigate the sensitivity analysis of aquifer parameters, which has been carried out in this paper. Also, some other new analytical solutions for steady-state conditions and sudden fall in streams level are provided as well. Some significant new analytical expressions describing the sudden fall in streams level as well as steady-state conditions are derived. The associated new analytical expressions of steady-state conditions for parts 1 to 4 are given in Eqs. (78)-(89). Equations (78), (81), (84) and (87) are very important new exact expressions and can be used to calculate the initial groundwater hydraulic head in an anisotropic aquifer when the initial levels of the surrounding streams are different. Moreover, the new analytical expressions of sudden fall in streams level for parts 1 to 4 are given in Eqs.
(90)-(93). Using these equations, the groundwater hydraulic head distribution for three scenarios of sudden fall in streams level is presented as well.
Several examples are provided to investigate the effects of anisotropy ( K y < K x ) on hydraulic head and flow rate. The following pertinent conclusions are drawn: • It is shown when the rates of the north and south streams are more than those of the other streams, the hydraulic head for isotropic aquifer is less than that for anisotropic aquifer and when the rates of the west and east streams are more than that of the other streams, the hydraulic head for isotropic aquifer is more than that for anisotropic aquifer. • The hydraulic head in an anisotropic aquifer is more sensitive to changes in rates of the east and west streams and less sensitive to changes in rates of the north and south streams than that in an isotropic aquifer. • Furthermore, it is demonstrated that an isotropic aquifer provides more water exchange at the interfaces than an anisotropic aquifer. An anisotropy ratio of = 3 reduced the average of the dimensionless values of flow rate at boundaries about 0.425, 0.09, 0.425 and 0.09 for cases 1 to 4, respectively. In other words, neglecting anisotropy of the aquifer leads to overestimate the flow rate at interfaces.
In this research, we showed the significance of consideration of the rates of the streams when dealing with anisotropic aquifers. Actually, the hydraulic conductivity is a controlling parameter in determination of flow direction in the aquifer. This parameter shows that the aquifer is how and how much sensitive to each stream. Finally, it can be stated that neglecting the effects of anisotropy causes to wrongly predict the flow path line, and as a result, the researchers will be misled in estimation of travel time and contaminant distribution in the aquifer.
The main applicability of the new analytical solutions is to investigate interactions between stream and aquifer. These new analytical solutions can also be used to evaluate aquifer response to gradual and sudden drop in stream stage. Also, the derived new analytical solutions could be used inversely to find the aquifer parameters. It can be mentioned that the presented new analytical solutions could further be used in many practical problems in stream-aquifer systems. Furthermore, it could be utilized for the validation of experimental and numerical models. Also, the results of the present new analytical solutions obtained will enable a better understanding regarding the modeling of the interaction between the river and the aquifer. Therefore, this research is a contribution to a better understanding of the fluxes between the river and the aquifer. Finally, the current study contributes to overcome common weaknesses of model applications, fulfills gaps in the existing literature and highlights the importance of the modeling process in planning sustainable management of groundwater resources.
The novelty of this paper is evaluating the anisotropy effects on groundwater hydraulic head as well as flow rate in a confined aquifer in contact with varying level boundaries. Therefore, a set of new analytical expressions are obtained by means of the Laplace and Fourier transforms and the solutions applicability is shown by the help of hypothetical examples.
Funding The authors received no specific funding for this work.