Explore the vast range of titles available.

Groundwater contaminant transport problems remain challenging with respect to their computing requirements. This often limits the exploration of the conceptual uncertainty that is mainly related to large-scale geological features – such as faults, fractures, and stratigraphic variations – and limited characterization. Here, to facilitate geological conceptual uncertainty exploration, we develop further the use of graph representation for geological models to approximate groundwater flow and transport. We consider a faulted multi-heterogeneous-layer medium to test our approach. The existing rank correlation between the shortest path distribution from a contaminant source to the model domain outlet and the cumulative mass distribution at the outlet enables us to perform scenario selection. The scenario selection approach relies on a metric combining the Jaccard dissimilarity and the Wasserstein distance to compare binary images. Among a set combining eight alternative scenarios, where three faults can act as either a flow barrier or a preferential path, we show that the use of graph approximations allows us to retain or reject scenarios with confidence, as well as to estimate the individual probability of a fault to act as a barrier or a path. This methodology framework opens up possibilities to explore more thoroughly conceptual geological uncertainty for processes affected by flow and transport.

The increasing catalogue of high-quality ice-penetrating radar data provides a unique insight in the internal layering architecture of the Greenland ice sheet. The stratigraphy, an indicator of past deformation, highlights irregularities in ice flow and reveals large perturbations without obvious links to bedrock shape. In this work, to establish a new conceptual model for the formation process, we analysed the radar data at the onset of the Petermann Glacier, North Greenland, and created a three-dimensional model of several distinct stratigraphic layers. We demonstrate that the dominant structures are cylindrical folds sub-parallel to the ice flow. By numerical modelling, we show that these folds can be formed by lateral compression of mechanically anisotropic ice, while a general viscosity contrast between layers would not lead to folding for the same boundary conditions. We conclude that the folds primarily form by converging flow as the mechanically anisotropic ice is channelled towards the glacier. A range of mechanisms has been proposed for large-scale folding in polar ice sheets. Here, using new three-dimensional reconstructions of such folds in the onset region of the Greenland Petermann Glacier, the authors show that these formed due to flow convergence and the high mechanical anisotropy of ice.

West Africa was subjected to deformation and exhumation in response to Gondwana break-up. The timing and extent of these events are recorded in the thermal history of the margin. This study reports new apatite fission track (AFT) data from Palaeoproterozoic basement along the primary NE-SW structural trend of the Bole-Nangodi shear zone in northwestern Ghana. The results display bimodality in AFT age (populations of ~210-180 Ma and ~115-105 Ma) and length distributions (populations of 12.2 ± 1.6 and 13.1 ± 1.4 µm), supported by differences in apatite chemistry (U concentrations). The bimodal AFT results and associated QTQt thermal history models provide evidence for multiple cooling phases. Late Triassic – Early Jurassic cooling is interpreted to be related with thermal relaxation after the emplacement of the Central Atlantic Magmatic Province (CAMP). Early to middle Cretaceous cooling is thought to be associated with exhumation during the Cretaceous onset of rifting between West Africa and Brazil. Late Cretaceous – Cenozoic cooling can be related with exhumation of the Ivory Coast – Ghana margin and NNW-SSE shortening through western Africa. Furthermore, our data record differential exhumation of the crust with respect to the Bole-Nangodi shear zone, preserving older (CAMP) cooling ages to the south and younger (rifting) cooling ages to the north of the shear zone, respectively. This suggests that the Palaeoproterozoic BN shear zone was reactivated during the Cretaceous as a result of deformation in the Equatorial Atlantic region of Africa.

