Integrating geophysical data and upscaling techniques to support regional groundwater flow modelling: A practical example of the Lagan Valley, Northern Ireland

When studying heterogeneous aquifer systems, especially at regional scale, a degree of generalization is anticipated. This can be due to sparse sampling regimes, complex depositional environments or lack of accessibility to measure the subsurface. This can lead to an inaccurate conceptualization which can be detrimental when applied to groundwater flow models. It is important that numerical models are based on observed and accurate geological information and do not rely on the distribution of artificial aquifer properties. This can still be problematic as data will be modelled at a different scale to which it was collected. It is proposed here that integrating geophysics and upscaling techniques can assist in a more realistic and deterministic groundwater flow model. In this study, the sedimentary aquifer of the Lagan Valley in Northern Ireland is chosen due to intruding sub-vertical dolerite dykes. These dykes are of a lower permeability than the sandstone aquifer. The use of airborne magnetics allows the delineation of heterogeneities, confirmed by field analysis. Permeability measured at the field scale is then upscaled to different levels using a correlation with the geophysical data, creating equivalent parameters that can be directly imported into numerical groundwater flow models. These parameters include directional equivalent permeabilities and anisotropy. Several stages of upscaling are modelled in finite element. Initial modelling is providing promising results, especially at the intermediate scale, suggesting an accurate distribution of aquifer properties. This deterministic based methodology is being expanded to include stochastic methods of obtaining heterogeneity location based on airborne geophysical data. This is through the Direct Sample method of Multiple-Point Statistics (MPS). This method uses the magnetics as a training image to computationally determine a probabilistic occurrence of heterogeneity. There is also a need to apply the method to alternate geological contexts where the heterogeneity is of a higher permeability than the host rock.

Integrating Multiple Point Statistics with Aerial Geophysical Data to assist Groundwater Flow Models

The process of accounting for heterogeneity has made significant advances in statistical research, primarily in the framework of stochastic analysis and the development of multiple-point statistics (MPS). Among MPS techniques, the direct sampling (DS) method is tested to determine its ability to delineate heterogeneity from aerial magnetics data in a regional sandstone aquifer intruded by low-permeability volcanic dykes in Northern Ireland, UK. The use of two two-dimensional bivariate training images aids in creating spatial probability distributions of heterogeneities of hydrogeological interest, despite relatively ‘noisy’ magnetics data (i.e. including hydrogeologically irrelevant urban noise and regional geologic effects). These distributions are incorporated into a hierarchy system where previously published density function and upscaling methods are applied to derive regional distributions of equivalent hydraulic conductivity tensor K. Several K models, as determined by several stochastic realisations of MPS dyke locations, are computed within groundwater flow models and evaluated by comparing modelled heads with field observations. Results show a significant improvement in model calibration when compared to a simplistic homogeneous and isotropic aquifer model that does not account for the dyke occurrence evidenced by airborne magnetic data. The best model is obtained when normal and reverse polarity dykes are computed separately within MPS simulations and when a probability threshold of 0.7 is applied. The presented stochastic approach also provides improvement when compared to a previously published deterministic anisotropic model based on the unprocessed (i.e. noisy) airborne magnetics. This demonstrates the potential of coupling MPS to airborne geophysical data for regional groundwater modelling.

Diagenetic Alterations and Reservoir Quality Evolution of Lower Cretaceous Fluvial Sandstones: Nubian Formation, Sirt Basin, North-Central Libya

Lower Cretaceous meandering and braided fluvial sandstones of the Nubian Formation form some of the most important subsurface reservoir rocks in the Sirt Basin, north-central Libya. Mineralogical, petrographical and geochemical analyses of sandstone samples from well BB6-59, Sarir oilfield, indicate that the meandering fluvial sandstones are fine- to very fine-grained subarkosic arenites (av. Q91F5L4), and that braided fluvial sandstones are medium- to very coarse-grained quartz arenites (av. Q96F3L1). The reservoir qualities of these sandstones were modified during both eodiagenesis (ca. <70oC; <2 km) and mesodiagenesis (ca. >70oC; >2km). Reservoir quality evolution was controlled primarily by the dissolution and kaolinitization of feldspars, micas and mud intraclasts during eodiagenesis, and by the amount and thicknessof grain-coating clays, chemical compaction and quartz overgrowths during mesodiagenesis. However, dissolution and kaolinitization of feldspars, micas and mud intraclasts resulted in the creation of intercrystalline micro- and mouldic macro-porosity and permeability during eodiagenesis, which were more widespread in braided fluvial than in meandering fluvial sandstones. This was because of the greater depositional porosity and permeability in the braided fluvial sandstones which enhanced percolation of meteoric waters. The development of only limited quartz overgrowths in the braided fluvial sandstones, in which quartz grains are coated by thick illite layers, retained high porosity and permeability (12-23 % and 30- 600 mD). By contrast, meandering fluvial sandstones underwent porosity loss as a result of quartz overgrowth development on quartz grains which lack or have thin and incomplete grain-coating illite (2-15 % and 0-0.1mD). Further loss of porosity in the meandering fluvial sandstones occurred as a result of chemical compaction (pressuredissolution) induced by the occurrence of micas along grains contacts. Otherdiagenetic alterations, such as the growth of pyrite, siderite, dolomite/ankerite and albitization, had little impact on reservoir quality. The albitization of feldspars may have had minor positive influence on reservoir quality throughthe creation of intercrystalline micro-porosity between albite crystals.The results of this study show that diagenetic modifications of the braided and meandering fluvial sandstones in the Nubian Formation, and resulting changes in reservoir quality, are closely linked to depositional porosity and permeability. They are also linked to the thickness of grain-coating infiltrated clays, and to variations in detrital composition, particularly the amounts of mud intraclasts, feldspars and mica grains as well as climatic conditions.

