Analytical solutions for tidal propagation in a rectangular basin

The tidal propagation in a rectangular basin of an estuary has been studied in this paper. The dynamic response of the tidal current and elevation inside the basin has been obtained using the concept of Kelvin waves and spectrum of Poincare waves. In this study, we have found that initial tidal elevations are amplified near the closed end of the basin or at the tidal barrier site.     

Digital computer solutions of groundwater flow

This paper examines critically the use of the alternating direction implicit method for solving groundwater flow problems. Comparisons with alternative solutions show that inaccuracies can arise if the time increments are too large. A method is suggested whereby the computer chooses an optimum time increment. The representation of recharge, infiltration and wells is considered and a detailed example compares the alternating direction implicit, the explicit and the resistance-capacitance network methods.     

Analytical solutions of the advection‐dispersion equation for transient groundwater flow. A numerical validation

An analytical transport‐model was developed to simulate the propagation of a contaminant in one‐ and two‐dimensional transient flow in groundwater. It is proved that the distribution of concentration at a given time and for a given discharge is identical to that obtained for a different discharge if the volumetric flux of water is the same in the two cases. The results of simulations have been compared with results obtained using the MT3DMS numerical model. There is good agreement when the calculated concentrations are flux‐weighted concentrations. On the other hand, there is a notable divergence when the resident mode is considered. Resident mode concentrations express the mass per unit volume whereas flux mode concentrations express the ratio of mass flux to fluid flux. The solutions presented in this paper can thus be a useful alternative to code MT3DMS when the objective is to simulate concentrations in transient flow according to a resident mode. Copyright © 2008 John Wiley & Sons, Ltd.     

Analytical stochastic solutions of saltwater/freshwater interface in coastal aquifers

Three analyitcal solutions of saltwater intrusion under uncertain hydrologic, hydrogeologic, and pumping conditions are presented. These solutions extends the existing deterministic, sharp interface solutions to stochastic ones. The randomness is represented in terms of statistical measures of mean, standard deviation and covariance. The analysis is based on perturbation using Taylor series expansion. Simulations based on probability distribution are conducted for verification.     

Analytical stochastic solutions of saltwater/freshwater interface in coastal aquifers

Three analyitcal solutions of saltwater intrusion under uncertain hydrologic, hydrogeologic, and pumping conditions are presented. These solutions extends the existing deterministic, sharp interface solutions to stochastic ones. The randomness is represented in terms of statistical measures of mean, standard deviation and covariance. The analysis is based on perturbation using Taylor series expansion. Simulations based on probability distribution are conducted for verification.     

Closed-form analytical solutions incorporating pumping and tidal effects in various coastal aquifer systems

Pumping wells are common in coastal aquifers affected by tides. Here we present analytical solutions of groundwater table or head variations during a constant rate pumping from a single, fully-penetrating well in coastal aquifer systems comprising an unconfined aquifer, a confined aquifer and semi-permeable layer between them. The unconfined aquifer terminates at the coastline (or river bank) and the other two layers extend under tidal water (sea or tidal river) for a certain distance L. Analytical solutions are derived for 11 reasonable combinations of different situations of the L-value (zero, finite, and infinite), of the middle layer’s permeability (semi-permeable and impermeable), of the boundary condition at the aquifer’s submarine terminal (Dirichlet describing direct connection with seawater and no-flow describing the existence of an impermeable capping), and of the tidal water body (sea and tidal river). Solutions are discussed with application examples in fitting field observations and parameter estimations.     

A two‐dimensional analytical solution for groundwater flow in a leaky confined aquifer system near open tidal water

Groundwater in coastal areas is commonly disturbed by tidal fluctuations. A two‐dimensional analytical solution is derived to describe the groundwater fluctuation in a leaky confined aquifer system near open tidal water under the assumption that the groundwater head in the confined aquifer fluctuates in response to sea tide whereas that of the overlying unconfined aquifer remains constant. The analytical solution presented here is an extension of the solution by Sun for two‐dimensional groundwater flow in a confined aquifer and the solution by Jiao and Tang for one‐dimensional groundwater flow in a leaky confined aquifer. The analytical solution is compared with a two‐dimensional finite difference solution. On the basis of the analytical solution, the groundwater head distribution in a leaky confined aquifer in response to tidal boundaries is examined and the influence of leakage on groundwater fluctuation is discussed. Copyright © 2001 John Wiley & Sons, Ltd.     

Semi-analytical solutions of groundwater flow in multi-zone (patchy) wedge-shaped aquifers

