Earth Model Research Articles

The need for ENVIRONMENTAL GEOPHYSICS: the case of abandoned complex behind UBA and FIDELITY banks near Greenhouse at Nsukka campus of University of Nigeria

Subsurface geologic structure around an abandoned building at the University of Nigeria, Nsukka campus was investigated to ascertain the cause or reason for abandonment. 2-D resistivity method was employed in the study. Three electrical resistivity tomography (ERT) profiles and one vertical electrical sounding (VES) profile were conducted around the building using SSR-MP-ATS resistivity meter. Winresist software and RES2DINV software were used to process the VES and ERT data respectively. The result of the analysis of the VES data showed a six layered earth model with the resistivity and depth respectively for each layer. The results of the ERT revealed that the area is covered with a conductive top-layer to a depth of about 15 m coupled with fault zone that passed through the building site. These results showed that the foundation of the building was laid without priori or adequate geologic information about the site before commencement of the building. Geologically, the areas with a conductive layer or fault zone have earth materials with low or poor engineering strength that could lead to collapse of engineering structure either by sinking or by fall and which is the case with this site we have studied.

Relative sea-level sensitivity in the Eurasian region to Earth and ice-sheet model uncertainty during the Last Interglacial

Fingerprinting the source and rate of the melt of polar ice sheets during the Last Interglacial is a key research challenge. This is reliant on high-quality relative sea-level constraints, and the correction of this data for the effects of glacial isostatic adjustment driven by ice sheet cover changes prior to the interglacial. However, both the spatial and temporal evolution of past ice sheets and the Earth’s rheological structure that serve as inputs to glacial isostatic adjustment predictions are significantly uncertain. This study sets out to determine the relative influence of each of these inputs on modelled values of Last Interglacial relative sea levels and how this influence varies spatially. To answer this question, we use a palaeo ice-sheet model and a gravitationally consistent glacial isostatic adjustment model. We develop new numerical tools to generate plausible ice-sheet extent and histories, quantify relative sea-level uncertainty, and perform a Sobol sensitivity analysis facilitated by the use of Gaussian process emulation. We find that Earth model parameters are the dominant contributors to relative sea-level uncertainty in most Eurasian regions, but that relative sea-level values in the Barents-Kara Sea are most influenced by ice-sheet loading, while the timing of the deglaciation has the greatest impact in the Baltic Sea. Our results show that the magnitude and rate of relative sea-level change is relatively insensitive to the specific timing of ice-sheet retreat, as well as the configuration of the far-field North American ice sheet. Overall, our work suggests that the coastlines of the southern North Sea and the English Channel are least influenced by relative sea-level uncertainty and are the most suitable for future data collection studies aiming to limit the influence of glacial isostatic adjustment.

Construction of Mechanical Earth Modeling for Mitigating Wellbore Instability in Tanuma Shaly Formation -Southern Iraqi Oil Field

Wellbore instability has become an increasing concern for horizontal and extended-reach wells, especially through shale cap rocks. This paper presents a comprehensive analysis of the factors that led to the wellbore being unstable in the Tanuma Formation, the Nasiriyah oilfield in southern Iraq. By constructing a mechanical earth model using open-hole wireline data from two drilled wells. Mohr-Coulomb, Mogi-Coulomb, and Modified Lade failure criteria are used to analyze the breakouts and determine the density of mud that is required for a safe wellbore wall. The results indicated that the Mogi-Coulomb criterion is more similar to the actual well failure than the other two criteria. The wellbore instability findings showed that the most stable wells are those that are vertical and have a low deviated angle (less than 10°), and it is advised that these wells be drilled with a mud weight of 11.3 to 14 ppg. The acceptable range for the weight of the mud for horizontal and greatly deviated wells is between 11.9 and 13.8 ppg. Consequently, the future wells development plan for the studied region is improved and leads to cost reductions.

Limited sensitivity of Antarctic GIA mass change estimates to lateral viscosity variations

