A Geostatistical Investigation into Changing Permeability of Sandstones During Weathering Simulations

The ability to predict the behavior of masonry materials is crucial to conserve building stone. Natural stone, such as sandstone, is not immune from the processes of weathering in the built environment and suffers from decay by granular disintegration, contour scaling, and multiple flaking. Spatial variation of rock properties is a major contributing factor to inconsistent responses to weathering. This has implications for moisture movement and salt input and output and storage, and results in unpredictability in the decay dynamics of masonry materials. This article explores the use of variography and kriging to investigate the spatial interactions between the trigger factors of stone decay, in particular, permeability and its effect on salt penetration. Sandstone blocks were used to represent fresh building stones from a weathering perspective and gave baseline characteristics for the interpretation of subsequent deterioration and decay pathways. Simulated weathering trials involved preloading a sandstone block with salt and subjecting a separate block to 20 cycles of a weathering trial designed to simulate a temperate weathering regime. Geostatistical analysis indicated differences in the spatial variation of permeability of the fresh rock and that subjected to the weathering regimes. Spatial prediction and visualization showed differences in the spatial continuity of permeability in a horizontal and vertical direction through the preloaded block after salt weathering. Continual wetting with salt and alternate heating increased permeability in a vertical direction, enabling the ingress and movement of salt and moisture more effectively through the stone.

Dynamical instability in surface permeability characteristics of building sandstones in response to salt accumulation over time

This paper explores how the surface permeability of sandstone blocks changes over time in response to repeated salt weathering cycles. Surface permeability controls the amount of moisture and dissolved salt that can penetrate in and facilitate decay. Connected pores permit the movement of moisture (and hence soluble salts) into the stone interior, and where areas are more or less permeable soluble salts may migrate along preferred pathways at differential rates. Previous research has shown that salts can accumulate in the near-surface zone and lead to partial pore blocking which influences subsequent moisture ingress and causes rapid salt accumulation in the near-surface zone.

Two parallel salt weathering simulations were carried out on blocks of Peakmoor Sandstone of different volumes. Blocks were removed from simulations after 2, 5, 10, 20 and 60 cycles. Permeability measurements were taken for these blocks at a resolution of 20 mm, providing a grid of 100 permeability values for each surface. The geostatistical technique of ordinary kriging was applied to the data to produce a smoothed interpolation of permeability for these surfaces, and hence improve understanding of the evolution of permeability over time in response to repeated salt weathering cycles.

Results illustrate the different responses of the sandstone blocks of different volumes to repeated salt weathering cycles. In both cases, after an initial subtle decline in the permeability (reflecting pore blocking), the permeability starts to increase — reflected in a rise in mean, maximum and minimum values. However, between 10 and 20 cycles, there is a jump in the mean and range permeability of the group A block surfaces coinciding with the onset of meaningful debris release. After 60 cycles, the range of permeability in the group A block surface had increased markedly, suggesting the development of a secondary permeability. The concept of dynamic instability and divergent behaviour is applied at the scale of a single block surface, with initial small-scale differences across a surface having larger scale consequences as weathering progresses.

After cycle 10, group B blocks show a much smaller increase in mean permeability, and the range stays relatively steady — this may be explained by the capillary conditions set up by the smaller volume of the stone, allowing salts to migrate to the ‘back’ of the blocks and effectively relieving stress at the ‘front’ face.

A commentary on climate change, stone decay dynamics and the ‘greening’ of natural stone buildings: new perspectives on ‘deep wetting’

Environmental controls on stone decay processes are rapidly changing as a result of changing climate. UKCP09 projections for the 2020s (2010–2039) indicate that over much of the UK seasonality of precipitation will increase. Summer dryness and winter wetness are both set to increase, the latter linked to projected precipitation increases in autumn and spring months. If so, this could increase the time that stone structures remain wet and possibly the depth of moisture penetration, and it appears that building stone in Northern Ireland has already responded through an increased incidence of algal ‘greening’.This paper highlights the need for understanding the effects of climate change through a series of studies of largely sandstone structures. Current and projected climatic trends are therefore considered to have aesthetic, physical and chemical implications that are not currently built into our models of sandstone decay, especially with respect to the role played by deep-seated wetness on sandstone deterioration and decay progression and the feedbacks associated with, for example surface algal growth. In particular,it is proposed that algal biofilms will aid moisture retention and further facilitate moisture and dissolved salt penetration to depth. Thus, whilst the outer surface of stone may continue to experience frequent wetting and drying associated with individual precipitation events, the latter is less likely to be complete, and the interiors of building blocks may only experience wetting/drying in response to seasonal cycling. A possible consequence of deeper salt penetration could be a delay in the onset of surface deterioration,but more rapid and effective retreat once it commences as decay mechanisms ‘tap into a reservoir of deep salt’.

Bioaccessibility of trace elements in soils in Northern Ireland

