Geological Applications Research Articles

Derivation of site-specific soil-water characteristic curve (SWCC) from extremely sparse experimental data by hierarchical Bayesian method with consideration of geotechnical sites similarity

Soil-water characteristic curves (SWCCs) play a crucial role in understanding soil behavior related to water movement and soil moisture effects, rendering them an essential tool in engineering geology and geotechnical engineering applications. Traditionally, SWCCs are determined through labor-intensive laboratory experiments involving varying levels of suction, a process that can take several months. Moreover, obtaining a high-quality SWCC from numerous measurements becomes particularly challenging when immediate site-specific data are required for design and analysis. This paper introduces a hierarchical Bayesian method for deriving site-specific SWCCs by integrating extremely sparse data (e.g., one or two measurements) for the site of interest with existing data from sites with similar geological and sedimentary characteristics. The SWCC parameters are estimated using a Bayesian framework and Markov chain Monte Carlo simulations. This approach not only enables the derivation of accurate SWCCs but also helps quantify the associated uncertainties. The effectiveness of the proposed method is demonstrated using both numerical and real-world data from different types of loess and unsaturated soils in the unsaturated soil database (UNSODA). The results show that site-specific SWCCs of unsaturated soils can be accurately estimated from sparse measurements by incorporating information from similar sites. This work offers an efficient and reasonably accurate approach for deriving SWCCs of unsaturated soils for geotechnical applications, especially when the number of site-specific SWCC measurement is extremely sparse and limited.

Experimental and numerical analyses of creep-fatigue behavior in sandstone for energy storage applications

Porous sandstone formations are recognized as suitable hosts for the geological storage of clean energy sources such as hydrogen and compressed air. From a geomechanical perspective, the deformation mechanisms of the host rock must be understood under various loading and unloading conditions resulting from daily or hourly injection/production cycles. This study presents a combined experimental and numerical analysis of sandstone behavior under cyclic loading conditions. We present the results of several multi-level, multi-frequency cyclic loading experiments conducted on sandstone specimens. Three distinct deformation regimes are identified in these tests: 1) instantaneous elastic deformation; 2) transient and steady-state strain accumulation associated with creep; and 3) accelerated deformation and fatigue failure at high stress levels caused by the progressive development of micro-cracks. Furthermore, the experiments show that fatigue failure is accompanied by a rapid stiffness degradation within the specimens. Based on these observations, we propose a viscoelastic-damage model to characterize the creep-fatigue behavior in sandstone. The model combines the standard Burger's viscoelasticity with an energry-driven energy-driven damage model. The performance of the model is verified against results obtained from laboratory experiments. Finally, a parametric analysis is presented to assess the effects of the frequency and magnitude of the cyclic loads on the fatigue life of sandstone. These findings indicate that increasing the magnitude and frequency of cyclic stress can accelerate fatigue failure. These insights can be used to optimize design parameters to maintain the stability of sandstone rock masses in geological energy storage applications.

Ultra-wide viewing angle holographic display system based on spherical diffraction

In the field of holographic displays, achieving a large viewing angle (LVA) has long been an essential and formidable challenge due to the limited diffraction angle of planar spatial light modulators (SLMs). Spherical holography, which allows observation from all perspectives, represents a practical approach to achieve the LVA. However, the physical limitations of commercial SLMs constrain direct implementation, making the realization of spherical holographic displays almost impractical and thus hindering technological advancement. This paper introduces an optical solution for spherical holographic display system, enabling the accurate reproduction of ultra-wide viewing angles and large objects. Our system utilizes a novel technique involving a convex parabolic mirror to convert an incident plane wave into a spherical wave. This innovative strategy permits the use of a planar SLM to create a 360°spherical holographic display. Additionally, we address the sampling issue for generating spherical holograms by constraining the point spread function (PSF), reducing the number of samples, and adapting it for visible light. Numerical simulations and optical experiments validate our success in achieving an ultra-wide viewing angle with a 360°horizontal angle and an 80°zenith angle range. Our system outperforms the state-of-the-art holographic-optical-elements approach both in computation speed and object size. We believe that our work significantly advances spherical holographic displays and offers a practical solution with vast potential applications in medicine, geology, and astronomy.