Alluvial fans are potential sites of potable groundwater in many parts of the world. Characteristics of alluvial fans sediments are changed radially from high energy coarse-grained deposition near the apex to low energy fine-grained deposition downstream so that patchy wedge-shaped aquifers with radial heterogeneity are formed. The hydraulic parameters of the aquifers (e.g. hydraulic conductivity and specific storage) change in the same fashion. Analytical or semi-analytical solutions of the flow in wedge-shaped aquifers are available for homogeneous cases. In this paper we derive semi-analytical solutions of groundwater flow to a well in multi-zone wedge-shaped aquifers. Solutions are provided for three wedge boundary configurations namely: constant head–constant head wedge, constant head–barrier wedge and barrier–barrier wedge. Derivation involves the use of integral transforms methods. The effect of heterogeneity ratios of zones on the response of the aquifer is examined. The results are presented in form of drawdown and drawdown derivative type curves. Heterogeneity has a significant effect on over all response of the pumped aquifer. Solutions help understanding the behavior of heterogeneous multi-zone aquifers for sustainable development of the groundwater resources in alluvial fans.     

Influence of variable salinity conditions in a tidal creek on riparian groundwater flow and salinity dynamics

While recent studies have revealed that tidal fluctuations in an estuary significantly affect groundwater flows and salt transport in the riparian zone, only seawater salinity in the estuary has been considered. A numerical study is conducted to investigate the influence of estuarine salinity variations on the groundwater flow and salt dynamics in the adjacent aquifer to extend our understanding of these complex and dynamic systems. Tidal salinity fluctuations (synchronous with estuary stage) were found to alter the magnitude and distribution of groundwater discharge to the estuary, which subsequently impacted on groundwater salinity patterns and residence times, especially in the riparian zone. The effects of salinity fluctuations were not fully captured by adopting a constant, time-averaged estuarine salinity. The modelling analysis also included an assessment of the impact of a seasonal freshwater flush in the estuary, similar to that expected in tropical climates (e.g. mean estuary level during flood significantly greater than average), on adjacent groundwater flow and salinity conditions. The three-month freshwater flushing event temporarily disrupted the salt distribution and re-circulation patterns predicted to occur under conditions of constant salinity and tidal water level fluctuations in the estuary. The results indicate that the salinity variations in tidal estuaries impact significantly on estuary–aquifer interaction and need to be accounted for to properly assess salinity and flow dynamics and groundwater residence times of riparian zones.     

Analytical solutions for water flow passing over a vegetal area

This study is aimed at investigating the vertical velocity profile of flow passing over a vegetal area by an analytical approach. The soil ground is considered as pervious and thus non-zero velocity at the ground surface can be estimated. The soil and vegetation layers are regarded as homogeneous and isotropic porous media. Therefore the solution of the flow can be obtained by applying the theory of turbulent flow and Biot’s theory of poroelasticity after dividing the flow field into three layers: homogenous water, vegetation and pervious soil. The velocity distribution is compared with the experimental data of [Rowiński PM, Kubrak J. A mixing-length model for predicting vertical velocity distribution on flows through emergent vegetation. J Hydrol Sci 2002;47(6):893–904] to show its validity. In addition, five dimensionless parameters denoting the variation of slope, permeability of soil, Reynolds stress, density of vegetation, and relative height of vegetation are proposed to reveal their effects on the surface water flow. The analytical solutions of flow velocity can also be simplified into simpler expressions to describe the flow passing over a non-vegetated area.     

Analytical solutions for flow to a well through a fault

We present explicit analytical solutions to problems of steady groundwater flow to a pumping well in an aquifer divided by an infinite, linear fault. The transmissivity of the aquifer is allowed to jump from one side of the fault to the other to model the juxtaposition of host rocks with different hydrologic properties caused by faulting. The fault itself is represented as a thin anisotropic inhomogeneity; this allows the fault to act as a combined conduit–barrier to groundwater flow, as is commonly described in the literature. We show that the properties of the fault may be represented exactly by two lumped parameters—fault resistance and fault conductance—and that the effects of the fault on flow in the adjacent aquifer is independent of the fault width. We consider the limiting cases of a purely leaky and a purely conductive fault where the fault domain may be replaced exactly by internal boundary conditions, and we investigate the effects of fault properties on the flow behavior in the adjacent aquifers. We demonstrate that inferring fault properties based on field observations of head in the aquifer is inherently difficult, even when the fault may be described by one of the two limiting cases. In particular, the effects of a leaky fault and a conductive fault on heads and discharges in the aquifer opposite the fault from the well, are shown to be identical in some cases.     

