Driving mechanisms for groundwater flow and salt transport in a subterranean estuary

This paper presents field measurements and numerical simulations of groundwater dynamics in the intertidal zone of a sandy meso-tidal beach. The study, focusing on vertical hydraulic gradients and pore water salinities, reveals that tides and waves provide important forcing mechanisms for flow and salt transport in the nearshore aquifer. Such forcing, interacting with the beach morphology, enhances the exchange between the aquifer and ocean. The spatial and temporal variations of vertical hydraulic gradients demonstrate the complexity and dynamic nature of the processes and the extent of mixing between fresh groundwater and seawater in a subterranean estuary''. These results provide evidence of a potentially important reaction zone in the nearshore aquifer driven by oceanic oscillations. Land-derived contaminants may undergo important biogeochemical transformations in this zone prior to discharge.

Three-dimensional model for multi-component reactive transport with variable density groundwater flow

PHWAT is a new model that couples a geochemical reaction model (PHREEQC-2) with a density-dependent groundwater flow and solute transport model (SEAWAT) using the split-operator approach. PHWAT was developed to simulate multi-component reactive transport in variable density groundwater flow. Fluid density in PHWAT depends not on only the concentration of a single species as in SEAWAT, but also the concentrations of other dissolved chemicals that can be subject to reactive processes. Simulation results of PHWAT and PHREEQC-2 were compared in their predictions of effluent concentration from a column experiment. Both models produced identical results, showing that PHWAT has correctly coupled the sub-packages. PHWAT was then applied to the simulation of a tank experiment in which seawater intrusion was accompanied by cation exchange. The density dependence of the intrusion and the snow-plough effect in the breakthrough curves were reflected in the model simulations, which were in good agreement with the measured breakthrough data. Comparison simulations that, in turn, excluded density effects and reactions allowed us to quantify the marked effect of ignoring these processes. Next, we explored numerical issues involved in the practical application of PHWAT using the example of a dense plume flowing into a tank containing fresh water. It was shown that PHWAT could model physically unstable flow and that numerical instabilities were suppressed. Physical instability developed in the model in accordance with the increase of the modified Rayleigh number for density-dependent flow, in agreement with previous research. (c) 2004 Elsevier Ltd. All rights reserved.

Imaging and quantifying salt-tracer transport in a riparian groundwater system by means of 3D ERT monitoring

Determining groundwater flow paths of infiltrated river water is necessary for studying biochemical processes in the riparian zone, but their characterization is complicated by strong temporal and spatial heterogeneity. We investigated to what extent repeat 3D surface electrical resistance tomography (ERT) can be used to monitor transport of a salt-tracer plume under close to natural gradient conditions. The aim is to estimate groundwater flow velocities and pathways at a site located within a riparian groundwater system adjacent to the perialpine Thur River in northeastern Switzerland. Our ERT time-lapse images provide constraints on the plume's shape, flow direction, and velocity. These images allow the movement of the plume to be followed for 35 m. Although the hydraulic gradient is only 1.43 parts per thousand, the ERT time-lapse images demonstrate that the plume's center of mass and its front propagate with velocities of 2x10(-4) m/s and 5x10(-4) m/s, respectively. These velocities are compatible with groundwater resistivity monitoring data in two observation wells 5 m from the injection well. Five additional sensors in the 5-30 m distance range did not detect the plume. Comparison of the ERT time-lapse images with a groundwater transport model and time-lapse inversions of synthetic ERT data indicate that the movement of the plume can be described for the first 6 h after injection by a uniform transport model. Subsurface heterogeneity causes a change of the plume's direction and velocity at later times. Our results demonstrate the effectiveness of using time-lapse 3D surface ERT to monitor flow pathways in a challenging perialpine environment over larger scales than is practically possible with crosshole 3D ERT.

Zum Transport von partikelgebundenen Schadstoffen in der ungesättigten Zone