Assessment of elevated concentrations of potentially toxic elements (PTE) in soils and the association with specific soil parent material have been the focus of research for a number of years. Risk-based assessment of potential exposure scenarios to identified elevated PTE concentrations has led to the derivation of site- and contaminant-specific soil guideline values (SGVs), which represent generic assessment criteria (GACs) to identify exceeded levels that may reflect an unacceptable risk to human health. A better understanding of the ‘bioavailable’ or ‘bioaccessible’ contaminant concentrations offers an opportunity to better refine contaminant exposure assessments. Utilizing a comprehensive soil geochemical dataset for Northern Ireland provided by the Tellus Survey (GSNI) in conjunction with supplementary bioaccessibility testing of selected soil samples following the Unified BARGE Method, this paper uses exploratory data analysis and geostatistical analysis to investigate the spatial variability of pseudo-total and bioaccessible concentrations of As, Cd, Co, Cr. Cu, Ni, Pb, U, V and Zn. The paper investigates variations in individual element concentrations as well as cross-element correlations and observed lithological/pedological associations. The analysis of PTE concentrations highlighted exceeded levels of GAC values for V and Cr and exceeded SGV/GAC values for Cd, Cu, Ni, Pb, and Zn. UBM testing showed that for some soil parent materials associated with elevated PTE concentrations e.g. the Antrim Lava Group with high Ni concentrations, the measured oral bioaccessible fraction was relatively low. For other soil parent materials with relatively moderate PTE concentrations, measured oral bioaccessible fraction was relatively high (e.g. the Gala Sandstone Group of the Southern Uplands-Down Longford Terrain). These findings have implications for regional human health risk assessments for specific PTEs.

How Porosity and Permeability Vary Spatially With Grain Size, Sorting, Cement Volume, and Mineral Dissolution In Fluvial Triassic Sandstones: The Value of Geostatistics and Local Regression

Although it is well known that sandstone porosity and permeability are controlled by a range of parameters such as grain size and sorting, amount, type, and location of diagenetic cements, extent and type of compaction, and the generation of intergranular and intragranular secondary porosity, it is less constrained how these controlling parameters link up in rock volumes (within and between beds) and how they spatially interact to determine porosity and permeability. To address these unknowns, this study examined Triassic fluvial sandstone outcrops from the UK using field logging, probe permeametry of 200 points, and sampling at 100 points on a gridded rock surface. These field observations were supplemented by laser particle-size analysis, thin-section point-count analysis of primary and diagenetic mineralogy, quantitiative XRD mineral analysis, and SEM/EDAX analysis of all 100 samples. These data were analyzed using global regression, variography, kriging, conditional simulation, and geographically weighted regression to examine the spatial relationships between porosity and permeability and their potential controls. The results of bivariate analysis (global regression) of the entire outcrop dataset indicate only a weak correlation between both permeability porosity and their diagenetic and depositional controls and provide very limited information on the role of primary textural structures such as grain size and sorting. Subdividing the dataset further by bedding unit revealed details of more local controls on porosity and permeability. An alternative geostatistical approach combined with a local modelling technique (geographically weighted regression; GWR) subsequently was used to examine the spatial variability of porosity and permeability and their controls. The use of GWR does not require prior knowledge of divisions between bedding units, but the results from GWR broadly concur with results of regression analysis by bedding unit and provide much greater clarity of how porosity and permeability and their controls vary laterally and vertically. The close relationship between depositional lithofacies in each bed, diagenesis, and permeability, porosity demonstrates that each influences the other, and in turn how understanding of reservoir properties is enhanced by integration of paleoenvironmental reconstruction, stratigraphy, mineralogy, and geostatistics.

The Influence of Aspect on the Biological Colonization of Stone in Northern Ireland.

The rate and type of biological colonization of stone is influenced by a wide array of environmental factors in addition to substrate characteristics. A series of experiments was designed to compare the rate and type of biological colonization of stone at varying locations over a 21-month time period. Exposure
trials were set up at nine different sites across Northern Ireland that covered a wide variety of environmental conditions. To determine aspect-related differences in colonization, blocks of Peakmoor sandstone and Portland limestone were placed on the north- and south-facing sides of purpose-designed exposure racks. Colorimetry and visual analysis were carried out on collected samples at increasing time intervals. Results showed significantly different rates of darkening and greening over time between north-facing and south-facing blocks, for both sandstone and limestone. This difference is likely to be representative of the fact that in Northern Ireland’s wet climate and northern-latitude position, the north face of a building will receive less direct sunlight. Therefore north-facing blocks, once wet, will remain damp for much longer than blocks on other façades. This slow-drying phenomenon is much more hospitable for biological colonization and continued growth than the hostile environment of rapid wetting and drying cycles experienced on the south face.

Changing climate, changing process: implications for salt transportation and weathering within building sandstones in the UK

