Lattice Boltzmann Modeling of Fluid Flow and Solute Transport in Karst Aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

Optimal aquifers and reservoirs for CCS and EOR in the Kingdom of Saudi Arabia : An overview

Peer reviewed

Table 2 & Appendix A. Properties of operations injecting CO2 into saline aquifers

The experience from CO2 injection at pilot projects (Frio, Ketzin, Nagaoka, US Regional Partnerships) and existing commercial operations (Sleipner, Snøhvit, In Salah, acid-gas injection) demonstrates that CO2 geological storage in saline aquifers is technologically feasible. Monitoring and verification technologies have been tested and demonstrated to detect and track the CO2 plume in different subsurface geological environments. By the end of 2008, approximately 20 Mt of CO2 had been successfully injected into saline aquifers by existing operations. Currently, the highest injection rate and total storage volume for a single storage operation are approximately 1 Mt CO2/year and 25 Mt, respectively. If carbon capture and storage (CCS) is to be an effective option for decreasing greenhouse gas emissions, commercial-scale storage operations will require orders of magnitude larger storage capacity than accessed by the existing sites. As a result, new demonstration projects will need to develop and test injection strategies that consider multiple injection wells and the optimisation of the usage of storage space. To accelerate large-scale CCS deployment, demonstration projects should be selected that can be readily employed for commercial use; i.e. projects that fully integrate the capture, transport and storage processes at an industrial emissions source.

ROCK OUTCROPS IN THE CANADIAN SHIELD: AN INVESTIGATION OF CONTAMINANT TRANSPORT FROM SURFACE SOURCES IN FRACTURED ROCK AQUIFERS

Flow, recharge and transport dynamics in fractured rock aquifers with low lying rock outcrops is a largely unexplored area of study in hydrogeology. The purpose of this thesis is to examine these topics in an agricultural area in Eastern Ontario. The study consists of a regional scale groundwater quality study, an infiltration experiment that considers bacteria transport from the ground surface to a well, and a numerical modelling study that tests the parameters that affect surface infiltration of a tracer from a rock outcrop to a deeper horizontal fracture. In the water quality study, approximately 65% of the samples contained total coliform, 16% contained E. coli, and 1% contained nitrate-N at greater than 5 mg/L. Occurrence of E. coli increased when considering seasonality, where wells were drilled on rock outcrops, and for shallow well intervals. Nitrate-N did not occur above the Guidelines for Canadian Drinking Water Quality (Health Canada, 2012) of 10 mg/L. Rapid arrival times were observed in the infiltration study for both the microspheres (30 minutes) and a dye tracer (45 minutes) in a well approximately 6.0 m in horizontal and 2.8 m in vertical distance from the tracer source. Transport velocities were approximately 38.9 m/day for the dye tracer and 115.2 m/day for the colloidal tracer. Results of the model runs indicate that overburden can provide an effective protective layer to transport in fractures, that high groundwater velocities occur in larger fracture apertures and higher gradients dilute tracer concentrations, and that lower groundwater velocities occur with smaller fracture apertures and lower gradients result in elevated tracer concentrations. Lower rainfall rates, larger fracture apertures, early tracer time, larger gradients, and lower water levels maintained unsaturated conditions for longer time periods such that tracer transport was delayed until saturated conditions were attained. The overall heterogeneity of this aquifer environment creates a source water protection conundrum where the water quality is generally good, while transport can occur very quickly in proximity to rock outcrops and in areas with limited overburden.

The influence of bedrock hydrogeology on catchment-scale nitrate fate and transport in fractured aquifers