IMAGE-BASED YOLOV4 ARCHITECTURE FOR DETECTING MINERAL IN SEDIMENTARY ROCKS THIN SECTION

Thin section analysis of sedimentary rocks is the basis for identifying minerals and textures. In general, quantitative analysis of thin sections of rock often requires many hours of work when done manually. In today's era, mineralogical interpretation and percentage calculations must be carried out automatically using more practical applications. The research method begins with the identification of 44 thin section samples in parallel plane polarized (PPL) and crossed polarized (XPL) conditions with thin section analysis then mineralogy detection is carried out using a computational approach, namely the use of image-based Deep Learning YoloV4 architecture with 2D RGB image objects from the thin section of sedimentary rock. The results of this study show the best values of Average Precision in Quartz, Feldspar, and lithic are 39.21% in the XPL model, 26.53% in the XPL model, and 15.75% in combined mode, according to the training and testing of YoloV4 Models for the identification of rock minerals in thin sections. Based on the complexity of the mineral types, the granularity of the detection, and the specific geological objectives, establishing a meaningful benchmark or baseline for comparison is always challenging. Additionally, consider discussing the trade-offs between precision and recall, as a higher precision may be more critical in some geological applications. It is expected that the application of this research can produce practical, fast and accurate interpretation of the determination of minerals in sedimentary rocks from all thin-section images of rocks and thus provide a complete understanding of geological views automatically.

Depth dimension in seepage detection: Insights for exploration and geological gas storage surveillance

The search for surface hydrocarbon seeps has historically served as a fundamental tool in the early stages of oil and gas exploration. However, as easily detectable surface hydrocarbon seeps become increasingly scarce, advancement in analytical capabilities have shifted the focus towards detecting ultra-trace levels of hydrocarbons, commonly referred to as microseepages. While it is critical to understand the distribution of hydrocarbon microseepage for assessing petroleum systems and optimizing exploration efforts, the role of sampling depth in the detection of microseepage remains poorly understood. In this study, we investigated the influence of sampling depth on microseepage detection, using passive soil-gas sampling techniques at depths ranging from 0.5 to 10 m at sites with hydrocarbon-bearing reservoirs and dry wells. Organic geochemical analysis, including thermal desorption-gas chromatography-mass spectrometry (GC-MS), was employed to identify and quantify hydrocarbon and nonhydrocarbon compounds.The study result demonstrates that sampling at the Minimal Interference Zone (MIZ) depths of 5 m or greater significantly mitigates the risk of interference signals, thereby enabling the identification and quantification of hydrocarbon microseepages at trace and ultra-trace concentrations. At these depths, the attenuation of non-hydrocarbon compounds was observed, making the hydrocarbon signal markedly more distinguishable. Geochemical characterization of the samples revealed spatial and vertical heterogeneity in the distribution of hydrocarbon compounds, indicative of variations in the microseepage mechanisms potentially influenced by factors such as sedimentary basin settings and secondary alterations associated with the vadose zone. By utilizing passive soil-gas sampling techniques coupled with GC-MS, this study emphasizes the importance of sampling depth to improve the detection of reliable microseepage signal while reducing interference and background noise. The proposed technique not only refines existing methodologies for petroleum system assessments but also paves the way for broader geological applications. Specifically, the approach offers critical insights for the development of high-resolution geochemical surveys adaptable for surface background profiling in CO2 sequestration projects to monitor for storage efficiency and seal integrity.

Technical note: Optimizing the in situ cosmogenic 36Cl extraction and measurement workflow for geologic applications