Semigroup solutions to stochastic unsteady groundwater flow subject to random parameters

Two methods for the solution of partial differential equations (PDE) for the general case of random in time physical parameters are presented and their application to the solution of unsteady regional groundwater flow equations are illustrated. The first method is the semigroup approach which directly offers a solution without resorting to closure approximations (hierarchy techniques), perturbation techniques, or Montecarlo simulation techniques. The semigroup approach can also handle the general stochastic problem when randomness also appears as initial conditions, boundary conditions or forcing terms. The second method is an approximation scheme to obtain the semigroup solution in complex cases and permits the solution of equations with more than one random coefficient.     

Analytical solutions for 2D topography-driven flow in stratified and syncline aquifers

Two analytical solution methods are presented for regional steady-state groundwater flow in a two-dimensional stratified aquifer cross section where the water table is approximated by the topographic surface. For the first solution, the surficial aquifer is represented as a set of dipping parallel layers with different, but piecewise constant, anisotropic hydraulic conductivities, where the anisotropy is aligned with the dip of the layered formation. The model may be viewed as a generalization of the solutions developed by [Tóth JA. A theoretical analysis of groundwater flows in small drainage basins. J Geophys Res 1963;68(16):4795–812; Freeze R, Witherspoon P. Theoretical analysis of regional groundwater flow 1) analytical and numerical solution to the mathematical model, water resources research. Water Resour Res 1966;2(4):641–56; Selim HM. Water flow through multilayered stratified hillside. Water Resour Res 1975;11:949–57] to an multi-layer aquifer with general anisotropy, layer orientation, and a topographic surface that may intersect multiple layers. The second solution presumes curved (syncline) layer stratification with layer-dependent anisotropy aligned with the polar coordinate system. Both solutions are exact everywhere in the domain except at the topographic surface, where a Dirichlet condition is met in a least-squared sense at a set of control points; the governing equation and no-flow/continuity conditions are met exactly. The solutions are derived and demonstrated on multiple test cases. The error incurred at the location where the layer boundaries intersect the surface is assessed.     

Effect of variable‐density groundwater flow on nitrate flux to coastal waters

A two‐dimensional variable‐density groundwater flow and transport model was developed to provide a conceptual understanding of past and future conditions of nitrate (NO₃) transport and estimate groundwater nitrate flux to the Gulf of Mexico. Simulation results show that contaminant discharge to the coast decreases as the extent of saltwater intrusion increases. Other natural and/or artificial surface waters such as navigation channels may serve as major sinks for contaminant loading and act to alter expected transport pathways discharging contaminants to other areas. Concentrations of NO₃ in the saturated zone were estimated to range between 30 and 160 mg?L^?1 as NO₃. Relatively high hydraulic vertical gradients and mixing likely play a significant role in the transport processes, enhancing dilution and contaminant migration to depth. Residence times of NO₃ in the deeper aquifers vary from 100 (locally) to about 300 years through the investigated aquifer system. NO₃ mass fluxes from the shallow aquifers (0 to 5.7 × 10⁴ mg?m^?2?day^?1) were primarily directed towards the navigation channel, which intersects and captures a portion of the shallow groundwater flow/discharge. Direct NO₃ discharge to the sea (i.e. Gulf of Mexico) from the shallow aquifer was very low (0 to 9.0 × 10¹ mg · m^?2?day^?1) compared with discharge from the deeper aquifer system (0 to 8.2 × 10³ mg?m^?2?day^?1). Both model‐calibrated and radiocarbon tracer‐determined contaminant flux estimates reveal similar discharge trends, validating the use of the model for density‐dependent flow conditions. The modelling approach shows promise to evaluate contaminant and nutrient loading for similar coastal regions worldwide. Copyright © 2015 John Wiley & Sons, Ltd.     

Use of time series analysis for the identification of tidal effect on groundwater in the coastal area of Kimje, Korea

The purpose of this study is to apply time series analysis to investigate whether the groundwater quality in the coastal area is affected by the tide. Continuous and regular in situ monitoring data of electrical conductivity (EC) and groundwater level, and tidal level data measured by the National Oceanographic Research Institute were used for the time series analysis. Through the time series analysis, it is known that EC and groundwater level conspicuously fluctuate with two periodicities (15.4 and 0.52-day), which is very similar to those of the tide. Also the behaviors of their fluctuations vary in accordance with the tidal period. These indicate that the groundwater quality has been mainly controlled by the tidal level, and the strength of tidal effect on the groundwater quality is different according to the tidal period.     

Analytic solutions for interface flow in combined confined and semi-confined,coastal aquifers

