Improved estimates of percolation and anisotropic permeability from 3-D X-ray microtomography using stochastic analyses and visualization

X-ray microtomography (micro-CT) with micron resolution enables new ways of characterizing microstructures and opens pathways for forward calculations of multiscale rock properties. A quantitative characterization of the microstructure is the first step in this challenge. We developed a new approach to extract scale-dependent characteristics of porosity, percolation, and anisotropic permeability from 3-D microstructural models of rocks. The Hoshen-Kopelman algorithm of percolation theory is employed for a standard percolation analysis. The anisotropy of permeability is calculated by means of the star volume distribution approach. The local porosity distribution and local percolation probability are obtained by using the local porosity theory. Additionally, the local anisotropy distribution is defined and analyzed through two empirical probability density functions, the isotropy index and the elongation index. For such a high-resolution data set, the typical data sizes of the CT images are on the order of gigabytes to tens of gigabytes; thus an extremely large number of calculations are required. To resolve this large memory problem parallelization in OpenMP was used to optimally harness the shared memory infrastructure on cache coherent Non-Uniform Memory Access architecture machines such as the iVEC SGI Altix 3700Bx2 Supercomputer. We see adequate visualization of the results as an important element in this first pioneering study.

Measuring root traits in barley (Hordeum vulgare ssp vulgare and ssp spontaneum) seedlings using gel chambers, soil sacs and X-ray microtomography

Root characteristics of seedlings of five different barley genotypes were analysed in 2D using gel chambers, and in 3D using soil sacs that were destructively harvested and pots of soil that were assessed non-invasively using X-ray microtomography. After 5 days, Chime produced the greatest number of root axes (similar to 6) and Mehola significantly less (similar to 4) in all growing methods. Total root length was longest in GSH01915 and shortest in Mehola for all methods, but both total length and average root diameter were significantly larger for plants grown in gel chambers than those grown in soil. The ranking of particular growth traits (root number, root angular spread) of plants grown in gel plates, soil sacs and X-ray pots was similar, but plants grown in the gel chambers had a different order of ranking for root length to the soil-grown plants. Analysis of angles in soil-grown plants showed that Tadmore had the most even spread of individual roots and Chime had a propensity for non-uniform distribution and root clumping. The roots of Mehola were less well spread than the barley cultivars supporting the suggestion that wild and landrace barleys tend to have a narrower angular spread than modern cultivars. The three dimensional analysis of root systems carried out in this study provides insights into the limitations of screening methods for root traits and useful data for modelling root architecture.

Tracking larval insect movement within soil using high resolution X-ray microtomography

1. In contrast to above-ground insects, comparatively little is known about the behaviour of subterranean insects, due largely to the difficulty of studying them in situ. 2. The movement of newly hatched (neonate) clover root weevil (Sitona lepidus L. Coleoptera: Curculinidae) larvae was studied non-invasively using recently developed high resolution X-ray microtomography. 3. The movement and final position of S. lepidus larvae in the soil was reliably established using X-ray microtomography, when compared with larval positions that were determined by destructively sectioning the soil column. 4. Newly hatched S. lepidus larvae were seen to attack the root rhizobial nodules of their host plant, white clover (Trifolium repens L.). Sitona lepidus larvae travelled between 9 and 27 mm in 9 h at a mean speed of 1.8 mm h(-1). 5. Sitona lepidus larvae did not move through the soil in a linear manner, but changed trajectory in both the lateral and vertical planes.

Dose optimization approach to fast X-ray microtomography of the lung alveoli

A basic prerequisite for in vivo X-ray imaging of the lung is the exact determination of radiation dose. Achieving resolutions of the order of micrometres may become particularly challenging owing to increased dose, which in the worst case can be lethal for the imaged animal model. A framework for linking image quality to radiation dose in order to optimize experimental parameters with respect to dose reduction is presented. The approach may find application for current and future in vivo studies to facilitate proper experiment planning and radiation risk assessment on the one hand and exploit imaging capabilities on the other.