Abstract. In situ cosmogenic 36Cl analysis by accelerator mass spectrometry (AMS) is routinely employed to date Quaternary surfaces and assess rates of landscape evolution. However, standard laboratory preparation procedures for 36Cl dating require the addition of large amounts of isotopically enriched chlorine spike solution; these solutions are expensive and increasingly difficult to acquire from commercial sources. In addition, the typical workflow for 36Cl dating involves measuring both 35Cl/37Cl and 36Cl/Cl concurrently on the high-energy (post-accelerator) end of the AMS system, but 35Cl/37Cl determinations using this technique can be complicated by isotope fractionation and system memory during measurement. The traditional workflow also does not provide 36Cl extraction laboratories with the data needed to calculate native Cl concentrations in advance of 36Cl/Cl measurements. In light of these concerns, we present an improved workflow for extracting and measuring chlorine in geologic materials. Our initial step is to characterize 35Cl/37Cl on sample aliquots of up to ∼1 g prepared in Ag(Cl, Br) matrices, which greatly reduces the amount of isotopically enriched spike solution required to measure native Cl content in each sample. To avoid potential issues with isotope fractionation through the accelerator, 35Cl/37Cl is measured on the low-energy, pre-accelerator end of the AMS line. Then, for 36Cl/Cl measurements, we extract Cl as AgCl or Ag(Cl, Br) in analytical batches with a consistent total Cl load across all samples; this step is intended to minimize source memory effects during 36Cl/Cl measurements and allows the preparation of AMS standards that are customized to match known Cl contents in the samples. To assess the efficacy of this extraction and measurement workflow, we compare chlorine isotope ratio measurements on seven geologic samples prepared using standard procedures and the updated workflow. Measurements of 35Cl/37Cl and 36Cl/Cl are consistent between the two workflows, and 35Cl/37Cl values measured using our methods have considerably higher precision than those measured following standard protocols. The chemical preparation and measurement workflow presented here (1) reduces the amount of isotopically enriched chlorine spike used per rock sample by up to 95 %; (2) identifies rocks with high native Cl concentrations, which may be lower priority for 36Cl surface exposure dating, at an early stage of analysis; and (3) allows laboratory users to maintain control over the total chlorine content within and across analytical batches. These methods can be incorporated into existing laboratory and AMS protocols for 36Cl analyses and will increase the accessibility of 36Cl dating for geologic applications.

Mixed Potential Sensors for the Detection of Hydrogen Leaks from Geologic Reservoirs

Large-scale hydrogen (H2) geologic storage presents a viable alternative to surface storage that is both safer and more economical. One option being considered by DOE is H2 storage in salt formations. These formations have been successfully used for gas storage for decades due to their extremely low permeability and may see significant investment as the US ramps up its R&D efforts in this field. New Mexico has a unique advantage in that bedded salt formations are located strategically in areas where excess energy production can be harnessed for hydrogen production. However, bedded salt is more permeable than domal salt, necessitating sensor deployment to monitor H2 containment. Robust sensors that can provide early leakage detection will be an essential component of the safety of H2 operating facilities. In a LANL-funded Mission Foundations Research (MFR) project, we prepared several zirconia-based mixed potential hydrogen sensors utilizing a tin-doped indium oxide electrode for H2 leak detection to test its stability in subsurface environments. Previous electrochemical sensor work has shown this type of sensor in a hydrogen safety sensor role represents a significant improvement over currently available sensors as it has an extremely stable baseline and limited cross sensitivity to other gaseous compounds.1 The feasibility of this technology will be demonstrated through a hydrogen permeation experiment to measure rates of underground hydrogen permeation between boreholes drilled underground at the DOE Waste Isolation Pilot Plant (WIPP) site. This experiment will (1) establish the use of this novel sensor in geologic H2 storage applications and (2) investigate the feasibility of H2 storage in bedded rock salt. Acknowledgement This research is supported by the Los Alamos National Laboratory, Laboratory Directed Research and Development (LDRD) Office, Missions Foundation Research Charter (Debora Wagner Program Manager). References E. L. Brosha, C. J. Romero, D. Poppe, T. L. Williamson, C. R. Kreller, R. Mukundan, R. S. Glass, and A. S. Wu, “Field Trials Testing of Mixed Potential Electrochemical Hydrogen Safety Sensors at Commercial California Hydrogen Filling Stations,” Journal of the Electrochemical Society 164 (13) B681-B689 (2017).

Modeling hydrogeological conditions of noble gases in water and aqueous electrolyte solutions