Salt weathering is a crucial process that brings about a change in stone, from the scale of landscapes to stone outcrops and natural building stone facades. It is acknowledged that salt weathering is controlled by fluctuations in temperature and moisture, where repeated oscillations in these parameters can cause re-crystallisation, hydration/de-hydration of salts, bringing about stone surface loss in the form of, for example, granular disaggregation, scaling, and multiple flaking. However, this ‘traditional’ view of how salt weathering proceeds may need to be re-evaluated in the light of current and future climatic trends. Indeed, there is considerable scope for the investigation of consequences of climate change on geomorphological processes in general. Building on contemporary research on the ‘deep wetting’ of natural building stones, it is proposed that (as stone may be wetter for longer), ion diffusion may become a more prominent mechanism for the mixing of molecular constituents, and a shift in focus from physical damage to chemical change is suggested. Data from ion diffusion cell experiments are presented for three different sandstone types, demonstrating that salts may diffuse through porous stone relatively rapidly (in comparison to, for example, dense concrete). Pore water from stones undergoing diffusion experiments was extracted and analysed. Factors controlling ion diffusion
relating to ‘time of wetness’ within stones are discussed, (continued saturation, connectivity of pores, mineralogy, behaviour of salts, sedimentary structure), and potential changes in system dynamics as a result of climate change are addressed. System inputs may change in terms of increased moisture input, translating into a greater depth of wetting front. Salts are likely to be ‘stored’ differently in stones, with salt being in solution for longer periods (during prolonged winter wetness). This has myriad implications in terms of the movement of ions by diffusion and the potential for chemical change in the stone (especially in more mobile constituents), leading to a weakening of the stone matrix/grain boundary cementing. The ‘output’ may be mobilisation and precipitation of elements leading to, for example, uneven cementing in the stone. This reduced strength of the stone, or compromised ability of the stone to absorb stress, is likely to make crystallisation a more efficacious mechanism of decay when it does occur. Thus, a delay in the onset of crystallisation while stonework is wet does not preclude exaggerated or accelerated material loss when it finally happens.

Application of factor analysis and geostatistics to identify processes controlling the distribution and bioaccessibility of trace elements in Northern Ireland soils

Geologic and environmental factors acting over varying spatial scales can control
trace element distribution and mobility in soils. In turn, the mobility of an element in soil will affect its oral bioaccessibility. Geostatistics, kriging and principal component analysis (PCA) were used to explore factors and spatial ranges of influence over a suite of 8 element oxides, soil organic carbon (SOC), pH, and the trace elements nickel (Ni), vanadium (V) and zinc (Zn). Bioaccessibility testing was carried out previously using the Unified BARGE Method on a sub-set of 91 soil samples from the Northern Ireland Tellus1 soil archive. Initial spatial mapping of total Ni, V and Zn concentrations shows their distributions are correlated spatially with local geologic formations, and prior correlation analyses showed that statistically significant controls were exerted over trace element bioaccessibility by the 8 oxides, SOC and pH. PCA applied to the geochemistry parameters of the bioaccessibility sample set yielded three principal components accounting for 77% of cumulative variance in the data
set. Geostatistical analysis of oxide, trace element, SOC and pH distributions using 6862 sample locations also identified distinct spatial ranges of influence for these variables, concluded to arise from geologic forming processes, weathering processes, and localised soil chemistry factors. Kriging was used to conduct a spatial PCA of Ni, V and Zn distributions which identified two factors comprising the majority of distribution variance. This was spatially accounted for firstly by basalt rock types, with the second component associated with sandstone and limestone in the region. The results suggest trace element bioaccessibility and distribution is controlled by chemical and geologic processes which occur over variable spatial ranges of influence.

Numerical assessment of potential impacts of hydraulically fractured Bowland Shale on overlying aquifers

Natural gas extracted from hydraulically fractured shale formations potentially has a big impact on the global energy landscape. However, there are concerns of potential environmental impacts of hydraulic fracturing of the shale formations, particularly those related to water quality. To evaluate the potential impact of hydraulically fractured shale on overlying aquifers, we conduct realizations of numerical modeling simulations to assess fluid flow and chloride transport from a synthetic Bowland Shale over a period of 11,000 years. The synthetic fractured shale was represented by a three-dimensional discrete fracture model that was developed by using the data from a Bowland Shale gas exploration in Lancashire, UK. Chloride mass exchange between fractures and the rock matrix was fully accounted for in the model. The assessment was carried out to investigate fluid and chloride mass fluxes before, during, and after hydraulic fracturing of the Bowland Shale. Impacts of the upward fracture height and aperture, as well as hydraulic conductivity of the multilayered bedrock system, are also included this assessment. This modeling revealed that the hydraulically fractured Bowland Shale is unlikely to pose a risk to its overlying groundwater quality when the induced fracture aperture is ≤200 µm. With the fracture aperture ≥1000 µm, the upward chloride flux becomes very sensitive to the upward fracture height growth and hydraulic conductivity of the multilayered bedrock system. In the extremely unlikely event of the upward fracture growth directly connecting the shale formation to the overlying Sherwood Sandstone aquifer with the fracture aperture ≥1000 µm, the upward chloride mass flux could potentially pose risks to the overlying aquifer in 100 years. The model study also revealed that the upward mass flux is significantly intercepted by the horizontal mass flux within a high permeable layer between the Bowland Shale and its overlying aquifers, reducing further upward flux toward the overlying aquifers.

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.