Seawater-fresh groundwater natural equilibrium conditions in the aquifers of the Sel&ccedil / uk sub-basin may be disturbed by the discharge occurred during the past 30 years in order to supply water for irrigation and domestic purpose usage. Two dimensional density dependent cross sectional saturated flow and solute transport simulations were carried out to determine whether seawater intrusion has occurred in the Sel&ccedil / uk sub-basin due to the imposed discharge or the determined salt-water in the western section of the aquifers represents natural interface equilibrium conditions. The numerical simulation model was calibrated using field measurements. The results of the simulations suggest that the seawater intrusion has occurred in the study area since the pumping activity increased in the region. Five scenarios were simulated for future predictions: (1) Present recharge and discharge conditions are maintained, (2) Discharge increases at a rate of the municipality need as a result of increasing population, (3) Discharge decreases by 12% from the present value, (4) Discharge decreases by 25% from the present value, and (5) No discharge occurs or recharge amount equals to the present discharge amount condition supplied to the aquifer. These scenario results suggest that seawater intrusion in the Sel&ccedil / uk sub-basin would progress in the following years unless proper management measures are taken into consideration. KEYWORDS: seawater intrusion, density dependent groundwater flow, solute transport, contamination, Sel&ccedil / uk sub-basin.

In this study, The Beach Dewatering System, a relatively recent technology to combat beach erosion, which is proposed as a practical alternative to more traditional shoreline stabilization methods, is investigated and an informative overview on the genesis, development and recent use of this technique is provided. On the basis of the link existing between the elevation of beach groundwater and erosional or accretionary trends at the beach face, a numerical model that simulates groundwater flow in a coastal aquifer under beach drainage is presented. In this model, the seaward boundary of the domain is considered to be tidally fluctuating in a large scale to represent the occurrence of seepage face significantly. The unsteady groundwater flow equation is solved numerically using the method of finite differences. The results clearly showed that the water table being lowered caused the reduction of the seepage face which is the main aim of Beach Dewatering projects. The positional design parameters, i.e. horizontal and vertical location of the drain, are also investigated by utilizing an efficiency index. It is observed that the system efficiency decreased as the drain is shifted landward. The results also indicated that, the efficiency slightly increased with the vertical drain elevation.

The Table Mountain group aquifer system has the potential to be an important supply of water. Although the aquifer system is used to some extent, a number of aspects relating to the aquifer system are poorly understood and unquantified. This study aimed to take into consideration the importance of differenct ecosytems, which is essential in predicting the effects of groundwater abstruction. However, the ecological requirements of systems that depend on groundwater are poorly understood. This project identified &quot / type areas&quot / for further detailed research into the impacts of large-scale groundwater abstraction from the Table Mountain group aquifer system based on the nature and functioning of ecosystems across groundwater dependent ecosystem boundaries of a regional scale.