Dissolved noble gases in water hold significant potential as tracers and precursors in various hydrogeological and geological applications, particularly in understanding surface water-groundwater interactions. These applications necessitate precise solubility predictions. Employing the fugacity-fugacity approach, our study models the solubilities of Neon, Argon, Krypton, Xenon, and Radon in both pure water and electrolyte solutions. Leveraging the Cubic-Plus-Association Equation of State, we achieve remarkable accuracy, with average deviations below 2.5 % for vapor pressure, liquid, and vapor-phase densities. Additionally, in water–gas binary systems, deviations in gas solubilities within pure water remain under 5.0 %. Extending our analysis to water–gas-salt ternary systems, our model demonstrates deviations smaller than 5.6 %. Alongside discussing the gas dissolution mechanism, this work emphasizes the adaptability of the model, informed by robust experimental and modeling data. The interpretations derived from modeling hydrogeological conditions are poised to enhance our understanding of hydrological processes, thereby assisting planners and water managers in achieving sustainable development and management of water resources.

Sorption Efficacy of Weathered Basalt Rock for Metal Ions of Nuclear Importance

Background:: A safe and long-term isolation of radionuclides is crucial for efficient waste management in the nuclear fuel cycle. Clay minerals of geological origin are considered candidate barrier materials for the sequestration of nuclear waste for geological disposal applications. Objective:: The present study encompasses a systematic characterization of weathered basalt rock and the influence of its constituent clay minerals (formed mainly due to the weathering process) in the uptake of metal ions of nuclear importance, such as cesium (Cs+) and strontium (Sr2+). Methods: The structural profile of the weathered basalt has been investigated using different analytical techniques, including polarizing microscope, XRD, FTIR and EDXRF. The sorption behaviour of the rock sample for Cs+ and Sr2+ has been investigated in a comprehensive batch mode by varying the experimental conditions. The analytical findings for structure and batch sorption performance of the material have been further correlated to understand the influence of different parameters on the uptake of metal ions and the underlying mechanism. Results:: Structural analyses confirmed the presence of clay minerals viz., kaolinite, illite and montmorillonite in the sample. A comprehensive sorption performance assessment carried out in batch mode at different experimental conditions revealed that the uptake of both the metal ions was rapid and dependent on initial metal ion concentration and solution pH. The uptake of Cs+ ions was found to be higher as compared to the Sr2+ ions. EDXRF analysis confirmed the loading of Cs+ and Sr2+ on the weathered basalt. Experimental batch sorption data presented a better agreement with the theoretical Freundlich isotherm pertaining to the heterogeneous nature of the sorbent. Conclusion:: The studies highlight that the clay minerals formed by structural alteration of basalt rock upon intense weathering could be very useful in fixing the nuclear fission waste components such as Cs+ and Sr2+.

Rheological measurements and transition to turbulence for moderate Reynolds number inertial suspensions

Particulate flows at moderate particle Reynolds numbers are important in critical engineering and geological applications. This experimental study explores neutrally buoyant suspensions in an outer-rotating coaxial rheometer for solid fractions, $\phi$ , from 0.1 to 0.5, and particle Reynolds number, $Re$ , from 0.5 to 800, covering laminar, transitional and turbulent regimes; $Re$ is defined in terms of the square of the particle diameter and the shear rate. For $0.1 < \phi < 0.4$ and $0.5 < Re <10$ , the direct torque measurements normalised by the laminar flow torque, $M/M_{lam}$ , are independent of $Re$ , but depend on $\phi$ . For the same range of $\phi$ and for $10< Re<100$ , the normalised torques depend on both $\phi$ and $Re$ , and show an increasing dependence on $Re$ . As $Re$ increases, the flow transitions to turbulence. Small particles delay the turbulent transition for $\phi \leqslant 0.3$ , while large particles augment the transition. A modified Reynolds number, $Re^\prime$ , that depends linearly on the particle diameter and the maximum velocity, $U_{o}$ , is introduced for both laminar and turbulent flows and shows a better correlation of the results as compared with $Re$ . For $\phi = 50\,\%$ , the normalised torque minus the torque at zero rotational speed is nearly independent of $Re^\prime$ . Rheological models based on $Re^\prime$ and the Krieger–Dougherty relative viscosity are proposed in the laminar regime for $10< Re^\prime <500$ ; in the turbulent regime, a correlation is proposed in terms of $Re^\prime$ and $\phi$ for $1000< Re^\prime < 6000$ .