Characterising catchment scale biogeochemical processes controlling nitrate fate in groundwater constitutes a fundamental consideration when applying programmes of measures to reduce risks posed by diffuse agricultural pollutants to water quality. Combining hydrochemical analyses with nitrate isotopic data and physical hydrogeological measurements permitted characterisation of biogeochemical processes influencing nitrogen fate and transport in the groundwater in two fractured bedrock aquifers with contrasting hydrogeology but comparable nutrient loads. Hydrochemical and isotopic analyses of groundwater samples collected from moderately fractured, diffusely karstified limestone indicated nitrification controlled dissolved nitrogen fate and delivery to aquatic receptors. By contrast nitrate concentrations in groundwater were considerably lower in a low transmissivity highly lithified sandstone and pyrite-bearing shale unit with patchy subsoil cover. Geophysical and hydrochemical investigations showed shallower intervals contained hydraulically active fractures where denitrification was reflected through lower nitrogen levels and an isotopic enrichment ratio of 1.7 between δ15N and δ18O. Study findings highlight the influence of bedrock hydrogeological conditions on aqueous nitrogen mobility. Investigation results demonstrate that bedrock conditions need to be considered when implementing catchment management plans to reduce the impact of agricultural practices on the quality of groundwater and baseflow in receiving rivers.

Isotope characterisation of deep aquifers in the Gwda catchment, northern Poland

We present the results of isotope measurements (δ 18 O, δ D, δ 13 C DIC and 14 C) and chemical analyses (TDS, TOC, HCO 3– , SO 42– , Cl – , NO 3– , NH 4+ , Ca 2+ , Mg 2+ Na + and K + ) conducted on groundwater samples collected from deep Cenozoic aqui- fers. These aquifers are the basic source of drinking water at numerous localities within the study area in northern Po- land. Most of the δ 18 O determinations are characterised by low variability (i.e., > 70 per cent of δ 18 O are between –9.5‰ and –9.2‰). In most cases tritium activity was not detected or its content slightly exceeded the uncertainty of measure- ment (from ±0.3 T.U. to ± 0.5 T.U.). On average, 14 C activity is twice higher than that under similar conditions and in hydrogeological systems. The δ 13 C DIC values fall within the –13.6‰ to –12.8‰ range. A slight variability is observed when considering all isotope and chemical data within the study area and under these hydrogeological conditions. In general, the results of isotope and chemical analyses seem to be homogeneous, indicating the presence of closely similar groundwaters in the system, irrespective of geological formation. It is likely that there is a significant hydraulic connec- tion between shallow and deep aquifers in the Gwda catchment, which indicates the potential for seepage of pollutants from shallow Pleistocene to deep Miocene aquifers. This can endanger the latter by e.g., high concentrations of NO 3– , SO 42– and Cl – ions from shallow aquifers within the Gwda catchment.

Friction loss correction in flowing well discharge tests

Artesian confined aquifers do not need pumping energy, and water from the aquifer flows naturally at the wellhead. This study proposes correcting the method for analyzing flowing well tests presented by Jacob and Lohman (1952) by considering the head losses due to friction in the well casing. The application of the proposed correction allowed the determination of a transmissivity (T = 411 m(2)/d) and storage coefficient (S = 3 x 10(-4)) which appear to be representative for the confined Guarani Aquifer in the study area. Ignoring the correction due to head losses in the well casing, the error in transmissivity evaluation is about 18%. For the storage coefficient the error is of 5 orders of magnitude, resulting in physically unacceptable value. The effect of the proposed correction on the calculated radius of the cone of depression and corresponding well interference is also discussed.

The introduction of the Salmonella/microsome mutagenicity assay in a groundwater monitoring program

The objective of this study was to verify the possible inclusion of the Salmonella/microsome mutagenicity assay in a groundwater monitoring program as a complementary assay to assess water quality. Groundwater samples belonging to seven wells from different types of aquifers were analyzed. Three different methods for sample preparation were used: membrane filtration; liquid-liquid and XAD-4 extraction. The filtered samples were tested using TA98, TA100, YG1041 and YG1042 and the water extracts only with TA98 and TA100. No mutagenic activity was observed in any of the 16 filtered samples tested. Out of the 10 samples analyzed using XAD-4 extraction, five showed mutagenic activity with potency ranging from 130 to 1500 revertants/L. Concerning the liquid-liquid extraction, from the 11 samples analyzed, 3 showed mutagenicity. The XAD-4 extraction was the most suitable sample preparation. TA98 without S9 was found to be the most sensitive testing condition. The wells presenting water samples with mutagenic activity belonged to unconfined aquifers, which are more vulnerable to contamination. The data suggest that Salmonella/microsome assay can be used as an efficient screening tool to monitor groundwater for mutagenic activity. (C) 2009 Elsevier B.V. All rights reserved.

