Quantitative discrimination between geological materials with variable density contrast by high resolution X-ray computed tomography: An example using amygdule size-distribution in ancient lava flows

Abstract

The bubble-size distribution in 2.7 billion year old lava flows can be used as a proof of concept illustrating a new set of techniques for measuring volumes of geological materials with variable density contrasts using high-resolution X-ray computed tomography. Such studies have been limited in the past to high-contrast situations such as vesicles devoid of secondary fill. We present a new dynamic thresholding method for computationally separating amygdules from their basaltic matrix in X-ray images that is based on a technique used in seismology. The technique is sensitive to the gradient of the gray-scale value, rather than an absolute threshold value often applied to an entire set of X-ray images. Additionally, we present statistical methods for extrapolating the volumetric measurement mean and standard deviation of amygdules in the measured samples to the entire population in the flow. To do so, we create additional amygdule sample sets from the original sample set in the process of ‘bootstrap’ resampling, and use the Central Limit Theorem to calculate the mean and standard deviation of the amygdule population from these sample sets. This suite of methods allows the extension of bubble-size distribution studies typically done on modern flows to the ancient rock record and potentially has many other uses in geosciences where quantitative discrimination between materials with a range of densities is required.