The Theoretical Calculation of the Cu Isotope Fractionation Effect in Solution/Hydrothermal Solution Systems.

Copper (Cu) is an important transition metal, and its isotopes have important applications in geology, environmental science, soil science, and other fields. Cu isotope fractionation can occur in many natural processes. However, the mechanism of Cu isotope fractionation in solution/hydrothermal solution systems is not very clear. In this study, the fractionation effects of complexes of Cu(I) and Cu(II) in solution/hydrothermal solution systems were systematically studied by means of an ab initio method based on first principles. In the simulation of an aqueous solution system, the theoretical treatment method used is the "water-droplet" method. The results show that the heavy Cu isotope (65Cu) enrichment capacity of the Cu-bearing complex solutions is greatly affected by the ligand types both for Cu(I) and Cu(II). For Cu(I) complex solutions, the heavy Cu isotope enrichment sequence is [Cu(HS)2]-·(H2O)42 > [Cu(HS)(H2O)]·(H2O)42 ≈ [Cu(HS)(H2S)]·(H2O)42 > [CuCl]·(H2O)42 > [CuCl2]-·(H2O)42 > [CuCl3]2-·(H2O)42. For the aqueous solutions of Cu(II) with an inorganic ligand (such as H2O, OH-, NO3-, SO42- and CN-), the order of heavy Cu isotope enrichment is as follows: [Cu(H2O)6]2+·(H2O)42 > [Cu(NO3)2]·(H2O)42 > [Cu(OH)2]·(H2O)42 > [CuSO4(H2O)3]·(H2O)42 > [CuNO3(H2O)4]+·(H2O)42 > [CuCN]+·(H2O)42. For the Cu(II) complex solutions with a halogen as ligands, the change order of 1000lnβ is [CuCl]+·(H2O)42 > [CuCl2]·(H2O)42 > [CuBr2]·(H2O)42 > [CuCl3]-·(H2O)42. The sequence of 1000lnβ for Cu(II) organic complex aqueous solutions is [Cu(HOC6H4COO)]+·(H2O)42 > [Cu(CH3CH2COO)]+·(H2O)42 > [Cu(COOHCOO)]+·(H2O)42. The calculation also found that for Cu(I) complex aqueous solutions, the difference in Cu isotope fractionation parameters (1000lnβ) between [CuCl2]-·(H2O)42 and [Cu(HS)2]-·(H2O)42 is relatively large. At 100 °C, the 1000lnβ of the two species are 1.14 and 1.55 (‱), respectively. The difference between the two could be reached up to 0.41 (‱). The Cu isotope fractionation parameter obtained with the "water droplet" method is also very different from the results of previous studies, which indicate that the Cu isotope fractionation behavior of the two is similar. At the same time, the exciting discovery is that the enrichment capacity of heavy Cu isotopes is significantly different between Cu(I) complex aqueous solutions and Cu(II) complex aqueous solutions. At 100 °C, the 1000lnβ of 6 Cu(I) complex aqueous solutions and 13 Cu(II) complex aqueous solutions ranged from 0.90 to 1.55 and 2.24 to 3.25(‱), respectively. It also shows that the REDOX reaction has a significant effect on the Cu isotope fractionation, especially in ore-forming fluids. Therefore, the ligand type is a factor that cannot be ignored when considering the mechanism of Cu isotope fractionation in solution/hydrothermal solution systems. Whether the solvation effect of an aqueous solution is considered or not has a great influence on the numerical values of the final Cu isotope fractionation factors. Hence, the solvation effect of an aqueous solution is an essential determinant in the theoretical calculation of the Cu isotope fractionation factors for Cu-bearing complex solutions.

Macro-micro mechanical behavior of saturated cemented sands during drained triaxial shearing