Bei einer Risikoabschätzung bezüglich einer Gefährdung des Schutzgutes Grundwasser müssen alle relevanten Transportpfade, auf denen Schadstoffe durch die Bodenzone bis ins Grundwasser verlagert werden, identifiziert und quantifiziert werden. Die Verlagerung von Schadstoffen gebunden an mobile Partikel im Sickerwasser wird dabei oft vernachlässigt. In dieser Arbeit wurden sowohl experimentelle Untersuchungen zum Partikeltransport in der Bodenzone als auch Szenarienmodellierungen hinsichtlich der Wechselwirkung Partikel/Schadstoff durchgeführt. Die experimentellen ungesättigten Säulenversuche wurden unter naturnahen stationären und instationären hydraulischen und hydrochemischen Bedingungen durchgeführt. Dabei wurde der Einfluss der Parameter Durchmesser Bodenmatrix, Partikelgröße, Beregnungsintensität, Oberflächenspannung und Hydrochemie auf den Transport von natürlichen und synthetischen Partikeln untersucht. Des Weiteren wurden Untersuchungen zur partikelgebundenen Verlagerung von Phenanthren durchgeführt. In einer numerischen Szenarienmodellierung mit dem Modell SMART wurde untersucht, unter welchen Randbedingungen der Transport von Partikeln gleichzeitig zu signifikanten partikelgebundenen Schadstoffkonzentrationen im Grundwasser führt. Dabei wurden die Parameter Lithologie Partikel/Boden, Hydrophobizität Schadstoff, Partikelkonzentration, Partikeldurchmesser sowie Körnung Bodenmatrix variiert. Die Ergebnisse dieser Arbeit zeigen, dass der partikelgebundene Schadstofftransportpfad in der ungesättigten Bodenzone in verschiedenen Szenarien den Anteil mobiler Schadstoffe, die mit dem Sickerwasser ins Grundwasser gelangen, signifikant erhöht. Auf Basis der experimentellen und theoretischen Untersuchungen wurde ein zweistufiges Bewertungsschema entwickelt, das bereits im Vorfeld einer Risikoabschätzung als Entscheidungshilfe hinsichtlich der Relevanz einer Mobilisierung, eines Transports und des Rückhalts von partikelgebundenen Schadstoffen in der ungesättigten Zone dient.

Model reduction of stochastic groundwater flow and transport processes

This work presents a comprehensive methodology for the reduction of analytical or numerical stochastic models characterized by uncertain input parameters or boundary conditions. The technique, based on the Polynomial Chaos Expansion (PCE) theory, represents a versatile solution to solve direct or inverse problems related to propagation of uncertainty. The potentiality of the methodology is assessed investigating different applicative contexts related to groundwater flow and transport scenarios, such as global sensitivity analysis, risk analysis and model calibration. This is achieved by implementing a numerical code, developed in the MATLAB environment, presented here in its main features and tested with literature examples. The procedure has been conceived under flexibility and efficiency criteria in order to ensure its adaptability to different fields of engineering; it has been applied to different case studies related to flow and transport in porous media. Each application is associated with innovative elements such as (i) new analytical formulations describing motion and displacement of non-Newtonian fluids in porous media, (ii) application of global sensitivity analysis to a high-complexity numerical model inspired by a real case of risk of radionuclide migration in the subsurface environment, and (iii) development of a novel sensitivity-based strategy for parameter calibration and experiment design in laboratory scale tracer transport.

THE TRANSPORT OF HEAT BY FLOWING GROUNDWATER IN A DISCRETE FRACTURE

In typical theoretical or experimental studies of heat migration in discrete fractures, conduction and thermal dispersion are commonly neglected from the fracture heat transport equation, assuming heat conduction into the matrix is predominant. In this study analytical and numerical models are used to investigate the significance of conduction and thermal dispersion in the plane of the fracture for a point and line sources geometries. The analytical models account for advective, conductive and dispersive heat transport in both the longitudinal and transverse directions in the fracture. The heat transport in the fracture is coupled with a matrix equation in which heat is conducted in the direction perpendicular to the fracture. In the numerical model, the governing heat transport processes are the same as the analytical models; however, the matrix conduction is considered in both longitudinal and transverse directions. Firstly, we demonstrate that longitudinal conduction and dispersion are critical processes that affect heat transport in fractured rock environments, especially for small apertures (eg. 100 μm or less), high flow rate conditions (eg. velocity greater than 50 m/day) and early time (eg. less than 10 days). Secondly, transverse thermal dispersion in the fracture plane is also observed to be an important transport process leading to retardation of the migrating heat front particularly at late time (eg. after 40 days of hot water injection). Solutions which neglect dispersion in the transverse direction underestimate the locations of heat fronts at late time. Finally, this study also suggests that the geometry of the heat sources has significant effects on the heat transport in the system. For example, the effects of dispersion in the fracture are observed to decrease when the width of the heat source expands.

Quantitative risk assessment of groundwater quality utilizing GIS technology and coupled groundwater models