The Gravity Recovery and Climate Experiment (GRACE) has revealed spatiotemporal mass changes in the Antarctic Ice Sheet. However, GRACE data must be corrected for the gravity changes due to glacial isostatic adjustment (GIA). Here we investigate the sensitivity of GIA-induced gravity changes in Antarctica to the lithospheric thickness and upper mantle viscosity using a one-dimensional (1-D) model that assumes a radially varying Earth structure. The sensitivity is assessed using several Antarctic ice-history models that have been widely used to correct GRACE data. The results indicate a trade-off between lithospheric thickness and upper mantle viscosity in evaluating the Antarctic GIA correction. This trade-off exists for all ice-history models; however, the reason for the trade-off differs among models. Furthermore, since there is a sharp contrast in the Earth structure between West and East Antarctica, the adopted ice histories are separated into West and East Antarctic components to examine their contributions to the Antarctic GIA correction. We consider 1-D Earth structures that are averaged from the seismically derived three-dimensional Earth structure for West and East Antarctica. These results indicate that the contributions of the East and West Antarctic loads do not significantly affect the GIA corrections for the West and East Antarctic regions, respectively, and that the trade-off between lithospheric thickness and upper mantle viscosity results in minimal divergence in the assessment of the Antarctic GIA correction between the averaged Earth models of West and East Antarctica. Therefore, the contrast in Earth structure beneath Antarctica may have a limited effect on the ice-mass change estimates for the entire Antarctic Ice Sheet.

Topography of the 410 and 660 km discontinuities beneath the Tibetan plateau and adjacent areas

To investigate the three-dimensional spatial distribution of subducted oceanic slab segments and their consequential effect on the thermal and hydrous composition of the mantle transition zone (MTZ) beneath the entire Tibetan Plateau, we conducted a comprehensive analysis of approximately 655,000 P-to-s receiver functions (RFs), obtained from 735 broadband seismic stations. These RFs were utilized to delineate the 410 km (d410) and 660 km (d660) seismic discontinuities, which represent the upper and lower boundaries of the MTZ, respectively. The RFs were grouped within circular bins with a radius of 1 degree and stacked to image the discontinuities. The mean apparent depths derived based on the 1-D IASP91 Earth model for the d410 and d660 across the entire study area are 412.2±8.3 km and 668.9±8.5 km respectively, and the MTZ thickness is 256.5±6.9 km. The observed apparent depths underwent subsequent correction utilizing multiple velocity models. Several areas in central Tibet exhibit a normal d410 and an anomalously deep d660, which can be attributed to the combined effect of the negative thermal anomaly and dehydration associated with subducted slab segments that have penetrated at least to the d660 depth. The anomalous thickening of the MTZ beneath the southeastern Tibetan Plateau surrounding the Tengchong volcanic field can be explained by the dehydration of the subducted Indian Slab. Significant thinning of the MTZ associated with the deepening of the d410 beneath the western Tian Shan may indicate active thermal upwelling originating from the MTZ.

Scaling-up dynamic elastic logs to pseudo-static elastic moduli of rocks using a wellbore stability analysis approach in the Marun oilfield, SW Iran

