Sources of hydrogen sulfide in groundwater on reclaimed land

Fisherman Islands is an area of reclaimed land at the mouth of the Brisbane River in Queensland, Australia. Ongoing groundwater monitoring has found elevated concentrations of hydrogen sulfide (H2S) in the groundwater on the island. The presence of H2S on Fisherman Islands is of concern because of its toxic nature, the potential for acid sulfate soil formation, and its noxious odor. The aim of this study was,to, identify the sources of H2S within the groundwater on Fisherman Islands. It was hypothesized that the H2S is being formed by sulfate reducing bacteria acting on sulfate from seawater, rather than the introduction of sulfide with the dredge sediments. Groundwater and soil samples were collected and analyzed for sulfide, sulfate, and organic carbon contents. Elevated concentrations of sulfides coincide with,elevated concentrations of sulfate in the groundwater and elevated concentrations of organic carbon in the sediments, supporting the hypothesis that sulfide formation is the result of heterotrophic, sulfate reducing organisms.

Application of a coupled ground-surface water flow model to simulate periodic groundwater flow influenced by a sloping boundary, capillarity and vertical flows

Theoretical developments as well as field and laboratory data have shown the influence of the capillary fringe on water table fluctuations to increase with the fluctuation frequency. The numerical solution of a full, partially saturated flow equation can be computationally expensive. In this paper, the influence of the capillary fringe on water table fluctuations is simplified through its parameterisation into the storage coefficient of a fully-saturated groundwater flow model using the complex effective porosity concept [Nielsen, P., Perrochet, P., 2000. Water table dynamics under capillary fringes: experiments and modelling. Advances in Water Resources 23 (1), 503-515; Nielsen, P., Perrochet, P., 2000. ERRATA: water table dynamics under capillary fringes: experiments and modelling (Advances in Water Resources 23 (2000) 503-515). Advances in Water Resources 23, 907-908]. The model is applied to sand flume observations of periodic water table fluctuations induced by simple harmonic forcing across a sloping boundary, analogous to many beach groundwater systems. While not providing information on the moisture distribution within the aquifer, this approach can reasonably predict the water table fluctuations in response to periodic forcing across a sloping boundary. Furthermore, he coupled ground-surface water model accurately predicts the extent of the seepage face formed at the sloping boundary. (C) 2005 Elsevier Ltd. All rights reserved.

The cost of uniqueness in groundwater model calibration

Calibration of a groundwater model requires that hydraulic properties be estimated throughout a model domain. This generally constitutes an underdetermined inverse problem, for which a Solution can only be found when some kind of regularization device is included in the inversion process. Inclusion of regularization in the calibration process can be implicit, for example through the use of zones of constant parameter value, or explicit, for example through solution of a constrained minimization problem in which parameters are made to respect preferred values, or preferred relationships, to the degree necessary for a unique solution to be obtained. The cost of uniqueness is this: no matter which regularization methodology is employed, the inevitable consequence of its use is a loss of detail in the calibrated field. This, ill turn, can lead to erroneous predictions made by a model that is ostensibly well calibrated. Information made available as a by-product of the regularized inversion process allows the reasons for this loss of detail to be better understood. In particular, it is easily demonstrated that the estimated value for an hydraulic property at any point within a model domain is, in fact, a weighted average of the true hydraulic property over a much larger area. This averaging process causes loss of resolution in the estimated field. Where hydraulic conductivity is the hydraulic property being estimated, high averaging weights exist in areas that are strategically disposed with respect to measurement wells, while other areas may contribute very little to the estimated hydraulic conductivity at any point within the model domain, this possibly making the detection of hydraulic conductivity anomalies in these latter areas almost impossible. A study of the post-calibration parameter field covariance matrix allows further insights into the loss of system detail incurred through the calibration process to be gained. A comparison of pre- and post-calibration parameter covariance matrices shows that the latter often possess a much smaller spectral bandwidth than the former. It is also demonstrated that, as all inevitable consequence of the fact that a calibrated model cannot replicate every detail of the true system, model-to-measurement residuals can show a high degree of spatial correlation, a fact which must be taken into account when assessing these residuals either qualitatively, or quantitatively in the exploration of model predictive uncertainty. These principles are demonstrated using a synthetic case in which spatial parameter definition is based oil pilot points, and calibration is Implemented using both zones of piecewise constancy and constrained minimization regularization. (C) 2005 Elsevier Ltd. All rights reserved.

