Non-stationary variogram models for geostatistical sampling optimisation: An empirical investigation using elevation data

A problem with use of the geostatistical Kriging error for optimal sampling design is that the design does not adapt locally to the character of spatial variation. This is because a stationary variogram or covariance function is a parameter of the geostatistical model. The objective of this paper was to investigate the utility of non-stationary geostatistics for optimal sampling design. First, a contour data set of Wiltshire was split into 25 equal sub-regions and a local variogram was predicted for each. These variograms were fitted with models and the coefficients used in Kriging to select optimal sample spacings for each sub-region. Large differences existed between the designs for the whole region (based on the global variogram) and for the sub-regions (based on the local variograms). Second, a segmentation approach was used to divide a digital terrain model into separate segments. Segment-based variograms were predicted and fitted with models. Optimal sample spacings were then determined for the whole region and for the sub-regions. It was demonstrated that the global design was inadequate, grossly over-sampling some segments while under-sampling others.

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.

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.

Increasing the accuracy of Nitrogen Dioxide (NO2) pollution mapping using geographically weighted regression (GWR) and geostatistics.

Nitrogen Dioxide (NO2) is known to act as an environmental trigger for many respiratory illnesses. As a pollutant it is difficult to map accurately, as concentrations can vary greatly over small distances. In this study three geostatistical techniques were compared, producing maps of NO2 concentrations in the United Kingdom (UK). The primary data source for each technique was NO2 point data, generated from background automatic monitoring and background diffusion tubes, which are analysed by different laboratories on behalf of local councils and authorities in the UK. The techniques used were simple kriging (SK), ordinary kriging (OK) and simple kriging with a locally varying mean (SKlm). SK and OK make use of the primary variable only. SKlm differs in that it utilises additional data to inform prediction, and hence potentially reduces uncertainty. The secondary data source was Oxides of Nitrogen (NOx) derived from dispersion modelling outputs, at 1km x 1km resolution for the UK. These data were used to define the locally varying mean in SKlm, using two regression approaches: (i) global regression (GR) and (ii) geographically weighted regression (GWR). Based upon summary statistics and cross-validation prediction errors, SKlm using GWR derived local means produced the most accurate predictions. Therefore, using GWR to inform SKlm was beneficial in this study.

Spatial and temporal ecology of Scots pine ectomycorrhizas

Spatial analysis was used to explore the distribution of individual species in an ectomycorrhizal (ECM) fungal community to address: whether mycorrhizas of individual ECM fungal species were patchily distributed, and at what scale; and what the causes of this patchiness might be. Ectomycorrhizas were extracted from spatially explicit samples of the surface organic horizons of a pine plantation. The number of mycorrhizas of each ECM fungal species was recorded using morphotyping combined with internal transcribed spacer (ITS) sequencing. Semivariograms, kriging and cluster analyses were used to determine both the extent and scale of spatial autocorrelation in species abundances, potential interactions between species, and change over time. The mycorrhizas of some, but not all, ECM fungal species were patchily distributed and the size of patches differed between species. The relative abundance of individual ECM fungal species and the position of patches of ectomycorrhizas changed between years. Spatial and temporal analysis revealed a dynamic ECM fungal community with many interspecific interactions taking place, despite the homogeneity of the host community. The spatial pattern of mycorrhizas was influenced by the underlying distribution of fine roots, but local root density was in turn influenced by the presence of specific fungal species.

Patterns of Holocene relative sea level change in the North of Britain and Ireland

Temporal and spatial patterns of relative sea level (RSL) change in the North of Britain and Ireland during the Holocene are examined. Four episodes, each defined by marked changes in the RSL trend, are identified. Each episode is marked by a rise to a culminating shoreline followed by a fall. Episode HRSL-1 dates from the Younger Dryas to early in the Holocene; HRSL-2 to HRSL-4 occurred later in the Holocene. There is extensive evidence for each episode, and on this basis the spatial distribution of the altitude data for three culminating shorelines and a shoreline formed at the time of the Holocene Storegga Slide tsunami (ca 8110 ± 100 cal. BP) is analysed. Ordinary Kriging is used to determine the general pattern, following which Gaussian Trend Surface Analysis is employed. Recognising that empirical measurements of RSL change can be unevenly distributed spatially, a new approach is introduced which enables the developing pattern to be identified. The patterns for the most widely occurring shorelines were analysed and found to be similar and common centre and axis models were developed for all shorelines. The analyses described provide models of the spatial pattern of Holocene RSL change in the area between ca 8100 cal. BP and ca 1000 cal. BP based on 2262 high resolution shoreline altitude measurements. These models fit the data closely, no shoreline altitude measurement lying more than −1.70 m or +1.82 m from the predicted value. The models disclose a similar pattern to a recently published Glacial Isostatic Adjustment model for present RSL change across the area, indicating that the overall spatial pattern of RSL change may not have varied greatly during the last ca 8000 years.

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.