X-ray microtomography analysis of the damage micromechanisms in 3D woven composites under low-velocity impact

3D woven composites reinforced with either S2 glass, carbon or a hybrid combination of both and containing either polyethylene or carbon z-yarns were tested under low-velocity impact. Different impact energies (in the range of 21–316 J) were used and the mechanical response (in terms of the impact strength and energy dissipated) was compared with that measured in high-performance, albeit standard, 2D laminates. It was found that the impact strength in both 2D and 3D materials was mainly dependent on the in-plane fiber fracture. Conversely, the energy absorption capability was primarily influenced by the presence of z-yarns, having the 3D composites dissipated over twice the energy than the 2D laminates, irrespective of their individual characteristics (fiber type, compaction degree, porosity, etc.). X-ray microtomography revealed that this improvement was due to the z-yarns, which delayed delamination and maintained the structural integrity of the laminate, promoting energy dissipation by tow splitting, intensive fiber breakage under the tup and formation of a plug by out-of-plane shear.

X-ray microtomography for studying 3D-textures of speleothems developed inside historic walls

En este artículo se muestra una aplicación de la microtomografía computerizada de Rayos X (microCT-RX) como técnica no destructiva útil para la caracterización del interior de estructuras sin necesidad de perder la muestra. Gracias a la sensibilidad de la técnica ha sido posible distinguir diferentes tipos de crecimiento espeleotémico dentro de una estalactita localizada en las bóvedas interiores de la muralla histórica de la isla de Nueva Tabarca.

Application of X-ray microtomography to numerical simulations of agglomerate breakage by distinct element method

Computer-aided tomography has been used for many years to provide significant information about the internal properties of an object, particularly in the medical fraternity. By reconstructing one-dimensional (ID) X-ray images, 2D cross-sections and 3D renders can provide a wealth of information about an object's internal structure. An extension of the methodology is reported here to enable the characterization of a model agglomerate structure. It is demonstrated that methods based on X-ray microtomography offer considerable potential in the validation and utilization of distinct element method simulations also examined.

X-ray microtomography of 410 million-year-old optic capsules from placoderm fishes

The structure of two small ossified optic capsules from mid-Palaeozoic placoderm fishes has been revealed in fine detail, by the use of Xray microtomography analysis and 3D visualisation software. These two specimens are 410 million-year-old; they were collected from an Early Devonian (Lochkovian) limestone in central New South Wales, and are the oldest known optic capsules from jawed fishes. The capsules show attachment areas for seven extrinsic eye muscles, rather than the six until recently deemed universal for gnathostomes. The analysis also revealed structures within the ossified cartilage which covered the medial surface of the eyeball, including nerve tracts, vascular canals, and possibly a choroid rete mirabile. (c) 2005 Elsevier Ltd. All rights reserved.

Assessment of physical and hydrological properties of biological soil crusts using x-ray microtomography and modelling

Biological soil crusts (BSCs) are formed by aggregates of soil particles and communities of microbial organisms and are common in all drylands. The role of BSCs on infiltration remains uncertain due to the lack of data on their role in affecting soil physical properties such as porosity and structure. Quantitative assessment of these properties is primarily hindered by the fragile nature of the crusts. Here we show how the use of a combination of non-destructive imaging X-ray microtomography (XMT) and Lattice Boltzmann method (LBM) enables quantification of key soil physical parameters and the modeling of water flow through BSCs samples from Kalahari Sands, Botswana. We quantify porosity and flow changes as a result of mechanical disturbance of such a fragile cyanobacteria-dominated crust. Results show significant variations in porosity between different types of crusts and how they affect the flow and that disturbance of a cyanobacteria-dominated crust results in the breakdown of larger pore spaces and reduces flow rates through the surface layer. We conclude that the XMT–LBM approach is well suited for study of fragile surface crust samples where physical and hydraulic properties cannot be easily quantified using conventional methods.

Macropores generated during shrinkage in two paddy soils using X-ray micro-computed tomography