A quantitative understanding of ground water flow characteristics in unconfined aquifers is important because of the prevalence of abstraction from, and pollution of these systems. The current understanding of ground water flow in unconfined aquifers is limited because of the dominance of horizontal flow modelling strategies used to represent unconfined flow processes. The application of horizontal flow principles leads to an ignorance of seepage-face formation and can not predict the complicated three-dimensional nature of the ground water flow that dominates at the ground water-surface water interface. This study aims to address some of these deficiencies by exploring the true three-dimensional nature of ground water flow including the formation of seepage faces at the ground water-surface water interface using numerical and laboratory techniques. A finite element model for simulating two-dimensional (vertical) variably saturated flow is developed and benchmarked against standard laboratory and field-scale solutions. The numerical features of the finite element model are explored and compared to a simple finite difference formulation. The comparison demonstrates how finite element formulations lead to a broader spatial averaging of material properties and a different method for the representation of specified flux boundaries. A detailed comparison analysis indicates that these differences in the finite element solution lead to an improved approximation to the partial differential equation governing two-dimensional (vertical) variably saturated flow. A laboratory analysis of unconfined ground water flow and associated solute transport characteristics was performed. The analysis focused upon unconfined flow towards a pumping well. The laboratory observations were reliably reproduced using a three-dimensional (axi-symmetric), variably saturated ground water flow model. The model was benchmarked against the ground water flow characteristics such as the seepage-face height and total flow rate. In addition, the model was shown to reliably reproduce the solute transport features such as travel times and streamline distributions. This is the first time that a numerical model has been used to reliably reproduce the solute transport characteristics near a seepage-face boundary where the three-dimensional flow effects are prevalent. The ability to reliably predict solute transport patterns in the seepage-face zone is important since this region is known to support vital microbially facilitated reactions that control nutrient cycling and contaminant attenuation. The three-dimensional travel time distribution near the seepage-face was compared to that predicted using a horizontal flow modelling approach derived from the basic Dupuit-Forchheimer equations. The Dupuit-Forchheimer based model indicated that horizontal flow modelling would under-estimate the total residence time near a seepage-face boundary, thereby introducing a considerable source of error in a solute transport analysis. For this analysis, a new analytical solution for the steady travel time distribution in an unconfined aquifer subject to a single pumping well was derived. The analytical model has identified, for the first time in the hydrogeology literature, the use of the imaginary error function. The imaginary error function is a standard transcendental function and an infinite series approach to evaluate the function was successfully proposed. The two-dimensional (vertical) ground water flow model was extended to handle the case where the flow is driven by density gradients near the ground water-surface water interface. The unsteady, two-dimensional, Galerkin finite element model of density-dependent ground water flow in variably saturated porous media is rigorously presented and partially benchmarked under fully saturated (confined) conditions. The partial benchmarking involved reproducing solutions to the standard Henry salt-water intrusion and the Elder salt-convection problems. The model was used in a standard density-coupled and a new density-uncoupled mode to elucidate the worthiness of the Henry and Elder problems as benchmark standards. A comparison of the coupled and uncoupled solutions indicates that the Henry salt-water intrusion problem has limited worthiness as a benchmark as the patterns of ground water flow are relatively insensitive to density-coupled effects. Alternatively, the Elder problem is completely dependent upon a correct representation of the density-coupled flow and solute transport processes. The coupled versus uncoupled comparison is proposed as a new test of the worthiness of benchmark standards. The Henry salt-water intrusion problem was further analysed in an attempt to alleviate some of the difficulties associated with this benchmark problem. The numerical model was tested against a re-evaluated version of Henry's semi-analytical solution for the coupled solute concentration distribution. The numerical model was used to propose a modified version of the Henry problem where the importance of density-coupled processes was increased. The modified problem was shown to have an improved worthiness as compared to the standard solution. The numerical model results were benchmarked against a new set of semi-analytical results for the modified problem. Certain advantages in using the modified problem as a test case for benchmarking the results of a numerical model of density-dependent ground water flow are identified. A numerical investigation of the patterns of density-driven ground water flow at the ground water-surface water interface was undertaken. The numerical model is shown to produce grid-independent results for a finely discretised domain. The pattern of discharge is controlled, in part, by two parameters. One describes the recharge applied to the aquifer, and the second describes the magnitude of the density differences between the fresh recharging fluid and the saline receiving fluid. The influence of dense intrusions upon the formation of seepage-face boundaries at the ground water-surface water interface under steady-state conditions was also investigated. Dense intrusions are shown to dominate the pattern of ground water flow only under mild recharge conditions, while seepage faces dominate the outflow pattern under strong recharge conditions. Therefore, the formation of seepage-face boundaries and dense intrusions are unlikely to coincide under the conditions examined in this study.

Groundwater studies face computational limitations when providing local detail within regional models. The researchers are concentrated on applying the numerical models to minimize the difference between the physical reality and the implemented numerical model by considering the minimum computational cost. This work consists of the study of line-elements (such as line-doublets, circles, polygons, fractures) using the Analytic Element Method (AEM) for groundwater flow. In this work, we consider the study of two-dimensional groundwater flow in fractured porous media by the Analytic Element Method. We develop a numerical solution based on a series expansion for a problem with more than one fracture. Each fracture has an influence that can be expanded in a series that satisfies Laplaces equation exactly. In the series expansion, the unknown coefficients are obtained from the discharge potentials of all other elements that are related to the expansion coefficients. Sizes, locations and conductivities for all inhomogeneities are selected arbitrarily. This work also discusses a matrix method obtained by imposing the intern boundary conditions for the Analytic Element Method. The convergence analysis of a Gauss-Seidel type iterative method is also discussed. / Estudos de águas subterrâneas enfrentam limitações computacionais ao fornecer detalhes locais em modelos regionais. Os pesquisadores estão concentrados na aplicação dos modelos numéricos para minimizar a diferença entre a realidade física e o modelo numérico implementado considerando o custo computacional mínimo. Este trabalho consiste no estudo de elementos de linha (como line-doublets, círculos, polígonos, fraturas) usando o Método de Elemento Analítico (AEM) para o fluxo de águas subterrâneas. Neste trabalho, consideramos o estudo do fluxo bidimensional de águas subterrâneas em meios porosos fraturados pelo Método dos Elementos Analíticos. Desenvolvemos uma solução numérica baseada em uma expansão em série para um problema com mais de uma fratura. Cada fratura tem uma influência que pode ser expandida em uma série que satisfaça exatamente a equação de Laplace. Na expansão da série, os coeficientes desconhecidos são obtidos a partir dos potenciais de descarga de todos os outros elementos que estão relacionados aos coeficientes de expansão. Tamanhos, locais e condutividades para todas as não-homogeneidades são arbitrariamente selecionados. Este trabalho também discute o método da matriz obtido impondo as condições de contorno do interno para o Método do Elemento Analítico. A análise de convergência de um método iterativo tipo Gauss-Seidel também é discutida.