Wellbore stability analysis is a critical component of petroleum engineering, evaluating the risks of sanding, reservoir compaction, and casing failures. Laboratory rock mechanical measurements must be scaled up to reservoir scales to achieve accurate results. One challenge lies in upscaling dynamic measurements from petrophysical logs to pseudo-static elastic properties, which has significant implications for oil and gas operations. We present a novel approach that combines laboratory rock mechanical measurements with well-log data to develop a mechanical earth model (MEM) for an Iranian oilfield with over 350 wells. We conducted static elastic property measurements on 40 core samples from various layers and depths of carbonate and sandstone rocks, demonstrating the practical application of our approach. By integrating these measurements with dynamic log data and static-dynamic correlations, we established a framework for evaluating the mechanical properties of different layers. Our findings indicate that the safe mud weight window ranges from 41.5 to 118.59 pcf, while the stable mud weight window ranges from 41.5 to 156 pcf. We demonstrate the importance of conducting parallel rock mechanical studies on cores and logs to reduce uncertainties, costs, and risks during oil and gas operations. We also propose a novel methodology combining lithological characteristics, abnormally high pressure, and borehole instability mechanisms to evaluate the stability of borehole walls. This framework provides a fresh perspective on wellbore stability analysis and offers practical solutions for the industry. Essential novel techniques include developing a geomechanical model that integrates laboratory rock mechanical measurements with well-log data to evaluate mechanical properties and calculate safe and stable mud-weight windows. Our study advances wellbore stability analysis by providing a new method for addressing this long-standing challenge. It offers valuable insights for petroleum engineers working in the oil and gas industry.

Long-range LBL underwater acoustic navigation considering Earth curvature and Doppler effect

The accuracy of underwater acoustic navigation is significantly influenced by geometry configuration and systematic errors coming from geophysical environment. Long baseline (LBL) systems exhibit superior positioning precision than other acoustic navigation systems, primarily attributable to the better Geometric Dilution of Precision (GDOP). However, the geometry configuration tends to degrade when the underwater vehicle is far away from the seafloor beacons. Factors such as acoustic ray bending error, Doppler effect and Earth curvature may seriously affect the navigation accuracy. To address these issues, we present a robust Kalman filter with systematic error compensation for long-range LBL systems. Firstly, a reversed acoustic ray tracking method is proposed for a unique phenomenon in acoustic wave propagation. Following that, we construct an acoustic positioning model with a time bias parameter, to weaken the impact of Doppler effect on acoustic ranging. This additional parameter was estimated with the state of the user vehicle in real time. At last, an observation equation of pressure gauge considering Earth curvature is presented for depth constraint. Long-range LBL experiments conducted in the South China Sea is employed to validate the performance of the proposed methods, taking the positioning results of inverted ultra-short baseline (iUSBL) system as reference. The results show that the proposed method outperform traditional method with a decrease of about 60 % in the three-dimensional root mean square (3DRMS) of navigation errors. With mentioned three improvements, the 3DRMS of the long-range LBL system with the longest distance of 20 km from the beacon center is less than 14 m. These findings can provide theoretical basis for other long-range LBL systems.

Dataspace Integration for Agrobiodiversity Digital Twins with RO-Crate

A severe impact of human-induced climate change is the global decrease in crop yields, threatening the United Nations sustainable development goal (SDG) to end hunger, achieve food security, improve nutrition, and promote sustainable agriculture (SDG 2: “Zero Hunger Goal”). Decision-making on appropriate mitigation strategies requires data integration from heterogeneous data-rich contexts including Earth observation (EO), biodiversity monitoring, gene banks as well as socio-economic data. Accordingly, large-scale initiatives for thematic data spaces, such as the European projects Destination Earth (DestinE) and the Green Deal Data Space, have been launched. Underutilized crops and crop wild relatives (CWR) can contribute to mitigating the impacts of climate change since their greater genetic diversity allows adaptation of agro-food sources to extreme weather events such as flooding and drought or emerging crop pests. Within the framework of Biodiversity Digital Twin (BioDT), a prototype digital twin for crop genetic resources modeling is under development. This prototype combines the scanning of genetic resources of CWRs and traditional cultivars for promising pest resistance and abiotic tolerance traits with advanced species distribution modeling and actual climate scenarios from EO. The aim is to generate habitat suitability maps and identify potential areas where relevant populations are potentially growing (Chala et al. 2024). The talk gives an overview of data models used to facilitate the required interoperability of both Earth models and environmental data within and across different data spaces such as the DestinE Data Lake (Duatis Juarez et al. 2023). Based on our findings from the use case of agrobiodiversity digital twins, we will discuss approaches for processing machine-actionable digital objects (Fig. 1) containing workflow information (i.e., in the Common Workflow Language format, Crusoe et al. 2022), and for sharing resulting research data in a FAIR (Findable, Accessible, Interoperable, and Reusable) manner by leveraging on so-called webby FAIR Digital Objects involving Research Object Crate (RO-Crate, Soiland-Reyes et al. 2022), FAIR Signposting (Soiland-Reyes et al. 2024) and schema.org and its Bioschemas extension respectively (Castro et al. 2020),