Soil shrinkage curve represents a decrease of total porosity or an increase of bulk density with water loss. However, our knowledge of the dynamics of pores and their geometry during soil shrinkage is scarce, partially due to lack of reliable methods for determining soil pores in relation to change in soil water. This study aimed to investigate the dynamics of macropores (>30 mu m) of paddy soils during shrinkage. Two, paddy soils, which were sampled from one paddy field cultivated for 20 years (YPF) and the other one for over 100 years (OPF), represented difference in crack geometry in the field. Macropore parameters (volume, connectivity, and orientation of pores) and soil shrinkage parameters were determined on the same undisturbed soil cores by X-ray microtomography and shrinkage curve, respectively. Macroporosity was on average four times larger in the YPF than in the OPF whereas the shrinkage capacity was lower in the YPF as compared to the OPF (0.09 vs. 0.15 COLE). Soil shrinkage increased the volume of pores by 3.7% in the YPF and by 1.6% in the OPF as well as their connectivity. The formation of macropores occurred mostly in the proportional shrinkage phase. As a result, the slope of the proportional shrinkage phase was smaller in the YPF (0.65) than in the OPF (0.89). New macropores were cracks and extended pre-existing pores in the range of 225-1215 pm size without any preferential orientation. This work provides image evidences that in paddy soils with high shrinkage capacity more macropores are generated in the soil presenting a smaller proportional shrinkage slope. (C) 2015 Elsevier B.V. All rights reserved.

Estimating root–soil contact from 3D X-ray microtomographs

Adequate contact with the soil is essential for water and nutrient adsorption by plant roots, but the determination of root–soil contact is a challenging task because it is difficult to visualize roots in situ and quantify their interactions with the soil at the scale of micrometres. A method to determine root–soil contact using X-ray microtomography was developed. Contact areas were determined from 3D volumetric images using segmentation and iso-surface determination tools. The accuracy of the method was tested with physical model systems of contact between two objects (phantoms). Volumes, surface areas and contact areas calculated from the measured phantoms were compared with those estimated from image analysis. The volume was accurate to within 0.3%, the surface area to within 2–4%, and the contact area to within 2.5%. Maize and lupin roots were grown in soil (<2 mm) and vermiculite at matric potentials of −0.03 and −1.6 MPa and in aggregate fractions of 4–2, 2–1, 1–0.5 and < 0.5 mm at a matric potential of −0.03 MPa. The contact of the roots with their growth medium was determined from 3D volumetric images. Macroporosity (>70 µm) of the soil sieved to different aggregate fractions was calculated from binarized data. Root-soil contact was greater in soil than in vermiculite and increased with decreasing aggregate or particle size. The differences in root–soil contact could not be explained solely by the decrease in porosity with decreasing aggregate size but may also result from changes in particle and aggregate packing around the root.

Imaging and modeling of pore scale processes in porous media using X-ray computed tomography and lattice Boltzmann solver