Spatial statistics to estimate peat thickness using airborne radiometric data

Soil carbon stores are a major component of the annual returns required by EU governments to the Intergovernmental Panel on Climate Change. Peat has a high proportion of soil carbon due to the relatively high carbon density of peat and organic-rich soils. For this reason it has become increasingly important to measure and model soil carbon stores and changes in peat stocks to facilitate the management of carbon changes over time. The approach investigated in this research evaluates the use of airborne geophysical (radiometric) data to estimate peat thickness using the attenuation of bedrock geology radioactivity by superficial peat cover. Remotely sensed radiometric data are validated with ground peat depth measurements combined with non-invasive geophysical surveys. Two field-based case studies exemplify and validate the results. Variography and kriging are used to predict peat thickness from point measurements of peat depth and airborne radiometric data and provide an estimate of uncertainty in the predictions. Cokriging, by assessing the degree of spatial correlation between recent remote sensed geophysical monitoring and previous peat depth models, is used to examine changes in peat stocks over time. The significance of the coregionalisation is that the spatial cross correlation between the remote and ground based data can be used to update the model of peat depth. The result is that by integrating remotely sensed data with ground geophysics, the need is reduced for extensive ground-based monitoring and invasive peat depth measurements. The overall goal is to provide robust estimates of peat thickness to improve estimates of carbon stocks. The implications from the research have a broader significance that promotes a reduction in the need for damaging onsite peat thickness measurement and an increase in the use of remote sensed data for carbon stock estimations.

Trace Element Abundance and Renal Disease

This research investigates the relationship between elevated trace elements in soils, stream sediments and stream water and the prevalence of Chronic Kidney Disease (CKD). The study uses a collaboration of datasets provided from the UK Renal Registry Report (UKRR) on patients with renal diseases requiring treatment including Renal Replacement Therapy (RRT), the soil geochemical dataset for Northern Ireland provided by the Tellus Survey, Geological Survey of Northern Ireland (GSNI) and the bioaccessibility of Potentially Toxic Elements (PTEs) from soil samples which were obtained from the Unified Barge Method (UBM). The relationship between these factors derives from the UKRR report which highlights incidence rates of renal impaired patients showing regional variation with cases of unknown aetiology. Studies suggest a potential cause of the large variation and uncertain aetiology is associated with underlying environmental factors such as the oral bioaccessibility of trace elements in the gastrointestinal tract.
As previous research indicates that long term exposure is related to environmental factors, Northern Ireland is ideally placed for this research as people traditionally live in the same location for long periods of time. Exploratory data analysis and multivariate analyses are used to examine the soil, stream sediments and stream water geochemistry data for a range of key elements including arsenic, lead, cadmium and mercury identified from a review of previous renal disease literature. The spatial prevalence of patients with long term CKD is analysed on an area basis. Further work includes cluster analysis to detect areas of low or high incidences of CKD that are significantly correlated in space, Geographical Weighted Regression (GWR) and Poisson kriging to examine locally varying relationship between elevated concentrations of PTEs and the prevalence of CKD.

Surrogate-based analysis and optimization for the design of heat sinks with jet impingement

Heat sinks are widely used for cooling electronic devices and systems. Their thermal performance is usually determined by the material, shape, and size of the heat sink. With the assistance of computational fluid dynamics (CFD) and surrogate-based optimization, heat sinks can be designed and optimized to achieve a high level of performance. In this paper, the design and optimization of a plate-fin-type heat sink cooled by impingement jet is presented. The flow and thermal fields are simulated using the CFD simulation; the thermal resistance of the heat sink is then estimated. A Kriging surrogate model is developed to approximate the objective function (thermal resistance) as a function of design variables. Surrogate-based optimization is implemented by adaptively adding infill points based on an integrated strategy of the minimum value, the maximum mean square error approach, and the expected improvement approaches. The results show the influence of design variables on the thermal resistance and give the optimal heat sink with lowest thermal resistance for given jet impingement conditions.