The discovery and identification of geophysical features from diverse gridded datasets play a pivotal role in understanding geological phenomena. Traditional tools tailored to identify specific signatures, such as lineaments in magnetic data, do not account for (1) the naturally occurring complexity of the causative geology which results from competing and overprinting processes and (2) often overlook important related patterns present in other forms of data (e.g., gradients in the associated gravity data). To address those limitations, we propose a two-step data-driven approach that integrates diverse gridded datasets and autonomously reveals inherent patterns without predefined assumptions in terms of feature geometry. Utilizing Haralick texture descriptors, we encode data patches from different geophysical datasets as points in a high-dimensional unified representation space, facilitating seamless data fusion. This representation is then nonlinearly projected into a two-dimensional space using t-distributed stochastic neighbor embedding (t-SNE), forming an interactive “t-SNE Atlas”. This atlas positions tiles from each dataset according to their position in the t-SNE while linking them back to their original geographical coordinates. Such a setup allows earth scientists to intuitively navigate through the data, exploring complex relationships between geophysical responses and geological structures, thus facilitating the discovery of new insights, the formulation of hypotheses, and the exploration of non-trivial connections. This is illustrated using magnetic and gravity data covering in the east Yilgarn Craton, Western Australia. Our methodology is adaptable and can be extended to include various other gridded datasets like gamma-ray spectrometry data, satellite imagery, and digital terrain models, thus broadening its applicability and enhancing geoscientific exploration. Our approach not only reveals subtle geophysical patterns but also offers practical benefits for exploration geologists. The atlas enables the identification of geological settings associated with mineral prospectivity and can serve as a pre-planning tool for pinpointing promising exploration sites, potentially accelerating discovery timelines.

3D thermo-mechanical models have been constructed to explore the influence of pre-existing faults on the temporal-spatial distribution of high-grade amphibolite-granulite facies metamorphic rocks during extension. The different dip amounts of three pre-existing faults (refer to the Bole-Nangodi shear zones, the Jirapa shear zone and the Bulenga shear zone in the study area) are studied in models. The results show the lower crust exhumes occur at a relatively low rate prior to attaining a Stretching Factor = 4.2%. The partially molten lower crustal rocks tend to move (up to 20 km) towards the center of the model and focus exhumation in regions where pre-existing faults intersect. The high-strain corridors in models are used to understand the loci of exhumation in the Bole-Bulenga domain of NW Ghana. Accordingly, in the eastern and western parts of the high-grade rock corridors in NW Ghana, partially molten rocks exhumed from the lower into middle-upper crustal levels are interpreted to have been dominantly facilitated by the km-scale high-strain corridors. In the central part of the Bole-Bulenga domain, the high-grade rocks are interpreted to have been exhumed because of a coupling between two mechanisms: (1) The exhumation of partially molten rocks between the Jirapa and Bole-Nangodi faults increases in spatially due to the reduction in space from north to south. (2) The exhumation of lower partially molten rocks in the central part, as a result of inherited orthogonal (E-W) structures.

We propose, test and apply a methodology integrating 1D magnetotelluric (MT) and magnetic data inversion, with a focus on the characterisation of the cover–basement interface. It consists of a cooperative inversion workflow relying on standalone inversion codes. Probabilistic information about the presence of rock units is derived from MT and passed on to magnetic inversion through constraints combining structural constraints with petrophysical prior information. First, we perform the 1D probabilistic inversion of MT data for all sites and recover the respective probabilities of observing the cover–basement interface, which we interpolate to the rest of the study area. We then calculate the probabilities of observing the different rock units and partition the model into domains defined by combinations of rock units with non-zero probabilities. Third, we combine these domains with petrophysical information to apply spatially varying, disjoint interval bound constraints (DIBC) to least-squares magnetic data inversion using the alternating direction method of multipliers (or ADMM). We demonstrate the proof-of-concept using a realistic synthetic model reproducing features from the Mansfield area (Victoria, Australia) using a series of uncertainty indicators. We then apply the workflow to field data from the prospective mining region of Cloncurry (Queensland, Australia). Results indicate that our integration methodology efficiently leverages the complementarity between separate MT and magnetic data modelling approaches and can improve our capability to image the cover–basement interface. In the field application case, our findings also suggest that the proposed workflow may be useful to refine existing geological interpretations and to infer lateral variations within the basement.

We present a novel methodology for performing experiments with subsurface structural models using a set of flexible and extensible Python modules. We utilize the ability of kinematic modelling techniques to describe major deformational, tectonic, and magmatic events at low computational cost to develop experiments testing the interactions between multiple kinematic events, effect of uncertainty regarding event timing, and kinematic properties. These tests are simple to implement and perform, as they are automated within the Python scripting language, allowing the encapsulation of entire kinematic experiments within high-level class definitions and fully reproducible results. In addition, we provide a link to geophysical potential-field simulations to evaluate the effect of parameter uncertainties on maps of gravity and magnetics. We provide relevant fundamental information on kinematic modelling and our implementation, and showcase the application of our novel methods to investigate the interaction of multiple tectonic events on a pre-defined stratigraphy, the effect of changing kinematic parameters on simulated geophysical potential fields, and the distribution of uncertain areas in a full 3-D kinematic model, based on estimated uncertainties in kinematic input parameters. Additional possibilities for linking kinematic modelling to subsequent process simulations are discussed, as well as additional aspects of future research. Our modules are freely available on github, including documentation and tutorial examples, and we encourage the contribution to this project.