The thesis presents a two-dimensional Risk Assessment Method (RAM) where the assessment of risk to the groundwater resources incorporates both the quantification of the probability of the occurrence of contaminant source terms, as well as the assessment of the resultant impacts. The approach emphasizes the need for a greater dependency on the potential pollution sources, rather than the traditional approach where assessment is based mainly on the intrinsic geo-hydrologic parameters. The risk is calculated using Monte Carlo simulation methods whereby random pollution events were generated to the same distribution as historically occurring events or a priori potential probability distribution. Integrated mathematical models then simulate contaminant concentrations at the predefined monitoring points within the aquifer. The spatial and temporal distributions of the concentrations were calculated from repeated realisations, and the number of times when a user defined concentration magnitude was exceeded is quantified as a risk. The method was setup by integrating MODFLOW-2000, MT3DMS and a FORTRAN coded risk model, and automated, using a DOS batch processing file. GIS software was employed in producing the input files and for the presentation of the results. The functionalities of the method, as well as its sensitivities to the model grid sizes, contaminant loading rates, length of stress periods, and the historical frequencies of occurrence of pollution events were evaluated using hypothetical scenarios and a case study. Chloride-related pollution sources were compiled and used as indicative potential contaminant sources for the case study. At any active model cell, if a random generated number is less than the probability of pollution occurrence, then the risk model will generate synthetic contaminant source term as an input into the transport model. The results of the applications of the method are presented in the form of tables, graphs and spatial maps. Varying the model grid sizes indicates no significant effects on the simulated groundwater head. The simulated frequency of daily occurrence of pollution incidents is also independent of the model dimensions. However, the simulated total contaminant mass generated within the aquifer, and the associated volumetric numerical error appear to increase with the increasing grid sizes. Also, the migration of contaminant plume advances faster with the coarse grid sizes as compared to the finer grid sizes. The number of daily contaminant source terms generated and consequently the total mass of contaminant within the aquifer increases in a non linear proportion to the increasing frequency of occurrence of pollution events. The risk of pollution from a number of sources all occurring by chance together was evaluated, and quantitatively presented as risk maps. This capability to combine the risk to a groundwater feature from numerous potential sources of pollution proved to be a great asset to the method, and a large benefit over the contemporary risk and vulnerability methods.

An integrated flow and transport model to study the impact of mercury remediation strategies for East Fork Poplar Creek watershed, Oak Ridge, Tennessee

An integrated flow and transport model using MIKE SHE/MIKE 11 software was developed to predict the flow and transport of mercury, Hg(II), under varying environmental conditions. The model analyzed the impact of remediation scenarios within the East Fork Poplar Creek watershed of the Oak Ridge Reservation with respect to downstream concentration of mercury. The numerical simulations included the entire hydrological cycle: flow in rivers, overland flow, groundwater flow in the saturated and unsaturated zones, and evapotranspiration and precipitation time series. Stochastic parameters and hydrologic conditions over a five year period of historical hydrological data were used to analyze the hydrological cycle and to determine the prevailing mercury transport mechanism within the watershed. Simulations of remediation scenarios revealed that reduction of the highly contaminated point sources, rather than general remediation of the contaminant plume, has a more direct impact on downstream mercury concentrations.

An Integrated Flow and Transport Model to Study the Impact of Mercury Remediation Strategies for East Fork Poplar Creek Watershed, Oak Ridge, Tennessee

An integrated flow and transport model using MIKE SHE/MIKE 11 software was developed to predict the flow and transport of mercury, Hg(II), under varying environmental conditions. The model analyzed the impact of remediation scenarios within the East Fork Poplar Creek watershed of the Oak Ridge Reservation with respect to downstream concentration of mercury. The numerical simulations included the entire hydrological cycle: flow in rivers, overland flow, groundwater flow in the saturated and unsaturated zones, and evapotranspiration and precipitation time series. Stochastic parameters and hydrologic conditions over a five year period of historical hydrological data were used to analyze the hydrological cycle and to determine the prevailing mercury transport mechanism within the watershed. Simulations of remediation scenarios revealed that reduction of the highly contaminated point sources, rather than general remediation of the contaminant plume, has a more direct impact on downstream mercury concentrations.

An investigation of localised soil heterogeneities on solute transport using a multisegment percolation system