Reactant Discovery with an Ab Initio Nanoreactor: Exploration of Astrophysical N-Heterocycle Precursors and Formation Pathways.

The incorporation of nitrogen atoms into cyclic compounds is essential for terrestrial life; nitrogen-containing (N-)heterocycles make up DNA and RNA nucleobases, several amino acids, B vitamins, porphyrins, and other components of biomolecules. The discovery of these molecules on meteorites with non-terrestrial isotopic abundances supports the hypothesis of exogenous delivery of prebiotic material to early Earth; however, there has been no detection of these species in interstellar environments, indicating that there is a need for greater knowledge of their astrochemical formation and destruction pathways. Here, we present results of simulations of gas-phase pyrrole and pyridine formation from an ab initio nanoreactor, a first-principles molecular dynamics simulation method that accelerates reaction discovery by applying non-equilibrium forces that are agnostic to individual reaction coordinates. Using the nanoreactor in a retrosynthetic mode, starting with the N-heterocycle of interest and a radical leaving group, then considering the discovered reaction pathways in reverse, a rich landscape of N-heterocycle-forming reactivity can be found. Several of these reaction pathways, when mapped to their corresponding minimum energy paths, correspond to novel barrierless formation pathways for pyridine and pyrrole, starting from both detected and hypothesized astrochemical precursors. This study demonstrates how first-principles reaction discovery can build mechanistic knowledge in astrochemical environments as well as in early Earth models such as Titan's atmosphere where N-heterocycles have been tentatively detected.

Spatial positioning method based on range-only measurement of multi-station radar

Aiming at the characteristics that affect the high ranging accuracy of shipborne radar, a multi station radar spatial positioning method based on ranging is proposed. This method is based on the ECEF (Earth Centered Earth Fixed) coordinate system, and adopts the method of intersecting the detection distances of multiple radars to the target, and directly calculates the spatial position of the target. This method only needs the geodetic coordinates of each radar station and the detection range of the target to calculate the geodetic coordinates of the target, which can effectively solve the coordinate conversion problem caused by the curvature of the earth, and does not involve higher-order equations and multi-value solutions, and it does not involve iterative initial value calculation. The positioning system can not only obtain high positioning accuracy, but also has a long operating distance, which can better solve the multi-radar range measurement intersection positioning problem.

Impact of long-term depletion on horizontal wellbore stability in tight reservoirs-including changes in petrophysical and geomechanical properties

The development of tight reservoirs using a horizontal well requires caution owing to the reservoir sensitivity to stress. Predicting the horizontal wellbore instability during reservoir depletion is critical for tight reservoirs. The present study investigated the wellbore failure criteria and stress orientation changes with reservoir depletion of a horizontal well in a tight reservoir, considering rock sensitivity to stress with pressure reduction. An integrated reservoir model coupled with geomechanics was used to predict pore pressure. A mechanical earth model (MEM) was constructed to examine the stability of horizontal wells. The unconfined compressive strength (UCS) was correlated with changes in porosity due to stress changes. The MEM shows that the failure criteria around the horizontal wellbore, including breakout, loss, and breakdown pressures, change considerably with reservoir depletion. A more significant response was observed in the stress-sensitive layer, which was characterised by a higher permeability. The safe mud weight window of the studied horizontal well narrowed significantly after production for five years, whereas the stability was completely absent after ten years. The deterioration of stability with depletion is governed mainly by permeability reduction, which may cause severe pressure reduction, whereas the stability improvement due to the increase in UCS caused by porosity reduction seems marginal. Relatively small changes in the direction of the stresses around the wellbore with depletion were observed, with the stress distribution concentrated in the predominant direction as the depletion continued. This study indicates that the sensitivity of petrophysical and geomechanical properties to stress amplifies the instability in tight reservoirs, particularly around wellbores.