Submarine groundwater discharge and associated chemical input to a coastal sea

This paper presents a theoretical model of flow and chemical transport processes in subterranean estuaries (unconfined brackish groundwater aquifers at the ocean-land interface). The model shows that groundwater circulation and oscillating flow, caused by wave setup and tide, may constitute up to 96% of submarine groundwater discharge (SGWD) compared with 4% due to the net groundwater discharge. While these local flow processes do not change the total amount of land-derived chemical input to the ocean over a long period (e.g., yearly), they induce fluctuations of the chemical transfer rate as the aquifer undergoes saltwater intrusion. This may result in a substantial increase in chemical fluxes to the ocean over a short period (e.g., monthly and by a factor of 20 above the averaged level), imposing a possible threat to the marine environment. These results are essentially consistent with the experimental findings of Moore [1996] and have important implications for coastal resources management.

Groundwater waves in a coastal aquifer: A new governing equation including vertical effects and capillarity

Groundwater waves, that is, water table fluctuations, are a natural phenomenon in coastal aquifers. They represent an important part of the interaction between the ocean and aquifer and affect the mass exchange between them. This paper presents a new groundwater wave equation. Because it includes the effects of vertical flows and capillarity, the new equation is applicable to both intermediate-depth aquifers and high-frequency waves. Compared with the wave equation derived by Nielsen ed al. [1997], the present equation provides a closer representation of groundwater waves. In particular, it predicts high-frequency water table fluctuations as observed in the field. A validation of the new equation has been carried out by comparing the analytical solutions to it with predictions from direct simulations using the numerical model SUTRA. The effects of various physical parameters and their relative importance are also discussed.

Beach water table fluctuations due to spring-neap tides: moving boundary effects

Tidal water table fluctuations in a coastal aquifer are driven by tides on a moving boundary that varies with the beach slope. One-dimensional models based on the Boussinesq equation are often used to analyse tidal signals in coastal aquifers. The moving boundary condition hinders analytical solutions to even the linearised Boussinesq equation. This paper presents a new perturbation approach to the problem that maintains the simplicity of the linearised one-dimensional Boussinesq model. Our method involves transforming the Boussinesq equation to an ADE (advection-diffusion equation) with an oscillating velocity. The perturbation method is applied to the propagation of spring-neap tides (a bichromatic tidal system with the fundamental frequencies wt and wt) in the aquifer. The results demonstrate analytically, for the first time, that the moving boundary induces interactions between the two primary tidal oscillations, generating a slowly damped water table fluctuation of frequency omega(1) - omega(2), i.e., the spring-neap tidal water table fluctuation. The analytical predictions are found to be consistent with recently published field observations. (C) 2000 Elsevier Science Ltd. All rights reserved.

A two-dimensional analytical solution of groundwater responses to tidal loading in an estuary and ocean

Previous studies on tidal dynamics of coastal aquifers have focussed on the inland propagation of oceanic tides in the cross-shore direction, a configuration that is essentially one-dimensional. Aquifers at natural coasts can also be influenced by tidal waves in nearby estuaries, resulting in a more complex behaviour of head fluctuations in the aquifers. We present an analytical solution to the two-dimensional depth-averaged groundwater flow equation for a semi-infinite aquifer subject to oscillating head conditions at the boundaries. The solution describes the tidal dynamics of a coastal aquifer that is adjacent to a cross-shore estuary. Both the effects of oceanic and estuarine tides on the aquifer are included in the solution. The analytical prediction of the head fluctuations is verified by comparison with numerical solutions computed using a standard finite-difference method. An essential feature of the present analytical solution is the interaction between the cross- and along-shore tidal waves in the aquifer area near the estuary's entry. As the distance from the estuary or coastline increases, the wave interaction is weakened and the aquifer response is reduced, respectively, to the one-dimensional solution for oceanic tides or the solution of Sun (Sun H. A two-dimensional analytical solution of groundwater response to tidal loading in an estuary, Water Resour Res 1997;33:1429-35) for two-dimensional non-interacting tidal waves. (C) 2000 Elsevier Science Ltd. All rights reserved.