To enhance the understanding of the macro-micro mechanical behavior of cemented sands, a series of drained triaxial shearing tests were conducted on saturated cemented sands with varying cement contents and confining stresses. These tests incorporated high-performance acoustic emission (AE) monitoring and an environmental scanning electron microscope (ESEM) technique. The results yielded quantitative interpretations of macro-mechanical response, microstructural variation, and AE behaviors. Specifically, increasing cement content led to an exponential reduction in porosity (or pore area ratio) and a power increase in pore number, possibly due to the production of hydration products that segment large pores. Relationships were established among peak strength (or AE hits), cement content, and confining stress. However, AE trends diverged from those observed in peak strength due to microstructural modifications under different cement contents and their resulting differences in micromechanical behavior during shearing. General relationships between AE rate and mechanical parameters were established during the initial shearing stage, regardless of cement content and confining stress. Moreover, AE evolution efficiently characterized subtle variations in failure patterns due to its high sensitivity to failure-induced structural modifications. These contributions have the potential to establish a non-destructive field evaluation framework and provide guidelines for the efficient utilization of cemented soils in geological and geotechnical engineering applications.

Knudsen effusion mass spectrometry: Current and future approaches.

Knudsen effusion mass spectrometry (KEMS) has been a powerful tool in physical chemistry since 1954. There are many excellent reviews of the basic principles of KEMS in the literature. In this review, we focus on the current status and potential growth areas for this instrumental technique. We discuss (1) instrumentation, (2) measurement techniques, and (3) selected novel applications of the technique. Improved heating methods and temperature measurement allow for better control of the Knudsen cell effusive source. Accurate computer models of the effusive beam and its introduction to the ionizer allow optimization of such parameters as sensitivity and removal of background signals. Computer models of the ionizer allow for optimized sensitivity and resolution. Additionally, data acquisition systems specifically tailored to a KEMS system permit improved quantity and quality of data. KEMS is traditionally utilized for thermodynamic measurements of pure compounds and solutions. These measurements can now be strengthened using first principles and model-based computational thermochemistry. First principles can be used to calculate accurate Gibbs energy functions (gefs) for improving third law calculations. Calculated enthalpies of formation and dissociation energies from ab initio methods can be compared to those measured using KEMS. For model-based thermochemistry, solution parameters can be derived from measured thermochemical data on metallic and nonmetallic solutions. Beyond thermodynamic measurements, KEMS has been used for many specific applications. We select examples for discussion: measurements of phase changes, measurement/control of low-oxygen potential systems, thermochemistry of ultrahigh-temperature ceramics, geological applications, nuclear applications, applications to organic and organometallic compounds, and thermochemistry of functional room temperature materials, such as lithium ion batteries. We present an overview of the current status of KEMS and discuss ideas for improving KEMS instrumentation and measurements. We discuss selected KEMS studies to illustrate future directions of KEMS.

Identification of complex crude oil categories in ultra-deep petroleum systems based on a new fluorescence analysis method: An example from the southwestern Tabei Uplift, Tarim Basin, China

Ultra-deep crude oils have become one of the most important targets for petroleum exploration. In ultra-deep petroleum systems, differential hydrocarbon generation of various source rocks and multistage oil accumulation and secondary alteration processes led to the coexistence of crude oils with diverse sources and formation mechanisms. How to effectively identify complicated crude oil categories is a difficult issue confronted with ultra-deep petroleum exploration. Because crude oils are generally fluorescent, fluorescence analysis methods exhibit great application potential for classifying ultra-deep crude oils. However, due to the lack of other restrictive parameters, previous methods based on fluorescence color and fluorescence spectra limit the geological application of fluorescence analysis methods in ultra-deep petroleum systems. Based on the combination of fluorescence spectral parameters and fluorescence lifetimes, this study innovatively provides a new fluorescence analysis method for identifying complex crude oil categories in ultra-deep petroleum systems. Moreover, the sources and formation mechanisms of the ultra-deep crude oil can be further identified by combining fluorescence parameters with physical and geochemical parameters. By applying the above method, this study effectively identified three categories of crude oils and revealed their sources and formation mechanisms in the ultra-deep petroleum system of the southwestern Tabei Uplift, Tarim Basin, China. The new fluorescence analysis method may have practical significance for revealing the sources and formation mechanisms of ultra-deep crude oils.

