Simulation of benzene transport and biodegradation during transient hydraulic conditions

Thesis (M.S.) University of Alaska Fairbanks, 2000MODFLOW and BIOMOC were used to simulate transport and biodegradation of benzene in the alluvial aquifer adjacent to the Chena River. MODFLOW was used to calculate ground water fluxes at the boundaries of the BIOMOC model, which was used to model transport and biodegradation of benzene. A benzene plume located 300 ft. southeast of the study site is superimposed onto the cross-sectional model of the study area. Only saturated zone processes were modeled. Anaerobic biodegradation was the only simulated biodegradation process. The simulation shows 0.003% of the theoretical benzene entering the saturated zone is biodegraded, 0.6% is adsorbed by solids, and 99.4% leaves the model boundaries. The simulation predicts theoretical concentrations of benzene are 2 to 8 ug/l when discharging into the river. Field data do not support this finding. Processes not simulated, such as aerobic degradation at the water table, may make significant contributions toward limiting benzene transport

A Conceptualized Groundwater Flow Model Development for Integration with Surface Hydrology Model

A groundwater system model was developed and calibrated in the study area of Lehman Creek watershed, eastern Nevada. The model development aims for integrating the surface hydrologic model - precipitation runoff modeling system (PRMS) model - with the three-dimensional (3D) finite-difference model MODFLOW. A two-layer groundwater model was developed with spatial discretization of 100 x 100 m grid. The water balance was estimated with inflows of gravity drainage and initial streamflow estimated from a calibrated PRMS model, and with outflows of spring discharges, boundary fluxes, and stream base flow. A steady-state model calibration was performed to estimate the hydraulic properties. The modeling results were able to represent the geographic relieves, simulate water balance components, and capture the hydrogeologic features. The preliminary results presented in this study provide insights into the local groundwater flow system and lay groundwork for future study of interactive influences of surface hydrologic variation

Error Control of Iterative Linear Solvers for Integrated Groundwater Models

An open problem that arises when using modern iterative linear solvers, such as the preconditioned conjugate gradient (PCG) method or Generalized Minimum RESidual method (GMRES) is how to choose the residual tolerance in the linear solver to be consistent with the tolerance on the solution error. This problem is especially acute for integrated groundwater models which are implicitly coupled to another model, such as surface water models, and resolve both multiple scales of flow and temporal interaction terms, giving rise to linear systems with variable scaling. This article uses the theory of 'forward error bound estimation' to show how rescaling the linear system affects the correspondence between the residual error in the preconditioned linear system and the solution error. Using examples of linear systems from models developed using the USGS GSFLOW package and the California State Department of Water Resources' Integrated Water Flow Model (IWFM), we observe that this error bound guides the choice of a practical measure for controlling the error in rescaled linear systems. It is found that forward error can be controlled in preconditioned GMRES by rescaling the linear system and normalizing the stopping tolerance. We implemented a preconditioned GMRES algorithm and benchmarked it against the Successive-Over-Relaxation (SOR) method. Improved error control reduces redundant iterations in the GMRES algorithm and results in overall simulation speedups as large as 7.7x. This research is expected to broadly impact groundwater modelers through the demonstration of a practical approach for setting the residual tolerance in line with the solution error tolerance.Comment: 13 pages and 1 figur

An example of aquifer heterogeneity simulation to modeling well-head protection areas

Groundwater management requires the definition of Well-Head Protection Areas (WHPA) for water supply wells. Italian law uses geometrical, chronological and hydrogeological criteria for WHPA identification, providing a groundwater travel time of 60 days for the definition of the Zone of Travel (ZOT). An exhaustive ZOT delineation must involve numerical modeling of groundwater flow together with simulation of the advective component of the transport process. In this context, the spatial variability of hydrogeological and transport parameters has to be critically estimated during numerical modeling implementation. In the present article, geostatistical simulation using a transition probability approach and groundwater numerical modeling were performed to delineate WHPAs for several supply wells in the middle Venetian Plain, taking into account the lithologic heterogeneity of the aquifer. The transition probability approach for the lithologic data was developed by T-PROGS software, while MODDLOW-2005 and PEST-ASP were used, respectively, to reproduce and calibrate site-specific hydraulic head data. Finally, a backward particle tracking analysis was performed with MODPATH to outline the 60-day ZOT

