A generalized analytical solution of groundwater head response to dual tide in a multilayered island leaky aquifer system

The hydraulic properties of coastal aquifer systems are relevant to various hydrogeological, hydro-ecological and engineering problems. This study presents an analytical solution for predicting groundwater head fluctuations induced by dual-tide in multilayered island aquifer systems, consisting of an unconfined aquifer on the top and any number of leaky aquifers below. The solution was derived via the methods of matrix differential calculus and separation of variables. It is more general than any existing analytical solutions for the tidal pressure propagation since the new solution can consider multilayered aquifer systems along with the effects of leakage and aquifer length. Using this solution, we illustrated potential errors that may occur due to neglecting one or more vital factors affecting groundwater fluctuations. Besides, we articulated the groundwater response to the dual-tide in complex coastal aquifers. Considering that some thin semipermeable layers may be ignored in practical field investigation, we also demonstrated the effects due to simplification of aquifer layers. The results showed that with the increase in the number of overlapped leaky layers, the tidal propagation in the bottom part of multilayered aquifer system approaches that in a single confined aquifer with the same transmissivity and storage.

Conceptualization and development of multi-layered groundwater model in transient condition

Effective management of water resource is essential in arid and semi-arid areas of India. In Bihar, for drinking purpose humans, livestock is dependent on the groundwater as well as in agricultural areas groundwater plays an important role in irrigation directly or indirectly. There is rise in the groundwater demand due to rapid population increase and fast industrialization. To meet this groundwater demand, excessive withdrawal of groundwater is a point of concern due to limited storage of it. Assessment of the groundwater was done by preparing a numerical model of the groundwater flow. This model is capable of solving large groundwater problems and associated complexity with it. In this study, a transient multi-layered groundwater flow model was conceptualized and developed for the Koshi River basin. In north Bihar plains, the Koshi River is one of the biggest tributaries of the Ganga River system. Koshi originates from the lower part of Tibet and joins the Ganga River in Katihar district, Bihar, India. After model development, calibration of the model was also done, by considering three model parameters, to represent the actual field conditions. For validation of the model, fifteen observation wells have been selected in the area. With the help of observation well data, computed and observed heads were compared. Comparison results have been found to be encouraging and the computed groundwater head matched with the observed water head to a realistic level of accuracy. Developed groundwater model is used to predict the groundwater head and flow budget in the concerned area. The study revealed that groundwater modeling is an important method for knowing the behavior of aquifer systems and to detect groundwater head under different varying hydrological stresses. This type of study will be beneficial for the hydrologist and water resource engineers to predict the groundwater flow behavior, before implementing any project or to implement a correction scheme.

Aquifer Response to Stream-Stage Fluctuations: Field Tests and Analytical Solution for a Case Study of the Yangtze River in Wuhan, China

In this paper, field hydrological monitoring tests and an analytical solution for the aquifer response (i.e., the groundwater head) to stream-stage fluctuations are presented through a case study of the Yangtze River in Wuhan, China. A sinusoidal equation for the stream water level over time is proposed and is verified by the monitoring results. Based on the classical 1-D solution of ground water flow equations, a new flow model was proposed to analyze the confined aquifer response to stream-stage fluctuations. The groundwater head of the confined aquifer can then be calculated by embedding the proposed sinusoidal equation of the stream-stage fluctuation into the new flow model. Field monitoring data of the three hydrological years from 2017 to 2019 were compiled to evaluate the proposed model. The results show that the predicted groundwater head of the confined aquifer matches well with the monitoring data. The prediction error of the proposed model (about 8%) is much lower than that of the previous model (about 26%). This paper presents an improved analytical model with a satisfying capability to predict the response of the groundwater head of the confined aquifer to stream-stage fluctuation. This model can be applied into groundwater flow models (such as MODFLOW, FEFLOW, etc.). The results of this study provide a valuable hydrological resource for underground practical engineering in the riverside area, especially in the first terrace of the Yangtze River, and for the dewatering design of a foundation pit in a riverside area.

Integrated Hydrological Analysis of Little Akaki Watershed Using SWAT-MODFLOW, Ethiopia