Fabrication of magnetic field sensor based on a Co-based Schottky heterojunction

We have prepared a reliable magnetic field sensor based on a Co-based Schottky heterojunction with the possibility of finding practical application when the current–voltage (I–V) characteristics of the heterojunctions were modified by the external magnetic fields. For this purpose, the heterojunctions based on ferromagnetic thin films Co were deposited on Si (100) substrate via ultra-thin Pt films. The structural and magnetic properties of the heterojunctions were demonstrated, the crystal structure and the saturation magnetization of heterojunctions significantly increased with the increase of annealing temperatures. The electrical properties of the Co/Pt/Si heterojunction were investigated using temperature-dependent I–V characteristics, the rectifying properties were observed in the whole temperature ranging from 233 K to 283 K. Based on a simplified energy band structure, the magnetic field transport properties of Co/Pt/Si heterojunctions were investigated by taking into account the effects of spin and magnetic field on the potential barrier height of Co/Pt/Si heterojunctions. The results showed that the experimental ratio of current measured in a magnetic field B and in zero magnetic field (I(B)/I(0)) is linearly with the magnetic field B and almost independent of the applied voltage. The sensitivity and linearity of the magnetic field sensor were measured. The sensitivity of the sensor can reach 72.09 μA/T at an applied voltage frequency of 2000 Hz, which indicates that the linearity and sensitivity of the sensor are improved with the increase of the applied voltage frequency. The proposed sensors would be easy to integrate with current semiconductor processes with simple structure and they could be suitable for geological applications such as energy exploration.

Fault-sketch: A framework for modeling geological faults and displacements

Geological faults, which result from rock fracturing and displacement in the subsurface, play an important role in shaping the Earth’s crust. Modeling fault surfaces and their impact on rock formations is, thus, useful for various geological applications. However, due to the complex configurations of faults, the modeling process can become difficult, requiring a framework with a specialized modeling toolkit.The primary challenge in fault modeling lies in properly correlating the geometry of faults with their effects on rock layers and other fault surfaces during displacement. To address this, we have developed a series of operators designed for fault modeling. In this manner, we present Fault-Sketch, a Sketch-Based Interface and Modeling (SBIM) framework that tackles the task of fault modeling and captures their influence on various geological structures. This framework enables the creation of surfaces as a result of rock configurations that have previously been largely overlooked in the existing literature.The collection of operators presented in Fault-Sketch are built upon a data structure fit for the purpose of consistently representing fault surfaces and their boundaries together with the influence they cause both on other surfaces and on the volumes they bound. Moreover, the operators are designed to mimic the geologic processes relevant to fault creation and displacement to facilitate their use by domain experts.

A Secular Solar System Resonance that Disrupts the Dominant Cycle in Earth’s Orbital Eccentricity (g 2 − g 5): Implications for Astrochronology

The planets’ gravitational interaction causes rhythmic changes in Earth’s orbital parameters (also called Milanković cycles), which have powerful applications in geology and astrochronology. For instance, the primary astronomical eccentricity cycle due to the secular frequency term (g 2−g 5) (∼405 kyr in the recent past) utilized in deep-time analyses is dominated by the orbits of Venus and Jupiter, i.e., long eccentricity cycle. The widely accepted and long-held view is that (g 2−g 5) was practically stable in the past and may hence be used as a “metronome” to reconstruct accurate geologic ages and chronologies. However, using state-of-the-art integrations of the solar system, we show here that (g 2−g 5) can become unstable over long timescales, without major changes in, or destabilization of, planetary orbits. The (g 2−g 5) disruption is due to the secular resonance σ 12 = (g 1 − g 2) + (s 1 − s 2), a major contributor to solar system chaos. We demonstrate that entering/exiting the σ 12 resonance is a common phenomenon on long timescales, occurring in ∼40% of our solutions. During σ 12-resonance episodes, (g 2−g 5) is very weak or absent and Earth’s orbital eccentricity and climate-forcing spectrum are unrecognizable compared to the recent past. Our results have fundamental implications for geology and astrochronology, as well as climate forcing, because the paradigm that the long eccentricity cycle is stable, dominates Earth's orbital eccentricity spectrum, and has a period of ∼405 kyr requires revision.