Contourite-like deposits suggest stronger-than-present circulation in the Plio-Pleistocene Red Sea

Red Sea Deep Water is presently slow-moving, but was this true of the earlier Plio-Pleistocene (PP)? In seismic reflection records, the PP deposits are distorted by halokinetic deformation of their underlying Miocene evaporites. However, if reflections are flattened to a prominent reflector representing the top of the Miocene, they reveal mounded deposits within the earlier PP along both sides of the sea. Off Egypt, a plastered drift occurs along a salt wall. In the central Red Sea, they are mounded drifts. Seismic reflections from these deposits change shape gradually upwards to the modern seabed, which is commonly flatter, suggesting a gradual change in depositional conditions. To explain their origins, we appeal to other evidence. DSDP cores from the Late Pleistocene contain the rigid aragonite cements formed by restricted conditions, but not the lower and middle PP. Furthermore, mid-PP sedimentary δ18O values are similar to global ocean δ18O for that time, not enhanced as expected from excess evaporation. These data suggest that there was a greater exchange of Red Sea waters with the Indian Ocean during the mid-PP. That exchange may have allowed waters densified by evaporation in shallow regions of the northern Red Sea to flow south vigorously (the mounds would then be contourites). Alternatively, as the Pliocene seabed was shallower, wind-driven eddies may have affected more of the water column. Overall, the results indicate for the first time that deep circulation was stronger in the earlier PP compared with the present day. That circulation needs to be considered when evaluating organism dispersions across the Red Sea, regional climate, and influence of Red Sea Outflow Water on Indian Ocean Intermediate Water. Plain language summaryDuring the Pliocene (about 3–5 million years ago), temperatures on Earth's surface were similar to those predicted in some models of Earth's future climate. This has led to researchers becoming interested in studying geological evidence of conditions during the Pliocene as a clue to Earth's future. In detail, however, the comparison may not be exact. In our article, we describe mounds of sediment that were originally formed on the bed of the Red Sea in the Pliocene. They appear similar to snow drifts, and like snow drifts, may have formed under steady currents. Sea water becomes dense with evaporation (making it more saline). At the present day, dense water created in the Gulf of Suez cascades down into the deep Red Sea as a slow current. Perhaps such movements were more vigorous in the Pliocene, helped by the water escaping through a deeper barrier in the south into the Indian Ocean? Alternatively, these mounds were formed by giant eddies, such as those on the surface of the modern Red Sea, which are moved by strong winds. In either case, the deep waters of the Pliocene Pliocene Red Sea were more vigorously moving compared to the quiescent deep waters of the modern-day.

A High-Precision Sub-Grid Parameterization Scheme for Clear-Sky Direct Solar Radiation in Complex Terrain—Part I: A High-Precision Fast Terrain Occlusion Algorithm

In atmospheric modeling, sub-grid parameterization is an important method for studying the topographic effects of solar radiation using high-resolution Digital Elevation Model (DEM) data. For reducing the amount of computation, some approximate methods that can lead to errors are used in existing sub-grid parameterization schemes for clear-sky direct solar radiation (SPS-CSDSR). The lack of a high-precision fast terrain occlusion algorithm (HPFTOA) remains one of the biggest constraints in this field. This study proposed an HPFTOA. It mainly uses two kinds of acceleration algorithms. One method is to use a dynamic, lossless, and fast occlusion search radius. Another way is to use the rectangular grid for calculations within the accuracy of DEM data to avoid coordinate projection conversions. The test results indicate that the HPFTOA can carry out large-scale computation based on DEM data with a resolution of 90 m. Because it rarely uses approximation algorithms and considers the curvature of the Earth, SPS-CSDSR can achieve unprecedented precision. The HPFTOA can also be used in the fields of mountain solar energy assessment, remote sensing, and telemetry, including terrain-obscuring the probe. As computer performance improves and algorithms and execution code are optimized, the application prospects will be very broad.