Development of a biocidal treatment regime to inhibit biological growths on cultural heritage: BIODAM

Existing chemical treatments to prevent biological damage to monuments often involve considerable amounts of potentially dangerous and even poisonous biocides. The scientific approach described in this paper aims at a drastic reduction in the concentration of biocide applications by a polyphasic approach of biocides combined with cell permeabilisers, polysaccharide and pigment inhibitors and a photodynamic treatment. A variety of potential agents were screened to determine the most effective combination. Promising compounds were tested under laboratory conditions with cultures of rock deteriorating bacteria, algae, cyanobacteria and fungi. A subsequent field trial involved two sandstone types with natural biofilms. These were treated with multiple combinations of chemicals and exposed to three different climatic conditions. Although treatments proved successful in the laboratory, field trials were inconclusive and further testing will be required to determine the most effective treatment regime. While the most effective combination of chemicals and their application methodology is still being optimised, results to date indicate that this is a promising and effective treatment for the control of a wide variety of potentially damaging organisms colonising stone substrates

The influence of substrate texture on early biological colonization

Biological colonization of stone is a major concern in the preservation and presentation of cultural heritage. Colonization is typically associated with unpleasant soiling, and varying degrees of biodeterioration. A better understanding of why organisms grow where they do, will aid in
developing preventative, and treatment methods for biosoiling of cultural heritage. Sandstone exposure trials were set up at nine different locations across Northern Ireland to investigate the influences of local climate, local environmental,and micro-climatic factors on the early stages (up to 21 months) of biological colonization.
Results showed that, green and yellow soiling occurred on tooled stone surfaces, whereas darkening occurred preferentially on smooth surfaces. It is likely that different populations of organisms occur on these surfaces with green algae occurring on tooled surfaces due to slower drying rates (i.e. prolonged moisture retention), and cyanobacteria and fungi thriving on smooth surfaces due to their ability to withstand moisture fluctuation.

The inventorying of building stone performance-in-use to assess susceptibility to decay

The Natural Stone Database for Northern Ireland was constructed to address the paucity of information available to stone conservation practitioners within the region. Almost 2000 listed buildings and monuments were surveyed over three years to produce an interactive GIS database. This contained information on stone sources, together with details of stone condition and decay processes. This paper uses elements of this GIS to investigate stone decay patterns across Northern Ireland. The results demonstrate that as the level of stone decay increases, so does the proportion of buildings with sandstone as the primary stone type. It appears that a
higher open porosity level combined with Northern Ireland’s wetter climate and maritime location leads to rapid wetting and drying cycles within sandstones. This is coupled with the ingress and crystallisation of marine and other salts within stone pores leading to considerably
higher rates of decay than for any other stone type.

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.

Evaporative Moisture Loss from Heterogeneous Stone: Material- Environment Interactions During Drying

The complexities of evaporation from structurally and mineralogically heterogeneous sandstone (Locharbriggs Sandstone) are investigated through a laboratory-based experiment in which a variety of environmental conditions are simulated. Data reported demonstrate the significance of material-environment interactions on the spatial and temporal variability of evaporative dynamics. Evaporation from porous stone is determined by the interplay between environmental, material and solution properties, which govern the rate and mode by which water is transmitted to, and subsequently removed from, an evaporating surface. Initially evaporation is marked by high rates of moisture loss controlled by external atmospheric conditions; then, when a critical level of surface moisture content is reached, hydraulic continuity between the stone surface and subsurface is disrupted and the drying front recedes
beneath the surface, evaporation rates decrease and are controlled by the ability of the material to transport water vapour to the surface. Pore size distribution and connectivity, as well as other material properties, control the timing of each stage of evaporation and the nature of the transition.

These experimental data highlight the complexity of evaporation, demonstrating that different regions of the same stone can exhibit varying moisture dynamics during drying and that the rate and nature of evaporative loss differs under different environmental conditions. The results identify the importance of material-environment interactions during drying and that stone micro-environmental conditions cannot be inferred from ambient data alone.
These data have significance for understanding the spatial distribution of stone surface weathering-related morphologies in both the natural and built environments where mineralogical and/or structural heterogeneity creates differences in moisture flux and hence variable drying rates. Such differences may provide a clearer explanation for the initiation and subsequent development of complex weathering responses where areas of significant deterioration can be found alongside areas that exhibit little or no evidence surface breakdown.