In Ethiopia, groundwater is the main source of freshwater to support human consumption and socio-economic development. Little Akaki watershed is located in Upper Awash basin, known for its high annual rainfall and considered as the potential groundwater recharge zone. On the contrary, urbanization and industrial expansion are increasing at an alarming rate in the area. This became a concern threatening the groundwater resources’ sustainability. To address these challenges, integrated analysis of groundwater recharge and groundwater numerical simulations were made. For groundwater recharge estimation, SWAT model was used. The result indicated that recharge in the watershed mostly occurs from July to October with maximum values in August. On average, the estimated annual catchment recharge was 179 mm. For the numerical simulation and prediction of the groundwater flow system, MODFLOW 2005 was used. The model simulations indicated that the groundwater head converges towards the main river and, finally, to the outlet of the watershed. The study indicated areas of interactions between the river and groundwater. The scenario examination result reveals increasing the present pumping rate by over fifty percent (by 50%, 100%, and 200%) will surely cause visible groundwater head decline near the outlet of the watershed, and substantial river baseflow reduction. The recharge reduction scenario also indicates the huge risk of groundwater sustainability in the area.

The Model of Land Surface Movements Induced by Groundwater Rebound in the Area of Former Mining Exploitation

<p>The objective of the research was to investigate the process of rock mass recompaction related to groundwater rebound induced by underground mining. Research has been conducted in the area of the closed copper ore mine (Konrad) as well as the anhydrite and gypsum mine (Lubichów) in south-eastern Poland.</p><p>The mining operation was carried out in the years 1944-2001 in the area of the Konrad mine and 1944-2015 in the area of the Lubichów mine. It resulted in substantial land subsidence of up to 1.4 m and drainage of the aquifer system. However, it is estimated that the subsidence caused by groundwater pumping during these periods was 0.3 m in total. Furthermore, the spatial extent of the depression cone in the aquifer system immediately after the cessation of exploitation significantly exceeded the limits of the mining areas. Following the closure of the mine, a continuous increase in the groundwater head and land uplift is observed.</p><p>Classical survey results and the Persistent Scatter Satellite Radar Interferometry (PSInSAR) method were used to determine land surface movements in the period from November 2015 to November 2020. The results of the research show in the area of the Lubichów mine closed in June 2015, vertical land uplift reached a maximum of approx. 92 mm in that period. At the same time, in the Konrad mine area, closed in March 2001, no significant land uplift was observed. However, the main part of the investigation concerned the development of a novel method of land uplifting prediction. As a result, an attempt was made to comparatively analyze the dynamics of land uplift associated with the life cycle of the mine and the increase in the groundwater head.</p><p>These analyzes allowed the time factor for the modelling of the land uplift to be determined. This time factor is approx. 5 months in the area of the Lubichów mine and indicates that there is a time lag between the start of the groundwater head increase and the land uplift occurrence. Also, the investigation revealed that land uplift will occur in the analyzed area for the next five years. However, the dynamics of such movements will gradually decline in the years to come.</p><p>The methodology developed could be applied to any post-mining area where groundwater rebound-related uplifts are observed. It may be an appropriate tool for estimating both the time during which the land uplift is expected to begin after the mine drainage has been stopped, as well as the total duration of the land uplift phenomena.</p>

A New Physically-Based Spatially-Distributed Groundwater Flow Module for SWAT+

Watershed models are used worldwide to assist with water and nutrient management under conditions of changing climate, land use, and population. Of these models, the Soil and Water Assessment Tool (SWAT) and SWAT+ are the most widely used, although their performance in groundwater-driven watersheds can sometimes be poor due to a simplistic representation of groundwater processes. The purpose of this paper is to introduce a new physically-based spatially-distributed groundwater flow module called gwflow for the SWAT+ watershed model. The module is embedded in the SWAT+ modeling code and is intended to replace the current SWAT+ aquifer module. The model accounts for recharge from SWAT+ Hydrologic Response Units (HRUs), lateral flow within the aquifer, Evapotranspiration (ET) from shallow groundwater, groundwater pumping, groundwater–surface water interactions through the streambed, and saturation excess flow. Groundwater head and groundwater storage are solved throughout the watershed domain using a water balance equation for each grid cell. The modified SWAT+ modeling code is applied to the Little River Experimental Watershed (LREW) (327 km2) in southern Georgia, USA for demonstration purposes. Using the gwflow module for the LREW increased run-time by 20% compared to the original SWAT+ modeling code. Results from an uncalibrated model are compared against streamflow discharge and groundwater head time series. Although further calibration is required if the LREW model is to be used for scenario analysis, results highlight the capabilities of the new SWAT+ code to simulate both land surface and subsurface hydrological processes and represent the watershed-wide water balance. Using the modified SWAT+ model can provide physically realistic groundwater flow gradients, fluxes, and interactions with streams for modeling studies that assess water supply and conservation practices. This paper also serves as a tutorial on modeling groundwater flow for general watershed modelers.