No description available.

This project implements a real application on the Spider II, which is a simulation of Soil Water Movement Model. The main objectives of this project were to develop a parallel and distributed algorithm for the Soil Water Model; implement the Soil Water Movement Simulation model on the Spider II distributed system and to evaluate the performance of simulating the Soil Water Movement Model on Spider II.

Masters of Science / The thesis investigates the potential of the Rietfontein and Spitskop aquifer to meet a demand of 1000 m3/d (12 Lis) as an alternative water resource for the Rietfontein and Spitskop Community. Increasing demand for clean and hygienic drinking water puts more pressure on one of our most valuable resources and supplying all communities with surface water is an extremely difficult and costly task in rural areas like Rietfontein and Spitskop in South Africa. Therefore it is necessary that interim water supplies be found from local aquifers and be utilized to address water supply challenges. Groundwater may serve as a short-term and an interim water supply which may be useful during future dry periods. Abstraction of groundwater is sensitive to recharge. Due to semi-arid conditions in Rietfontein and Spitskop area, there is high rainfall variation and disparity each year. During the research, hydrocensus was carried out. Water samples for chemistry analysis were taken. Literature review and pumping test data was utilized from the previous studies done by different consultants (Botha, 2000, Vivier and Pretorius, 2003). A numerical groundwater flow model for the local aquifers in the area and surroundings was constructed focusing on recharge and abstraction scenarios for the water supply from the local aquifers. For Rietfontein and Spitskop area, the mean annual precipitation (MAP) is 617mm/year. To be assured and rational in determining aquifers ability to meet the required demand a recharge with ninety-five (95th ) percentile was recommended, which estimates the MAP of 308mm/year which is 50% lower than the average MAP of 617mm/year. For a period of ninety six (96) years, the data indicates a severe drought occurred four (4) times where the rainfall was even lower than 95th percentile level of assurance of recharge estimated. This gives a comprehensible indication that average mean is not ideal or steadfast stature when building a water supply numerical groundwater flow modelling. These aquifers can only be exploited and managed if a reliable method can be obtained to estimate their long-term sustainable potential, since the sustainable potential of these aquifers to supply the communities is dependent upon the recharge from rainfall. The results from a numerical groundwater flow model indicated that a combined potential from the local aquifer from eighteen (18) boreholes is sufficient to meet the required demand and a total of 2600m3 Id can potentially be abstracted from the aquifer. With ninety-five (95th ) percentile recharge rate of 308mm/year a numerical groundwater flow model shows that the rate of abstraction is 80% far less than recharge, which gives high level of assurance in terms of local aquifer water supply demand. The abstraction of the boreholes confirmed by a numerical groundwater flow model shows the least impact on the surrounding aquifer for an extended period of time. In the event of drought, the boreholes will see a decline in water level after two (2) months of pumping local aquifer. The water level will decline steadily from two (2) months to two (2) years with a change in water levels of up to 40m. The impact of the drought is minimal compared to recharge rate, which verifies less depletion of the aquifer. The local aquifer shows the potential of 3MLld can be supplied to the communities with an assurance level of 95th percentile of rainfall. Reliable quantification of groundwater recharge rate remains the main challenges the hydrologist experienced and further research is essential for improvement of groundwater management for the area concerned.