Analytic solutions are presented for steady interface flow in aquifers consisting of a confined and a semi-confined section. The total discharge is constant in the confined section and is directed towards the semi-confined section, which is bounded on top by a leaky layer that separates the aquifer from the sea. All solutions are based on the Dupuit approximation. The interface position is fully determined by two dimensionless parameters. The first parameter is the product of the uniform gradient towards the sea and the dimensionless leakage factor divided by the dimensionless density difference. The second dimensionless parameter is the length of the semi-confined section divided by the leakage factor. Four types of flow are distinguished. For cases I and II the interface does not reach the end of the semi-confined section, while it does for cases III and IV. For cases I and III the interface extends into the confined section, while for cases II and IV it remains entirely in the semi-confined section. Analytic solutions are presented for the position of the interface for all four cases. Diagrams are presented to determine the type of flow and the position of the interface toe based on the two dimensionless parameters. The pre-development position of the interface along the Georgia coast is computed as an illustration of the practical application of the presented formulas. Another practical application is the estimation of how far a numerical model of seawater intrusion should be extended into the sea for accurate simulations.     

Analytical studies on transient groundwater flow induced by land reclamation using different fill materials

Land reclamation may have a significant influence on groundwater regimes. Analytical solutions have been developed in the past to study the impact of land reclamation on a steady‐state groundwater flow and transient flow in fill materials, assuming that the reclamation site consists of a single zone of uniform hydraulic parameters. In this paper, we derive analytical solutions to describe the transient water table change in response to multi‐stage land reclamation where the fill material is uniform in each stage but the hydraulic conductivity of the fill material varies from stage to stage. By introducing the method of separation of variables, we develop a transient analytical solution to study the impact of land reclamation consisting of fill material with different hydraulic properties on groundwater dynamics. The results show that the water table first increases significantly into the reclaimed zone following the fill material deposition, and then the increase gradually propagates into the original aquifer. The change of water table in the original aquifer mainly depends on the value of hydraulic conductivity of the fill materials. Examples in this paper illustrate how the aquifer system experiences a long time unsteady‐state flow as a result of the reclamation, and it takes at least tens of years for the system to approach a new equilibrium. It is suggested that for a large‐scale reclamation project, the response of the groundwater regime to reclamation should be carefully studied. Copyright © 2013 John Wiley & Sons, Ltd.     

New analytical solutions for groundwater flow in wedge-shaped aquifers with various topographic boundary conditions

The hydraulic head distribution in a wedge-shaped aquifer depends on the wedge angle and the topographic and hydrogeological boundary conditions. In addition, an equation in terms of the radial distance with trigonometric functions along the boundary may be suitable to describe the water level configuration for a valley flank with a gentle sloping and rolling topography. This paper develops a general mathematical model including the governing equation and a variety of boundary conditions for the groundwater flow within a wedge-shaped aquifer. Based on the model, a new closed-form solution for transient flow in the wedge-shaped aquifer is derived via the finite sine transform and Hankel transform. In addition, a numerical approach, including the roots search scheme, the Gaussian quadrature, and Shanks’ method, is proposed for efficiently evaluating the infinite series and the infinite integral presented in the solution. This solution may be used to describe the head distribution for wedges that image theory is inapplicable, and to explore the effects of the recharge from various topographic boundaries on the groundwater flow system within a wedge-shaped aquifer.     

Spectral responses of water level in tidal river and groundwater

The water level of five river stages and seven groundwater wells in the Taipei Basin were analysed by spectral analysis in the frequency domain. The diurnal, semi‐diurnal and quarter‐diurnal tidal components of the Tanshui River appear to relate closely to astronomical tides as K₁, M₂ and M₄, respectively. It is also found that the diurnal component reveals a reversed phase angle in the middle section of the Tanshui River; the phase of the quarter‐diurnal component is also found to be reversed at stations upstream in the Tanshui River and Hsintien Stream. It is believed that these phenomena could be caused by local variation in the river channel topography. The autospectrum and cross‐spectrum between groundwater elevation and nearby river stage were observed to correlate highly with the frequency of the astronomical tides K₁, M₂ and M₄. From the study of the phase shift and time lag of water level fluctuations at river stages and groundwater wells, it was found that the tidal effects of diurnal, semi‐diurnal, and quarter‐diurnal components were significantly different. The relationships between phase and the fluctuated range of atmospheric pressure and water level imply that change in atmospheric pressure does not affect water level fluctuation in the river stage and groundwater well. Copyright © 1999 John Wiley & Sons, Ltd.