PRIME‐SH: A Data‐Driven Probabilistic Model of Earth's Magnetosheath

AbstractA data‐driven model of Earth's magnetosheath is developed by training a recurrent neural network (RNN) with probabilistic outputs to reproduce Magnetospheric MultiScale (MMS) measurements of the magnetosheath plasma and magnetic field using measurements from the Wind spacecraft upstream of Earth at the first Earth‐Sun Lagrange point (L1). This model, called Probabilistic Regressor for Input to the Magnetosphere Estimation‐magnetosheath (PRIME‐SH) in reference to its progenitor algorithm PRIME, is shown to predict spacecraft observations of magnetosheath conditions accurately in a statistical sense with a continuous rank probability score of 0.227σ (dimensionless standard deviation units). PRIME‐SH is shown to be more accurate than many current analytical models of the magnetosheath. Furthermore, PRIME‐SH is shown to reproduce physics not explicitly enforced during training, such as field line draping, the dayside plasma depletion layer, the magnetosheath flow stagnation point, and the Rankine‐Hugoniot MHD shock jump conditions. PRIME‐SH has the additional benefits of being computationally inexpensive relative to global MHD simulations, being capable of reproducing difficult‐to‐model physics such as temperature anisotropy, and being capable of reliably estimating its own uncertainty to within 3.5%.

Thermal Earth model for the conterminous United States using an interpolative physics-informed graph neural network

This study presents a data-driven spatial interpolation algorithm based on physics-informed graph neural networks used to develop a thermal Earth model for the conterminous United States. The model was trained to approximately satisfy Fourier’s Law of conductive heat transfer by simultaneously predicting subsurface temperature, surface heat flow, and rock thermal conductivity. In addition to bottomhole temperature measurements, we incorporated other spatial and physical quantities as model inputs, such as depth, geographic coordinates, elevation, sediment thickness, magnetic anomaly, gravity anomaly, gamma-ray flux of radioactive elements, seismicity, electrical conductivity, and proximity to faults and volcanoes. With a spatial resolution of 18km2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$18 \\ km^2$$\\end{document} per grid cell, we predicted heat flow at surface as well as temperature and rock thermal conductivity across depths of 0-7km\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$0-7 \\ km$$\\end{document} at an interval of 1km\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$1 \\ km$$\\end{document}. Our model showed temperature, surface heat flow and thermal conductivity mean absolute errors of 6.4∘C\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$6.4^\\circ C$$\\end{document}, 6.9mW/m2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$6.9 \\ mW/m^2$$\\end{document} and 0.04W/m-K\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$0.04 \\ W/m-K$$\\end{document}, respectively. This thorough modeling of the Earth’s thermal processes is crucial to understanding subsurface phenomena and exploiting natural underground resources. Our thermal Earth model is available as web application at https://stm.stanford.edu, feature layers on ArcGIS at https://arcg.is/nLzzT0, and tabulated data on the Geothermal Data Repository at https://gdr.openei.org/submissions/1592.

From the Spherical Earth Model to the Globe: The Effectiveness of a Planetary Model-Building Intervention

From the Spherical Earth Model to the Globe: The Effectiveness of a Planetary Model-Building Intervention

Meshfree modelling of magnetotelluric and controlled-source electromagnetic data for conductive earth models with complex geometries

Meshfree modelling of magnetotelluric and controlled-source electromagnetic data for conductive earth models with complex geometries

Towards a new standard for seismic moment tensor inversion containing 3-D earth structure uncertainty