Parsimonious numerical modelling of deep geothermal reservoirs

Numerical modelling has been undertaken to help improve understanding of a deep geothermal system being considered for development in the vicinity of Eastgate (Weardale, County Durham, UK). A parsimonious numerical modelling approach is used, which allows the possibility to develop a workable formal framework, rigorously testing evolving concepts against data as they become available. The approach used and results presented in this study are valuable as a contribution to a wider understanding of deep geothermal systems. This modelling approach is novel in that it utilises the mass transport code MT3DMS as a surrogate representation for heat transport in mid-enthalpy geothermal systems. A three-dimensional heat transport model was built, based on a relatively simple conceptual model. Results of simulation runs of a geothermal production scenario have positive implications for a working geothermal system at Eastgate. The Eastgate Geothermal Field has significant exploitation potential for combined heat and power purposes; it is anticipated that this site could support several tens of megawatts of heat production for direct use and many megawatts of electrical power using a binary power plant

Mobility and removal of nitrate in heterogeneous Eocene aquifers

n a study of the principles of lateral and vertical nitrate mobility and removal in unconfined heterogeneous fine sandy Eocene aquifers by long-term observation, no significant nitrate reduction could be observed over short distances (90 m) in the lateral flow direction. Decreasing oxygen contents and limited flow velocity in a downward direction (caused by clay lenses and layers) resulted in a more efficient but incomplete nitrate removal with increasing depth. Chemo-organotrophic and chemo-lithotrophic denitrification coexist in the aquifers. Recharge input of organic matter and the reactivity of sedimentary organic sources, as well as the amounts of pyrite and Fe2+-bearing minerals, control which microbiologically catalysed process finally occurs. Sharp boundaries between different redox zones do not exist due to locally changing availability of organic and inorganic electron donors used for nitrate reduction. Furthermore, preferential flow paths result in a wide spread occurrence of low concentrations of nitrate below the main denitrification zone

Impact of recent climate variability on an aquifer system in north Belgium

The last decade it has been realised that climate and global change can (and very likely will affect) groundwater reserves. This will have impact on both the economical (groundwater exploitation) and ecological (ecosystems) aspects of aquifer systems. A key issue is the intrinsic variability in hydrodynamics due to natural fluctuations of meteorological conditions. In this paper the results of a modelling study are presented that reconstructs the hydrodynamic evolution of the important Neogene aquifer system in North-East Belgium, between 1833 and 2005. Boundary conditions are defined on a monthly basis. The results show that besides the yearly seasonal fluctuations also multi-year to decadal variations occur. These are especially important in the topographic higher regions like on the Campine plateau and the Campine cuesta, and can severely affect local groundwater dependent ecosystems as found in some protected areas and nature reserves. The same cyclicity is recognized in the intensity and extent of seepage of deep Neogene groundwater in the Nete basin

Bringing the OpenMI to LIFE Progress Report No. 4 - 31st March 2008 – 30th September 2008

The Water Framework Directive demands an integrated approach to water management. This requires the ability to predict how catchment processes will behave and interact in response to the activities of water managers and others. In most contexts, it is not feasible to build a single predictive model that adequately represents all the processes; therefore a means of linking models of individual processes is required. This is met by the FP5 HarmonIT project’s Open Modelling Interface and Environment (the OpenMI). The purpose of this project is to transform the OpenMI from a research output to a sustainable operational Standard. It will build the capacity to use the OpenMI and will demonstrate it under operational conditions. It will also develop, test and demonstrate the future support organisation for the OpenMI. Finally, information about the OpenMI will be disseminated to users

Numerical simulation of radionuclides migration in the far field of a geological repository

Safety conditions associated with geological repositories must be guaranteed also in the case of radionuclides migrating from the near field to the far field of a geological repository and to the external environment. For this reason the migration process of radionuclides and the factors affecting the process patterns have a crucial importance. In the present article, in order to simulate the migration process of radionuclides in the far-field of a geological repository, the groundwater simulation code PMWIN (Processing Modflow) is used, following a methodology applied by the same authors in a previous work. The present case study refers to a non-uniform groundwater flow field and shows the influence of two important parameters, the distribution coefficient and the hydraulic gradient. The results are compared with the ones previously obtained for the uniform flow case and the different scenarios are discussed