RupMC: a ray-unit parallel marching cubes algorithm on CPU/GPU heterogeneous architectures

ABSTRACT The marching cubes (MC) algorithm is widely used for extracting isosurfaces from volume data and 3D visualizations because of its effectiveness and robustness but require extensive memory and computing time for large-scale applications. Additionally, MC isosurfaces lack topologic information, making them difficult to use in some geologic applications. To overcome these limitations, this study proposes an enhanced MC using CPU/GPU heterogeneous architecture called the ray-unit parallel MC (rupMC) algorithm. First, ray units form the basic voxel to determine how the surface intersects to reduce repeated computations and enhance efficiency. Then, rupMC uses multiple computing processes and threads on a CPU/GPU heterogeneous architecture to process points concurrently. Finally, the unique surface intersection indices are preserved to compose the surface triangles, and the topological surface information is directly embedded in the triangle compositions. Experiments on five stratum datasets of varying sizes demonstrated that, rupMC achieved approximately dozens of times faster than other serial MC and 4 times faster than a parallel DMC. rupMC demonstrated high scalability and adaptability to various CPUs/GPUs and datasets of various sizes. rupMC has remarkable capabilities for efficiently and feasibly extracting precise surface intersections and triangles, making it well-suited for large-scale and high-density applications.

Assessing and Improving the Robustness of Bayesian Evidential Learning in One Dimension for Inverting Time-Domain Electromagnetic Data: Introducing a New Threshold Procedure

Understanding the subsurface is of prime importance for many geological and hydrogeological applications. Geophysical methods offer an economical alternative for investigating the subsurface compared to costly borehole investigations. However, geophysical results are commonly obtained through deterministic inversion of data whose solution is non-unique. Alternatively, stochastic inversions investigate the full uncertainty range of the obtained models, yet are computationally more expensive. In this research, we investigate the robustness of the recently introduced Bayesian evidential learning in one dimension (BEL1D) for the stochastic inversion of time-domain electromagnetic data (TDEM). First, we analyse the impact of the accuracy of the numerical forward solver on the posterior distribution, and derive a compromise between accuracy and computational time. We also introduce a threshold-rejection method based on the data misfit after the first iteration, circumventing the need for further BEL1D iterations. Moreover, we analyse the impact of the prior-model space on the results. We apply the new BEL1D with a threshold approach on field data collected in the Luy River catchment (Vietnam) to delineate saltwater intrusions. Our results show that the proper selection of time and space discretization is essential for limiting the computational cost while maintaining the accuracy of the posterior estimation. The selection of the prior distribution has a direct impact on fitting the observed data and is crucial for a realistic uncertainty quantification. The application of BEL1D for stochastic TDEM inversion is an efficient approach, as it allows us to estimate the uncertainty at a limited cost.

Introduction to this special section: Gravity, electrical, and magnetic methods

Geophysics plays an important role in all aspects of geologic analysis, spanning from regional exploration and tectonic mapping to local prospect-level studies and energy transition projects. Out of a variety of geophysical techniques, seismic methods often play the most central role, while less expensive nonseismic methods are less frequently applied. One of the objectives of SEG's Gravity and Magnetics Committee is to promote nonseismic geophysical methods and showcase their value in various geologic applications. As members of that committee, we, the editors of this special section, assert that gravity, magnetic, and electrical methodologies are powerful yet often undervalued tools. We present this special section focused on nonseismic geophysical methods and the impact they can make on various geoscience projects. Included here are a regional tectonic study in Antarctica, a local mining exploration mapping project in Canada, and an analytical methodology capable of providing critical information about subsurface geology and helping to determine the depths of important geologic structures. Although our scope is to provide a broader discussion of different types of nonseismic techniques, all of the papers presented in this special section focus on magnetic surveying.