In der Erdöl– und Gasindustrie sind bildgebende Verfahren und Simulationen auf der Porenskala im Begriff Routineanwendungen zu werden. Ihr weiteres Potential lässt sich im Umweltbereich anwenden, wie z.B. für den Transport und Verbleib von Schadstoffen im Untergrund, die Speicherung von Kohlendioxid und dem natürlichen Abbau von Schadstoffen in Böden. Mit der Röntgen-Computertomografie (XCT) steht ein zerstörungsfreies 3D bildgebendes Verfahren zur Verfügung, das auch häufig für die Untersuchung der internen Struktur geologischer Proben herangezogen wird. Das erste Ziel dieser Dissertation war die Implementierung einer Bildverarbeitungstechnik, die die Strahlenaufhärtung der Röntgen-Computertomografie beseitigt und den Segmentierungsprozess dessen Daten vereinfacht. Das zweite Ziel dieser Arbeit untersuchte die kombinierten Effekte von Porenraumcharakteristika, Porentortuosität, sowie die Strömungssimulation und Transportmodellierung in Porenräumen mit der Gitter-Boltzmann-Methode. In einer zylindrischen geologischen Probe war die Position jeder Phase auf Grundlage der Beobachtung durch das Vorhandensein der Strahlenaufhärtung in den rekonstruierten Bildern, das eine radiale Funktion vom Probenrand zum Zentrum darstellt, extrahierbar und die unterschiedlichen Phasen ließen sich automatisch segmentieren. Weiterhin wurden Strahlungsaufhärtungeffekte von beliebig geformten Objekten durch einen Oberflächenanpassungsalgorithmus korrigiert. Die Methode der „least square support vector machine” (LSSVM) ist durch einen modularen Aufbau charakterisiert und ist sehr gut für die Erkennung und Klassifizierung von Mustern geeignet. Aus diesem Grund wurde die Methode der LSSVM als pixelbasierte Klassifikationsmethode implementiert. Dieser Algorithmus ist in der Lage komplexe geologische Proben korrekt zu klassifizieren, benötigt für den Fall aber längere Rechenzeiten, so dass mehrdimensionale Trainingsdatensätze verwendet werden müssen. Die Dynamik von den unmischbaren Phasen Luft und Wasser wird durch eine Kombination von Porenmorphologie und Gitter Boltzmann Methode für Drainage und Imbibition Prozessen in 3D Datensätzen von Böden, die durch synchrotron-basierte XCT gewonnen wurden, untersucht. Obwohl die Porenmorphologie eine einfache Methode ist Kugeln in den verfügbaren Porenraum einzupassen, kann sie dennoch die komplexe kapillare Hysterese als eine Funktion der Wassersättigung erklären. Eine Hysterese ist für den Kapillardruck und die hydraulische Leitfähigkeit beobachtet worden, welche durch die hauptsächlich verbundenen Porennetzwerke und der verfügbaren Porenraumgrößenverteilung verursacht sind. Die hydraulische Konduktivität ist eine Funktion des Wassersättigungslevels und wird mit einer makroskopischen Berechnung empirischer Modelle verglichen. Die Daten stimmen vor allem für hohe Wassersättigungen gut überein. Um die Gegenwart von Krankheitserregern im Grundwasser und Abwässern vorhersagen zu können, wurde in einem Bodenaggregat der Einfluss von Korngröße, Porengeometrie und Fluidflussgeschwindigkeit z.B. mit dem Mikroorganismus Escherichia coli studiert. Die asymmetrischen und langschweifigen Durchbruchskurven, besonders bei höheren Wassersättigungen, wurden durch dispersiven Transport aufgrund des verbundenen Porennetzwerks und durch die Heterogenität des Strömungsfeldes verursacht. Es wurde beobachtet, dass die biokolloidale Verweilzeit eine Funktion des Druckgradienten als auch der Kolloidgröße ist. Unsere Modellierungsergebnisse stimmen sehr gut mit den bereits veröffentlichten Daten überein.

Ethnic differences in the BMI-%BF relationships between young Japanese and Australian-Caucasian males living in Australia using dual-energy X-ray absorptiometry

The Body Mass Index (BMI) has been used worldwide as an indicator of fatness. However, the universal cut-off points by the World Health Organisation (WHO) classification may not be appropriate for every ethnic group when consider the relationship with their actual total body fatness(%BF). The application of population-specific classifications to assess BMI may be more relevant to public health. Ethnic differences in the BMI%BF relationship between 45 Japanese and 42 Australian-Caucasian males were assessed using whole body dual-energy X-ray absorptiometry (DXA) scan and anthropometry using a standard protocol. Japanese males had significantly (p<0.05) greater %BF at given BMI values than Australian males. When this is taken into account the newly proposed Asia-Pacific BMI classification of BMI 23 as overweight and 25 as obese may better assess the level of obesity that is associated increased health risks for this population. To clarify the current findings, further studies that compare the relationships across other Japanese populations are recommended.