High-grade (amphibolite–granulite facies) tectono-metamorphic domains in the Sefwi terrane of Ghana are separated from adjacent lower-grade (greenschist facies) greenstone belts by two main shear zones. The high-grade rocks presumably exhumed along the sinistral shear zones during the D2 ENE-WSW transtension (~2073 Ma). To better understand the role boundary conditions and the spatial relationship of faults play in the exhumation of partially molten lower crust in the Sefwi terrane, ten 3D thermomechanical models have been constructed. The results show that the normal component of velocity boundary conditions mainly controls the exhumation (8–10 km) of the lower crust along pre-existing faults, while the exhumation in the relay zones between faults is controlled by the obliquity between the applied extensional velocity vector and the vertical wall on which it is applied. The strike of the exhumation belt made of partially molten lower crust rocks in the relay zone is sub-orthogonal to the horizontal maximum stretching axis. The isostatic compensation from low-density upper mantle to overlying crust (thinning) is higher under transtension than under extension. The lower crust exhumation influenced by inherited shear zones (ductile) can be used to better understand the loci of the high-grade rocks in the Sefwi terrane. We suggest that the Kukuom-Juaboso domain composed of amphibolite–migmatite facies rocks probably resulted from the concentration of partially molten rocks in the relay zone between the Ketesso and Kenyase shear zones during the D2 ENE-WSW transtension. The two shear zones probably underwent two main stages for growth and maturation from the D1 to D2 phases. The regional exhumation of the high-grade rocks in the Sefwi terrane probably occurred within < 5 Ma.

Australian commonwealth, state and territory geological surveys possess information on over 3 million drillhole logs. In addition to mineral exploration drilling, extensive drillhole datasets exist from oil and gas exploration and hydrogeological studies. Other countries no doubt have similar data holdings. Together these legacy drillhole datasets have the potential to significantly enhance constraints on regional 3D geological models and improve our understanding of subsurface architecture, but have limited use in their current form as many if not most drill logs lack stratigraphic information, containing only lithological descriptions. This study develops open-source codes and methodologies for stratigraphy recovery (determining the ordered sequence of stratigraphic units) from drillhole lithological data by introducing a search algorithm that systematically explores all geologically plausible stratigraphic orderings for individual drillholes, combined with a solution correlation algorithm that compares the topological relationships of stratigraphic units across multiple drillholes to identify geologically consistent solutions and reduce uncertainty. The algorithms combine constraints from lithological descriptions with stratigraphic relationships automatically derived from regional maps. In addition, the method quantifies uncertainty by generating multiple plausible stratigraphic interpretations, providing critical insights for resource estimation, scenario analysis, and data acquisition strategies. The application of our method to a dataset of 52 drillholes from South Australia demonstrated its ability to make useful predictions of stratigraphic solutions and quantifying associated uncertainties. These results not only validate our approach but also highlight opportunities to refine current stratigraphic descriptions and provide a valuable new source for regional 3D geological modelling.

In this contribution we introduce LoopStructural, a new open-source 3D geological modelling Python package (http://www.github.com/Loop3d/LoopStructural, last access: 15 June 2021). LoopStructural provides a generic API for 3D geological modelling applications harnessing the core Python scientific libraries pandas, numpy and scipy. Six different interpolation algorithms, including three discrete interpolators and 3 polynomial trend interpolators, can be used from the same model design. This means that different interpolation algorithms can be mixed and matched within a geological model allowing for different geological objects, e.g. different conformable foliations, fault surfaces and unconformities to be modelled using different algorithms. Geological features are incorporated into the model using a time-aware approach, where the most recent features are modelled first and used to constrain the geometries of the older features. For example, we use a fault frame for characterising the geometry of the fault surface and apply each fault sequentially to the faulted surfaces. In this contribution we use LoopStructural to produce synthetic proof of concepts models and a 86 km × 52 km model of the Flinders Ranges in South Australia using map2loop.