The middle Holocene cataclysmic eruption of Mount Mazama blanketed Walker Rim, in south central Oregon, with 270 cm to 300 cm of pumice, causing capture of surface water systems by groundwater, stream relocation, and the formation of biologically diverse fens and seasonal wetlands. The pumice aquifer at Round Meadow, an 8.6 km2 basin, hosts both a fen and seasonally ponded wetlands. The Round Meadow watershed lies within a closed basin between the upper Klamath and Deschutes river basins. As the highest meadow at Walker Rim, it is a relatively well-constrained system to study the effects of hydrological disruption. A water budget was calculated for the basin, hydraulic conductivity was evaluated for the three main sediment layers in the meadow, recharge sources and evaporative trends were studied using stable isotope analysis, and aquifer residence times were estimated using CFC tracer water age dating. Water year 2014 was a drought year and observation of the system under stressed conditions allowed discrimination of four independently functioning components of the hydrogeologic system. These were the meadow, which is by far the largest component in terms of water storage, the fen where iron cementation and up to 1 m of peat holds water in a berm above the meadow, three springs which are sourced from deeper groundwater hosted in the bedrock which underlies the pumice deposit, and the outflow area. In all cases, the aquifer material is pumice, but the influence of the pre-eruption landscape and post-eruption modifications of the aquifer material have resulted in partial isolation of the components. The water budget analysis indicated that the basin lost 44 cm of water storage during WY 2014. Hydraulic conductivity values of 1x10-6, 2x10-2, and 4x10-5 cm/s, were determined for the diatomaceous silt underlying the pumice, the Plinian pumice fall aquifer, and for the diatomaceous silt overlying the pumice, respectively. The pumice is characterized as a perched, weakly confined aquifer and residence times in the pumice are much longer (decades) than for water near the surface of the meadow. Water discharging at the springs is isotopically different (lighter) than either the surface water or groundwater in the pumice aquifer. The fen at Round Meadow appears dependent on seasonal precipitation to recharge water, and responds to fluctuations in annual precipitation. The wetland meadows are volumetrically the main water-storing features at Round Meadow, and these are not homogenous features, but a combination of discrete components.

Groundwater interactions at a regional scale are of great importance to characterize subsurface flow processes. Extensive researches have been conducted previously to determine the main factors controlling the regional implications on groundwater flux circulation. Groundwater circulation occurs due to variation in the groundwater table (hydraulic gradient) across the spatial scale. Previous research highlighted the correlation between groundwater table with both topography variation and the recharge from precipitation. This study aims to highlight the impact of these boundary conditions. Five catchments located across different regions of Sweden with different topographical, hydrological, and meteorological properties considered for this study: Bodalsån, Forsmarksån, Tullstorpsån, Sävaån, and Krycklan. Relevant data were collected and numerical models were set up in steady- state conditions for each of these catchments, using 3D Multiphysics COMSOL. Models were set up for both of the boundary conditions, using 10 m grid resolution.Groundwater flux profiles along the depth of the catchments were obtained as a result, in which significant differences were observed. This was associated predominantly with the difference in the nature of the topography, the slope and soil permeability in these regions. The data thus collected and the models so established have increased the understanding in these regions from a research perspective. / En ökad förståelse för hur grundvatten interagerar med ytvatten är av stor betydelse för att karakterisera underjordiska flödesprocesser. Omfattande undersökningar har tidigare genomförts för att bestämma de viktigaste faktorerna som styr de regionala konsekvenserna för cirkulationen av grundvattenflöde. Grundvattencirkulation uppstår på grund av variationer i grundvattentabellen (hydraulisk gradient) över den rumsliga skalan. Tidigare forskning belyste korrelationen mellan grundvattentabellen med både topografivariation och laddning från nederbörd. Denna studie syftar till att belysa effekterna av dessa gränsförhållanden. Fem avrinningsområden placerade över olika regioner i Sverige med olika topografiska, hydrologiska och meteorologiska egenskaper som beaktas för denna studie: Bodalsån, Forsmarksån, Tullstorpsån, Sävaån och Krycklan. Relevanta data samlades in och numeriska modeller sattes upp under steady-state-förhållanden för vart och ett av dessa avrinningsområden med användning av 3D Multiphysics COMSOL. Modeller konfigurerades för båda gränsförhållandena med 10 m nätupplösning. Grundvattenflödesprofiler längs avrinningsdjupet erhölls som resultat, där signifikanta skillnader observerades. Detta var främst förknippat med skillnaden i topografins natur och lutningen i dessa regioner. De data som samlats in och de så etablerade modellerna har ökat skapat en värdefull grund för vidare hydrologisk forskning i dessa regioner.