A multisegment percolation system (MSPS) consisting of 25 individual collection wells was constructed to study the effects of localised soil heterogeneities on the transport of solutes in the vadose zone. In particular, this paper discusses the transport of water and nutrients (NO3-, Cl-, PO43-) through structurally stable, free-draining agricultural soil from Victoria, Australia. A solution of nutrients was irrigated onto the surface of a large undisturbed soil core over a 12-h period. This was followed by a continuous irrigation of distilled water at a fate which did not cause pending for a further 18 days. During this time, the volume of leachate and the concentration of nutrients in the leachate of each well were measured. Very significant variation in drainage patterns across a small spatial scale was observed. Leaching of nitrate-nitrogen and chloride from the core occurred two days after initial application. However, less than 1% of the total applied phosphate-phosphorus leached from the soil during the 18-day experiment, indicating strong adsorption. Our experiments indicate considerable heterogeneity in water flow patterns and solute leaching on a small spatial scale. These results have significant ramifications for modelling solute transport and predicting nutrient loadings on a larger scale.

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.

Modelling solute transport in structured soils: Performance evaluation of the ADR and TRM models

The movement of chemicals through the soil to the groundwater or discharged to surface waters represents a degradation of these resources. In many cases, serious human and stock health implications are associated with this form of pollution. The chemicals of interest include nutrients, pesticides, salts, and industrial wastes. Recent studies have shown that current models and methods do not adequately describe the leaching of nutrients through soil, often underestimating the risk of groundwater contamination by surface-applied chemicals, and overestimating the concentration of resident solutes. This inaccuracy results primarily from ignoring soil structure and nonequilibrium between soil constituents, water, and solutes. A multiple sample percolation system (MSPS), consisting of 25 individual collection wells, was constructed to study the effects of localized soil heterogeneities on the transport of nutrients (NO3-, Cl-, PO43-) in the vadose zone of an agricultural soil predominantly dominated by clay. Very significant variations in drainage patterns across a small spatial scale were observed tone-way ANOVA, p < 0.001) indicating considerable heterogeneity in water flow patterns and nutrient leaching. Using data collected from the multiple sample percolation experiments, this paper compares the performance of two mathematical models for predicting solute transport, the advective-dispersion model with a reaction term (ADR), and a two-region preferential flow model (TRM) suitable for modelling nonequilibrium transport. These results have implications for modelling solute transport and predicting nutrient loading on a larger scale. (C) 2001 Elsevier Science Ltd. All rights reserved.

Tidal effects on groundwater dynamics in unconfined aquifers

The variation of seawater level resulting from tidal fluctuations is usually neglected in regional groundwater flow studies. Although the tidal oscillation is damped near the shoreline, there is a quasi-steady-slate rise in the mean water-table position, which may have an influence on regional groundwater flow. In this paper the effects of tidal fluctuations on groundwater hydraulics are investigated using a variably saturated numerical model that includes the effects of a realistic mild beach slope, seepage face and the unsaturated zone. In particular the impact of these factors on the velocity field in the aquifer is assessed. Simulations show that the tidal fluctuation has substantial consequences for the local velocity field in the vicinity of the exit face, which affects the nearshore migration of contaminant in coastal aquifers. An overheight in the water table as a result of the tidal fluctuation is observed anti this has a significant effect on groundwater discharge to the sea when the landward boundary condition is a constant water level. The effect of beach slope is very significant and simplifying the problem by considering a vertical beach face causes serious errors in predicting the water-table position and the groundwater flux. For media with a high effective capillary fringe, the moisture retained above the water table is important in determining the effects of the tidal fluctuations. Copyright (C) 2001 John Wiley & Sons, Ltd.

Project, modelling and simulation of air stripping columns for the treatment of groundwater contaminated with volatile organic compounds

Volatile organic compounds are a common source of groundwater contamination that can be easily removed by air stripping in columns with random packing and using a counter-current flow between the phases. This work proposes a new methodology for the column design for any particular type of packing and contaminant avoiding the necessity of a pre-defined diameter used in the classical approach. It also renders unnecessary the employment of the graphical Eckert generalized correlation for pressure drop estimates. The hydraulic features are previously chosen as a project criterion and only afterwards the mass transfer phenomena are incorporated, in opposition to conventional approach. The design procedure was translated into a convenient algorithm using C++ as programming language. A column was built in order to test the models used either in the design or in the simulation of the column performance. The experiments were fulfilled using a solution of chloroform in distilled water. Another model was built to simulate the operational performance of the column, both in steady state and in transient conditions. It consists in a system of two partial non linear differential equations (distributed parameters). Nevertheless, when flows are steady, the system became linear, although there is not an evident solution in analytical terms. In steady state the resulting system of ODE can be solved, allowing for the calculation of the concentration profile in both phases inside the column. In transient state the system of PDE was numerically solved by finite differences, after a previous linearization.