Measurement and modelling of controlled beach groundwater levels under wave action

New laboratory scale experimental data are presented on the forcing of beach groundwater levels by wave run-up. The experimental setup simulates a coastal barrier dividing the ocean from a relatively constant back beach water level, conditions approximating a closed off lagoon system or beach aquifer. The data are critically compared to an advanced numerical model for simulating wave and beach groundwater interaction in the coastal zone, and provide the first experimental verification of such a model. Overall model-data comparisons are good, but some systematic discrepancies are apparent, and reasons for these are discussed.

Theory and application of inverse problems in groundwater: numerical, laboratory and field studies.

Aquifers are a vital water resource whose quality characteristics must be safeguarded or, if damaged, restored. The extent and complexity of aquifer contamination is related to characteristics of the porous medium, the influence of boundary conditions, and the biological, chemical and physical processes. After the nineties, the efforts of the scientists have been increased exponentially in order to find an efficient way for estimating the hydraulic parameters of the aquifers, and thus, recover the contaminant source position and its release history. To simplify and understand the influence of these various factors on aquifer phenomena, it is common for researchers to use numerical and controlled experiments. This work presents some of these methods, applying and comparing them on data collected during laboratory, field and numerical tests. The work is structured in four parts which present the results and the conclusions of the specific objectives.

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.

A solar-powered reverse osmosis system for high recovery of freshwater from saline groundwater

Desalination of groundwater is essential in arid regions that are remote from both seawater and freshwater resources. Desirable features of a groundwater desalination system include a high recovery ratio, operation from a sustainable energy source such as solar, and high water output per unit of energy and land. Here we propose a new system that uses a solar-Rankine cycle to drive reverse osmosis (RO). The working fluid such as steam is expanded against a power piston that actuates a pump piston which in turn pressurises the saline water thus passing it through RO membranes. A reciprocating crank mechanism is used to equalise the forces between the two pistons. The choice of batch mode in preference to continuous flow permits maximum energy recovery and minimal concentration polarisation in the vicinity of the RO membrane. This study analyses the sizing and efficiency of the crank mechanism, quantifies energy losses in the RO separation and predicts the overall performance. For example, a system using a field of linear Fresnel collectors occupying 1000 m2 of land and raising steam at 200 °C and 15.5 bar could desalinate 350 m3/day from saline water containing 5000 ppm of sodium chloride with a recovery ratio of 0.7.

Stand-alone groundwater desalination system using reverse osmosis combined with a cooled greenhouse for use in arid and semi-arid zones of India

In many areas of northern India, salinity renders groundwater unsuitable for drinking and even for irrigation. Though membrane treatment can be used to remove the salt, there are some drawbacks to this approach e.g. (1) depletion of the groundwater due to over-abstraction, (2) saline contamination of surface water and soil caused by concentrate disposal and (3) high electricity usage. To address these issues, a system is proposed in which a photovoltaic-powered reverse osmosis (RO) system is used to irrigate a greenhouse (GH) in a stand-alone arrangement. The concentrate from the RO is supplied to an evaporative cooling system, thus reducing the volume of the concentrate so that finally it can be evaporated in a pond to solid for safe disposal. Based on typical meteorological data for Delhi, calculations based on mass and energy balance are presented to assess the sizing and cost of the system. It is shown that solar radiation, freshwater output and evapotranspiration demand are readily matched due to the approximately linear relation among these variables. The demand for concentrate varies independently, however, thus favouring the use of a variable recovery arrangement. Though enough water may be harvested from the GH roof to provide year-round irrigation, this would require considerable storage. Some practical options for storage tanks are discussed. An alternative use of rainwater is in misting to reduce peak temperatures in the summer. An example optimised design provides internal temperatures below 30EC (monthly average daily maxima) for 8 months of the year and costs about €36,000 for the whole system with GH floor area of 1000 m2 . Further work is needed to assess technical risks relating to scale-deposition in the membrane and evaporative pads, and to develop a business model that will allow such a project to succeed in the Indian rural context.