SUMMARY Moment tensor (MT) inversion is a classical geophysical inverse problem that infers a force-equivalent model of a seismic source from seismological observations. Like other inverse problems, the accuracy of the inversion depends on the reliability of the forward problem simulating waveforms from the source location through an Earth structural model. Apart from errors in data, the error in forward waveform simulation, also known as theory error, is a significant source of error contributing to the misfit function between the predicted and observed waveforms. Here, we set up numerical experiments to comprehensively probe the sensitivity of the linearized MT inversion to 3-D regional earth model errors, a known predominant factor of the theory error. Using the Monte Carlo method, we estimate the empirical structural covariance matrices to characterize the waveform mismatch due to the imperfect knowledge of Earth's structure. First, although the inversion accuracy deteriorates with increasing model errors, incorporating the structural covariance matrices into the misfit function improves the accuracy of inversion results for all theorized error distributions. Secondly, we propose a slightly modified form of the structural covariance matrix, which further enhances the inversion outcome. Lastly, as the true structural errors are likely spatially correlated, we highlight the importance of adequately treating the correlation into the MT inversion because of its significant impact on inversion. Overall, as a preliminary effort in quantifying 3-D structural errors on MT inversion, this study proves the computational feasibility by means of numerical experiments and will hopefully provide a way forward for future work on this topic.

Seismic modeling using pseudo-impedance derived from physical models

Building accurate numerical models of the earth for simulation is both challenging and time consuming. The subsurface complexity must be captured in sufficient detail for the results to be of practical use for geophysical analysis. Physical models provide an efficient way to add realistic detail to numerical models. Physical models, for example, are scaled simulations of large-scale geologic systems that provide great insight into understanding the interplay between sedimentation, salt tectonics, and faulting. The integration of geophysical workflows with sandbox models has been the subject of previous research, primarily to create synthetic data to compare with observed field data. However, capturing the detail from the physical model is challenging. A novel workflow is presented that uses pixel intensity values from photographic images of the analogue model to create a pseudo-impedance perturbation, mimicking the acoustic impedance contrasts that are encountered in the earth. This solution enables the detailed sedimentary architecture from the sandbox to be transferred to the digital earth model. Acoustic wave equation forward modeling and reverse time migration imaging have been applied using the generated pseudo-impedance model for a Gulf of Mexico-style salt canopy model. The images produced show clear definition of the salt/sediment interface and capture the sedimentary stratigraphy and faulting observed in the original sandbox model.

Generating extremely low-dimensional representation of subsurface earth models using vector quantization and deep Autoencoder

Geological model compression is crucial for making large and complex models more manageable. By reducing the size of these models, compression techniques enable efficient storage, enhance computational efficiency, making it feasible to perform complex simulations and analyses in a shorter time. This is particularly important in applications such as reservoir management, groundwater hydrology, and geological carbon storage, where large geomodels with millions of grid cells are common. This study presents a comprehensive overview of previous work on geomodel compression and introduces several autoencoder-based deep-learning architectures for low-dimensional representation of modified Brugge-field geomodels. The compression and reconstruction efficiencies of autoencoders (AE), variational autoencoders (VAE), vector-quantized variational autoencoders (VQ-VAE), and vector-quantized variational autoencoders 2 (VQ-VAE2) were tested and compared to the traditional singular value decomposition (SVD) method. Results show that the deep-learning-based approaches significantly outperform SVD, achieving higher compression ratios while maintaining or even exceeding the reconstruction quality. Notably, VQ-VAE2 achieves the highest compression ratio of 667:1 with a structural similarity index metric (SSIM) of 0.92, far surpassing the 10:1 compression ratio of SVD with a SSIM of 0.9. The result of this work shows that, unlike traditional approaches, which often rely on linear transformations and can struggle to capture complex, non-linear relationships within geological data, VQ-VAE's use of vector quantization helps in preserving high-resolution details and enhances the model's ability to generalize across varying geological complexities.