The use of aero-magnetics to enhance a numerical groundwater model of the Lagan Valley aquifer, Northern Ireland

Peer reviewedPublisher PD

Effect of volcanic dykes on coastal groundwater flow and saltwater intrusion : a field-scale multiphysics approach and parameter evaluation

Acknowledgments This research was primarily based on research grant‐aided by the Irish Department of Communications, Energy and Natural Resources under the National Geoscience Programme 2007–2013. It also benefited from complementary funding from the Scottish Alliance for Geoscience, Environment and Society (SAGES). We acknowledge the contribution in data acquisition of the MSc students in Environmental Engineering at Queen's University Belfast, the landowner for access to the inland fields and the Department of Geography, Archaeology and Paleoecology at QUB for provision of the tidal model of Belfast Lough. The data used are listed in the references, tables, and figures and are available from the corresponding author upon demand. We acknowledge the constructive comments by the Associate Editor and three reviewers, which helped in improving the final manuscript.Peer reviewedPublisher PD

Effect of sedimentary heterogeneities in the sealing formation on predictive analysis of geological CO₂ storage

Numerical models of geologic carbon sequestration (GCS) in saline aquifers use multiphase fluid flow-characteristic curves (relative permeability and capillary pressure) to represent the interactions of the non-wetting CO2 and the wetting brine. Relative permeability data for many sedimentary formations is very scarce, resulting in the utilisation of mathematical correlations to generate the fluid flow characteristics in these formations. The flow models are essential for the prediction of CO2 storage capacity and trapping mechanisms in the geological media. The observation of pressure dissipation across the storage and sealing formations is relevant for storage capacity and geomechanical analysis during CO2 injection. This paper evaluates the relevance of representing relative permeability variations in the sealing formation when modelling geological CO2 sequestration processes. Here we concentrate on gradational changes in the lower part of the caprock, particularly how they affect pressure evolution within the entire sealing formation when duly represented by relative permeability functions. The results demonstrate the importance of accounting for pore size variations in the mathematical model adopted to generate the characteristic curves for GCS analysis. Gradational changes at the base of the caprock influence the magnitude of pressure that propagates vertically into the caprock from the aquifer, especially at the critical zone (i.e. the region overlying the CO2 plume accumulating at the reservoir-seal interface). A higher degree of overpressure and CO2 storage capacity was observed at the base of caprocks that showed gradation. These results illustrate the need to obtain reliable relative permeability functions for GCS, beyond just permeability and porosity data. The study provides a formative principle for geomechanical simulations that study the possibility of pressure-induced caprock failure during CO2 sequestration

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

Elimination of the reaction rate 'scale effect': application of the Lagrangian reactive particle-tracking method to simulate mixing-limited, field-scale biodegradation at the Schoolcraft (MI, USA) site

This is the peer reviewed version of the following article: [Ding, D., Benson, D. A., Fernàndez‐Garcia, D., Henri, C. V., Hyndman, D. W., Phanikumar, M. S., & Bolster, D. (2017). Elimination of the reaction rate “scale effect”: Application of the Lagrangian reactive particle‐tracking method to simulate mixing‐limited, field‐scale biodegradation at the Schoolcraft (MI, USA) site. Water Resources Research, 53, 10,411–10,432. https://doi.org/10.1002/2017WR021103], which has been published in final form at https://doi.org/10.1002/2017WR021103. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving.Measured (or empirically fitted) reaction rates at groundwater remediation sites are typically much lower than those found in the same material at the batch or laboratory scale. The reduced rates are commonly attributed to poorer mixing at the larger scales. A variety of methods have been proposed to account for this scaling effect in reactive transport. In this study, we use the Lagrangian particle-tracking and reaction (PTR) method to simulate a field bioremediation experiment at the Schoolcraft, MI site. A denitrifying bacterium, Pseudomonas Stutzeri strain KC (KC), was injected to the aquifer, along with sufficient substrate, to degrade the contaminant, carbon tetrachloride (CT), under anaerobic conditions. The PTR method simulates chemical reactions through probabilistic rules of particle collisions, interactions, and transformations to address the scale effect (lower apparent reaction rates for each level of upscaling, from batch to column to field scale). In contrast to a prior Eulerian reaction model, the PTR method is able to match the field-scale experiment using the rate coefficients obtained from batch experiments.Peer ReviewedPostprint (author's final draft

CO2 capture and storage (CCS) selected papers – Technical annotated bibliography

Bureau of Economic Geolog

Numerical modelling in support of a conceptual model for groundwater flow and geochemical evolution in the southern Outaouais Region, Quebec, Canada

A two-dimensional vertical-section numerical model for groundwater flow and transport using age, tritium and chloride was used to help validate a conceptual model of geochemical evolution within a representative regional-scale hydrogeological system in the Outaouais Region, Quebec, Canada. The flow system includes up to 30 m of Quaternary sediments and marine clays overlying fractured silicate rock of the Canadian Shield. Calibration of the regional flow model using observed piezometric levels and tritium concentrations showed that shallow groundwater flow is dominated by local flow systems limited to 30–40 m depth, 1–5 km long, and with groundwater residence times of 10–50 years. Intermediate systems, on the order of 5–15 km long, are less extensive than initially thought and are characterised by maximum depths of about 100 m and residence times of 200–6000 years. A model-calibrated hydraulic conductivity of 8 × 10−5 m.s−1 was required in the upper 50 m of the fractured bedrock. The active flow zone was inferred to extend to depths of about 100–150 m, with any deeper regional flow essentially negligible. Differences between tritium-based ages and simulated mean residence times were attributed to mixing of groundwater in open boreholes. Concentrations of 4He could be explained by diffusive transport from deeper and older groundwater, exacerbated by sampling. With new insight from the numerical modelling, the conceptual flow model has been updated to now include only a weak component of regional flow combined with significant local- and intermediate-scale flow systems connected to the upper fractured bedrock. The simulated flow system is also consistent with the geochemical evolution of the region, which is dominated by young Ca-HCO3-type waters in the unconfined aquifer and by older Cl− signatures from the remnant Champlain Sea seawater

Simulation of the spatio-temporal extent of groundwater flooding using statistical methods of hydrograph classification and lumped parameter models

This article presents the development of a relatively low cost and rapidly applicable methodology to simulate the spatio-temporal occurrence of groundwater flooding in chalk catchments. In winter 2000/2001 extreme rainfall resulted in anomalously high groundwater levels and groundwater flooding in many chalk catchments of northern Europe and the southern United Kingdom. Groundwater flooding was extensive and prolonged, occurring in areas where it had not been recently observed and, in places, lasting for 6 months. In many of these catchments, the prediction of groundwater flooding is hindered by the lack of an appropriate tool, such as a distributed groundwater model, or the inability of models to simulate extremes adequately. A set of groundwater hydrographs is simulated using a simple lumped parameter groundwater model. The number of models required is minimized through the classification and grouping of groundwater level time-series using principal component analysis and cluster analysis. One representative hydrograph is modelled then transposed to other observed hydrographs in the same group by the process of quantile mapping. Time-variant groundwater level surfaces, generated using the discrete set of modelled hydrographs and river elevation data, are overlain on a digital terrain model to predict the spatial extent of groundwater flooding. The methodology is applied to the Pang and Lambourn catchments in southern England for which monthly groundwater level time-series exist for 52 observation boreholes covering the period 1975–2004. The results are validated against observed groundwater flood extent data obtained from aerial surveys and field mapping. The method is shown to simulate the spatial and temporal occurrence of flooding during the 2000/2001 flood event accurately