SEARCH · Search NASA
Results for “geologic structure”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Composite Analysis for Low-Level Waste Disposal in the Hanford Site Central Plateau (FY 2020) (Rev. 2)
This document provides an updated Hanford Site composite analysis (CA). A Hanford Site CA was prepared and issued in 1998 (PNNL-11800) with an addendum provided in 2002 (PNNL-11800 Addendum 1). The CA was approved in 2002 (Frei, 2002) and has been maintained (DOE/RL-2000-29 and subsequent revisions) to support low-level waste disposal performance assessments (PAs) and disposal authorizations for facilities at the Hanford Site, including the following: Continued operation of the Environmental Restoration Disposal Facility (ERDF) and the 200 East and 200 West Low-Level Burial Grounds; Construction of the Integrated Disposal Facility (IDF); Forthcoming closure of tank residual waste systems such as Waste Management Area (WMA) C. The CA maintenance program resulted in a determination in 2015 (DOE/RL-2015-66) that the Hanford Site CA needed an update for the following reasons: While the initial Hanford Site CA has been maintained since 2001, the accumulation of basis changes reported in the annual summary reports over the succeeding 14 years merit evaluation in an updated analysis; The U.S. Department of Energy (DOE) Headquarters requested in a memorandum in 2015 (Gilbertson and Marcinowski, 2015) that “as soon as the relevant PAs are complete, the CA will be revised to account for all of the new information.” This updated Hanford Site CA accounts for the following new information: 1. Inclusion of a detailed Hanford Site baseline disposition that projects remedial activities through site closure. There have been significant changes through decision making in the Comprehensive Environmental Response, Compensation, and Liability Act of 1980 (CERCLA) process that were not available when the original CA was produced; 2. Inclusion of an updated inventory basis, new modeling capabilities, and new decisions reached in the associated record of decision (ROD) that was provided by issuance of DOE/EIS-0391 in fiscal year (FY) 2013. Development of a Hanford Site groundwater model from the baseline provided in a technical transfer of models for the Final Tank Closure and Waste Management Environmental Impact Statement (DOE/EIS-0391) commenced in FY 2014 and resulted in the plateau-to-river (P2R) groundwater flow model that is used in this updated Hanford Site CA; 3. Accounting for pump and treat systems, which were not evaluated in the initial CA. Pump and treat systems have had significant impact on groundwater flow system behavior, contaminant transport, and contaminant removal from Hanford Site groundwater. This process is accounted for in the use of the P2R model for this updated Hanford Site CA; 4. Inclusion of water-level data collected since the initial CA was completed. Data collected as the unconfined aquifer water levels continue to recede since the cessation of large liquid discharges in the late 1990s have led to marked improvement in understanding of the flow system for future conditions, particularly regarding northward flow potential in the critical Gable Gap area. The calibration in the latest version of the P2R model accounts for this information; 5. Inclusion of results from the revised ERDF PA completed in 2013 (WCH-5209) that account for updated inventory and expansion of the ERDF facility to about twice the size that was evaluated in the original CA; 6. Use of updated geoframeworks to provide the structural basis for numerical fate and transport models in the groundwater pathway of this updated CA. The geologic basis for groundwater models has continued to improve with additional data collection and interpretation with the creation and maintenance of the Hanford South Geoframework and the Central Plateau Vadose Zone Geoframework tools; 7. Incorporation of updated tank residual inventory estimates. Tank residual inventory estimates have improved with the incorporation of tank retrieval inventory data for those tanks that have completed retrieval. The CA inventory data package includes this updated information; 8. Incorporation of WMA C PA results. Two additional PAs for tank farm closure decisions are in preparation during the period required to prepare an updated Hanford Site CA: WMA C (FY 2016) and WMA A-AX (in preparation). The WMA C PA results are incorporated into the updated Hanford Site CA, and its grouted residuals model is used as the basis for a release model to account for the other tank farm systems modeling in this updated CA; 9. An update to the IDF PA was submitted in FY 2017 and has been reviewed and approved. The results of the IDF PA are incorporated into this updated CA; 10. Updated risk assessment scenarios. The risk assessment scenarios currently in use for Hanford Site CERCLA and Resource Conservation and Recovery Act of 1976 (RCRA) analyses differ from those evaluated in the initial Hanford Site CA. The representative person exposure scenario evaluated in this updated Hanford Site CA is consistent with recent PAs and CERCLA and RCRA analyses as the Hanford Site. This updated Hanford Site CA provides the following: A comparison of the updated Hanford Site CA all-pathway dose results with the performance measures during the compliance period, which is assumed to begin with site closure in calendar 2070 with the last scheduled disposal action and continue for 1,000 years postclosure (to calendar year 3069); A comparison of Hanford Site CA all-pathway dose results with the performance measures during the postcompliance period to address potential peaks beyond the compliance period. This is accomplished by evaluating dose in the period for 9,000 years following the compliance period (i.e., from calendar year 3070 to 12070).
Optimizing Alabama’s CO 2 Storage in Shelby County (Project OASIS): Task 6.0 Deliverable – CarbonSAFE Phase III Readiness
The OASIS CarbonSAFE Phase II Project OASIS (Optimizing Alabama’s CO 2 Storage in Shelby County) is a geologic and reservoir characterization study designed to evaluate deep saline formations for potential long-term carbon dioxide (CO 2 ) storage in central Alabama near the National Carbon Capture Center (NCCC) and Alabama Power’s Plant Gaston. The project centers on understanding the potential of the Cambro-Ordovician Knox Group and underlying strata such as the Conasauga and Rome Formations for geologic storage of CO 2 . These formations were investigated as part of SECARB-USA (DE-FE0031830) and Project OASIS through the drilling of two stratigraphic test wells to obtain electronic well logs, core, and sidewall core plugs. These data provide direct measurements of porosity, permeability, and lithologic variability critical for reservoir characterization. Complementing the well program, a limited 2D seismic survey was conducted to help select the site for Westover #2 as part of SECARB-USA (DE-FE0031830), and a more regional Seismic Exchange (SEI) seismic survey was licensed and interpreted to define structural and stratigraphic frameworks in a new Static Earth Model (SEM), map reservoir continuity, and to identify potential sealing intervals. Integrated with geologic and reservoir modeling, these datasets form the basis for evaluating storage capacity, injectivity, and containment. While this document endeavors to provide readers with a high-level overview of Project OASIS activities and its suitability for subsequent CarbonSAFE Phases, such as a Phase III effort. Other project deliverables and milestones will provide more details on individual subjects.
Optimizing Alabama’s CO 2 Storage in Shelby County (Project OASIS): Task 4.0 Deliverable – Geologic Analysis Report
Project OASIS (Optimizing Alabama’s CO 2 Storage in Shelby County) is a geologic and reservoir characterization study designed to evaluate deep saline formations for potential long-term carbon dioxide (CO 2 ) storage in central Alabama near the National Carbon Capture Center (NCCC) and Alabama Power’s Plant Gaston. The project centers on the Cambro-Ordovician Knox Group and underlying strata such as the Conasauga and Rome Formations, which were investigated through the drilling of two stratigraphic test wells to obtain electronic well logs, core, and sidewall core plugs. These data provide direct measurements of porosity, permeability, and lithologic variability critical for reservoir characterization. Complementing the well program, a limited 2D seismic survey was acquired to help select the site for Westover #2, and a more regional Seismic Exchange (SEI) seismic survey was licensed and interpreted to define structural and stratigraphic frameworks in a new Static Earth Model (SEM), map reservoir continuity, and identify potential sealing intervals. Integrated with geologic and reservoir modeling, these datasets form the basis for evaluating storage capacity, injectivity, and containment.
RANGERS: Modeling Report on Integrity and Performance Assessment of Engineered Barrier Systems in a Salt Repository for HLW/SNF
The Engineered Barrier System (EBS) plays an important role in ensuring the long-term safety and containment of high-level waste (HLW) and spent nuclear fuel (SNF) in deep geological repositories in salt formation. As part of a multi-barrier system, the EBS works alongside the natural barrier, which is the salt formation itself and the technical barrier comprising the disposal casks. The primary function of the EBS is to maintain containment during a defined period until the backfill used in the repository made of crushed salt, develops its sealing capacity through compaction. Over the time, the backfill eventually compacts to a state of low porosity and permeability, acting as a long-term seal. However, until this process is complete, the EBS must retain its structural and functional integrity. Regulatory guidelines in Germany currently require the EBS to remain effective for up to next ice age, that is expected in 50,000 years. The significant hydro-geological and topographic changes expected during an ice age could make it impossible to accurately predict the hydro-chemical conditions within the repository system at that time. In response to these challenges, BGE TECHNOLOGY GmbH (BGE TEC) and Sandia National Laboratories (SNL) have jointly developed a comprehensive methodology for the design and safety assessment of engineered barrier systems within the scope of the RANGERS project. This methodology is tailored for repositories in salt formations. The developed methodology provides a structured approach for designing and assessing the performance of the EBS in salt-based repositories. It begins with defining a sealing concept based on the geological characteristics of the selected site and the overall repository design. The entire repository system, comprising the geological site, repository infrastructure, and EBS, is then subjected to a Features, Events, and Processes (FEP) analysis, focusing solely on those FEPs that affect the EBS. The derived FEPs help identify the loads and stresses acting on the EBS, which serve as the foundation for conducting an integrity assessment. This analysis helps predict the EBS’s evolution and performance over the regulatory time frame, feeding into integrated performance assessment simulations.
Multimodal Mass Spectrometry Imaging (MSI) of Archean and Jurassic Geologic Samples
Current geochemical analytical techniques such as gas chromatography-mass spectrometry (GC-MS and $G\acute{C}GC-MS$) can effectively determine the composition and structure of the organic biosignatures within a sample but cannot resolve the spatial distribution of organic biosignatures within rocks or sediments. The spatial distribution can be used to determine the character of the organic biosignatures [1], which can be indigenous (deposited with host rocks), non-indigenous (incorporated after deposition via fluid migration), or contaminant. Determining the character of organic biosignatures is critical for ancient samples [2] and will be critical if organics are observed in samples returned from Mars [3].Fs-LDPI-MS has the ability to map organic compounds across the surface of samples at 2 - 10 μm lateral resolution. Furthermore, fs-LDPI-MS can be used to carry out multiple analyses in the same location for micron-scale analysis of previously buried material [4]. We used fs-LDPI-MS to examine a ~164 million year old organic-rich mudstone from SW England (14.2 wt. % total organic carbon) [5] and determined the spatial distribution of likely indigenous components buried below potentially contaminated surface layers.We then used ToF-SIMS to re-analyze the mudstone to directly compare ToF-SIMS to fs-LDPI-MS datasets, and to analyze a series of ~2.7 billion year old geologic samples from Timmins, ON, CA [6]. Ref. [6] previously analyzed the Archean samples and observed a series of archaeal biomarkers, as well as hopanes and steranes. Using ToF-SIMS analysis, we were able to detect hopanes and steranes in the mudstone samples, but not in the Archean samples indicating that the hopanes and steranes are most likely contaminants in the Archean samples. We also did not observe the archaeal biomarkers; however, this may be from the lack of molecular ion preservation during ToF-SIMS.
Uranium (VI) Adsorbate Structures on Portlandite [Ca(OH)2] Type Surfaces Determined by Computational Modelling and X-ray Absorption Spectroscopy
Portlandite [Ca(OH) 2 ] is a potentially dominant solid phase in the high pH fluids expected within the cementitious engineered barriers of Geological Disposal Facilities (GDF). This study combined X-ray Absorption Spectroscopy with computational modelling in order to provide atomic-scale data which improves our understanding of how a critically important radionuclide (U) will be adsorbed onto this phase under conditions relevant to a GDF environment. Such data are fundamental for predicting radionuclide mass transfer. Surface coordination chemistry and speciation of uranium with portlandite [Ca(OH) 2 ] under alkaline groundwater conditions (ca. pH 12) were determined by both in situ and ex situ grazing incidence extended X-ray absorption fine structure analysis (EXAFS) and by computational modelling at the atomic level. Free energies of sorption of aqueous uranyl hydroxides, [UO 2 (OH) n ] 2–n (n = 0–5) with the (001), (100) and (203) or (101) surfaces of portlandite are predicted from the potential of mean force using classical molecular umbrella sampling simulation methods and the structural interactions are further explored using fully periodic density functional theory computations. Although uranyl is predicted to only weakly adsorb to the (001) and (100) clean surfaces, there should be significantly stronger interactions with the (203/101) surface or at hydroxyl vacancies, both prevalent under groundwater conditions. The uranyl surface complex is typically found to include four equatorially coordinated hydroxyl ligands, forming an inner-sphere sorbate by direct interaction of a uranyl oxygen with surface calcium ions in both the (001) and (203/101) cases. In contrast, on the (100) surface, uranyl is sorbed with its axis more parallel to the surface plane. The EXAFS data are largely consistent with a surface structural layer or film similar to calcium uranate, but also show distinct uranyl characteristics, with the uranyl ion exhibiting the classic dioxygenyl oxygens at 1.8 Å and between four and five equatorial oxygen atoms at distances between 2.28 and 2.35 Å from the central U absorber. These experimental data are wholly consistent with the adsorbate configuration predicted by the computational models. These findings suggest that, under the strongly alkaline conditions of a cementitious backfill engineered barrier, there would be significant uptake of uranyl by portlandite to inhibit the mobility of U(VI) from the near field of a geological disposal facility.
Data for The Value of Reversible Carbon Storage in a Zero-Emissions World
Atmospheric carbon dioxide removal (CDR) is required to stabilize global temperature. CDR can be achieved via ecosystem-based approaches that are cost-effective but reversible (e.g., soil and forest management) or by more durable but expensive approaches (e.g., direct air capture coupled with geologic storage). Here, we examine trade-offs between these approaches, focusing on timing, climate impacts, and cost. We simulated reversible carbon accrual for a range of CDR contract structures using a general minimalist model of ecosystem carbon cycling, and parameterized it to simulate US agricultural soil management─specifically cover cropping─as a case study. We then quantified the resulting impact on atmospheric carbon and global temperature using a climate model emulator. We find that maintaining a patchwork of reversible CDR projects by replacing lapsed projects with new projects can reduce warming by 22–195 μ°C in 2100 and that the magnitude of this cooling effect depends on how effectively the patchwork is maintained. Long-term maintenance of reversible CDR projects requires institutional stability that cannot be guaranteed over multiple decades. Consequently, effective CDR ultimately requires replacing reversible projects with durable projects. To address this problem, we modeled the cost of replacing reversible agricultural soil CDR with geologic CDR. We found that using reversible CDR as a bridge to durable CDR is potentially more cost-effective as a global cooling strategy (0.20–0.81 billion USD per μ°C avoided) than perpetual maintenance of reversible CDR (0.32–1.31 billion USD per μ°C avoided) or an immediate transition to durable CDR (1.37–2.19 billion USD per μ°C avoided). However, we emphasize that institutional commitments to maintain reversible CDR projects cannot be guaranteed. Reliance on reversible CDR as a bridge to durable CDR therefore carries an unknown amount of risk and will only function if efforts to maintain reversible CDR are robust.
Corrosion integrity of oil cement modified by environment responsive microspheres for CO2 geologic sequestration wells
In order to ensure the safe and effectiveness of CO{sub 2} geological storage, an environmental responsive polymer microspheres (ERPM) was prepared to improve the corrosion resistance of cement stone. The structure and environmental response characteristics of ERPM were characterized, and then the anti-corrosion performance and anti-corrosion mechanism for ERPM were discussed. The results shown that ERPM had good environmental response characteristics and temperature resistance. ERPM could effectively suppressed the corrosion rate and the damage degree of mechanical properties of cement stone. The corrosion integrity of cement stone was improved due to the acidic response characteristic and form polymer films of ERPM, which effectively shield the direct contact between corrosion medium and hydration products.
Probabilistic Evaluation of Geoscientific Hypotheses with Geophysical Data: Application to Electrical Resistivity Imaging of a Fractured Bedrock Zone
As climate changes and populations grow, groundwater sustainability is becoming increasingly important. Groundwater models, based on a conceptual understanding of the subsurface structure, are crucial tools for making sustainable management decisions. Conceptual models of the subsurface are based on knowledge of geological processes, and, frequently, observations from geophysical data. A frequent problem in groundwater model development occurs when multiple geological phenomena could explain a single subsurface observation. Uncertainty in geophysical data makes it even more difficult to discern which explanations are consistent with the geophysics. Here, we present a framework for testing geological when a geological feature is observed in geophysical data, but its physical characteristics are uncertain. The framework builds on Popper-Bayes methods developed in previous work, and is applied to study a fractured bedrock zone in a mountainous watershed in southwest Colorado. First, we outline six hypotheses based on the geological history of the watershed. Then, using the proposed Popper-Bayes approach, we demonstrate that three of the six hypotheses are inconsistent with measured electrical resistivity data, even after accounting for uncertainty. Finally, we discuss the importance of the prior model, and how this framework for handling geophysical uncertainty can be applied in other settings.
MST e-News Spring 2022
For Jeremy Mitchell, a scientific career grounded in a passion for geology has morphed into a decades-long fascination with plutonium. Just as a geologist investigates a rock’s structure and its history to unlock its secrets, so too does Mitchell gather facts and weigh evidence to crack the code of a material not readily accessible to most scientists—and at the core of the Lab’s national security mission.
Modeling supercritical CO 2 flow and mineralization in reactive host rocks with PFLOTRAN v7.0
Understanding the flow and reactivity of CO 2 injected into geological reservoirs is important for many subsurface applications including secure geologic carbon storage (GCS), critical mineral extraction, enhanced geothermal systems (EGS), and enhanced oil recovery (EOR). Traditionally, subsurface CO 2 injection for GCS applications has focused on geologic formations with favorable subsurface configurations for CO 2 migration and trapping through non-reactive mechanisms such as structural, solubility, and petrophysical trapping. Recently, CO 2 -reactive rocks such as mafic and ultramafic basalts have been investigated for their potential to react with injected CO 2 in situ to simultaneously dissolve host rock minerals and mineralize CO 2 as carbonates. Engineering rapid CO 2 mineralization in the subsurface is attractive because of the increased density of stored CO 2 , the additional safety factors associated with solidification, and the potential to extract valuable critical minerals. Here we present recent developments in the parallel flow and reactive transport simulator PFLOTRAN to model coupled CO 2 -brine flow and reactive transport for a wide range of injection and production applications involving reactive CO 2 -brine systems. These developments are based on the well established and trusted CO 2 flow capabilities in the STOMP-CO 2 simulator. New capabilities added to PFLOTRAN include new CO 2 -brine equations of state with optional thermal coupling, several new constitutive relationships like capillary pressure smoothing and scanning path hysteresis, a fully implicit well model, and native linkage with PFLOTRAN's well-established reactive transport libraries. A series of benchmarks between PFLOTRAN and STOMP-CO 2 verify the newly developed CO 2 -brine flow capabilities. Demonstrations of coupled CO 2 -brine flow modeling and reactive transport show how CO 2 mineralization can be engineered in reactive host rocks. Finally, an example use case involving copper leaching by CO 2 and critical mineral extraction is presented to showcase the strengths of this new implementation. Several limitations still remain, including limited availability of field data to parameterize models. Future work should constrain the evolution of mineral surface area during mineralization and the temperature and/or pH dependence of geochemical reactions for specific systems of interest.
LiDAR-Aided Interior Orientation Parameters Refinement Strategy for Consumer-Grade Cameras Onboard UAV Remote Sensing Systems
Unmanned aerial vehicles (UAVs) are quickly emerging as a popular platform for 3D reconstruction/modeling in various applications such as precision agriculture, coastal monitoring, and emergency management. For such applications, LiDAR and frame cameras are the two most commonly used sensors for 3D mapping of the object space. For example, point clouds for the area of interest can be directly derived from LiDAR sensors onboard UAVs equipped with integrated global navigation satellite systems and inertial navigation systems (GNSS/INS). Imagery-based mapping, on the other hand, is considered to be a cost-effective and practical option and is often conducted by generating point clouds and orthophotos using structure from motion (SfM) techniques. Mapping with photogrammetric approaches requires accurate camera interior orientation parameters (IOPs), especially when direct georeferencing is utilized. Most state-of-the-art approaches for determining/refining camera IOPs depend on ground control points (GCPs). However, establishing GCPs is expensive and labor-intensive, and more importantly, the distribution and number of GCPs are usually less than optimal to provide adequate control for determining and/or refining camera IOPs. Moreover, consumer-grade cameras with unstable IOPs have been widely used for mapping applications. Therefore, in such scenarios, where frequent camera calibration or IOP refinement is required, GCP-based approaches are impractical. To eliminate the need for GCPs, this study uses LiDAR data as a reference surface to perform in situ refinement of camera IOPs. The proposed refinement strategy is conducted in three main steps. An image-based sparse point cloud is first generated via a GNSS/INS-assisted SfM strategy. Then, LiDAR points corresponding to the resultant image-based sparse point cloud are identified through an iterative plane fitting approach and are referred to as LiDAR control points (LCPs). Finally, IOPs of the utilized camera are refined through a GNSS/INS-assisted bundle adjustment procedure using LCPs. Seven datasets over two study sites with a variety of geomorphic features are used to evaluate the performance of the developed strategy. The results illustrate the ability of the proposed approach to achieve an object space absolute accuracy of 3–5 cm (i.e., 5–10 times the ground sampling distance) at a 41 m flying height.
Effects of Supercritical CO 2 Injection on the Shale Pore Structures and Mass Transport Rates
Characterizing the pore structures and transport properties of low-permeability shales is critical for evaluating these formations as potential seals or storage sites for geological CO 2 sequestration. Here, we use low-pressure gas adsorption in conjunction with nuclear magnetic resonance (NMR) to characterize the pore-size distribution of shales before and after injection of supercritical CO 2 . Nitrogen gas was used as the detecting phase for the adsorption experiments and pentane liquid was used for the NMR experiments. We also performed time-resolved NMR and gravimetric microbalance measurements to observe mass transport during desaturation. We use these data to estimate the self-diffusion coefficient of pentane and changes in the saturation state of the pore network. We analyzed samples with a range of compositions from the Wolfcamp shale before and after exposure to supercritical CO 2 for 3 days. Integrating the gas adsorption and NMR data shows how supercritical CO 2 injection alters the pore-size distribution for pore sizes <1 nm to 1 mm. Finally, our results provide insights on how the pore structure and mass transport properties of different shale lithologies may evolve during storage of supercritical CO 2 .
Using Complementary Methods of Synchrotron Radiation Powder Diffraction and Pair Distribution Function to Refine Crystal Structures with High Quality Parameters—A Review
Determination of the atomic-scale structures of certain fine-grained minerals using single-crystal X-ray diffraction (XRD) has been challenging because they commonly occur as submicron and nanocrystals in the geological environment. Synchrotron powder diffraction and scattering techniques are useful complementary methods for studying this type of minerals. In this review, we discussed three example studies investigated by combined methods of synchrotron radiation XRD and pair distribution function (PDF) techniques: (1) low-temperature cristobalite; (2) kaolinite; and (3) vernadite. Powder XRD is useful to determine the average structure including unit-cell parameters, fractional atomic coordinates, occupancies and isotropic atomic displacement parameters. X-ray/Neutron PDF methods are sensitive to study the local structure with anisotropic atomic displacement parameters (ADP). The results and case studies suggest that the crystal structure and high-quality ADP values can be obtained using the combined methods. The method can be useful to characterize crystals and minerals that are not suitable for single-crystal XRD.
Joint Inversion of Surface Electrical Resistivity Tomography and Seismic Refraction Data between the 200 Areas
Geologic stratigraphy on the Hanford Site influences groundwater and contaminant migration through the aquifer system and the vadose zone. The current geologic framework model (GFM) relies heavily on a sparse distribution of borehole data in some locations to map geologic contacts and hydrologic properties in the subsurface. Non-invasive geophysical methods such as electrical resistivity tomography (ERT), transient electromagnetic surveying, and seismic imaging are being used at Hanford to map subsurface structure in areas with limited well observations. This is to develop and mature the capability of geophysical methods to aid in GFM refinement, to identify regions of subsurface complexity, and for optimal well siting. A joint inversion of co-located seismic refraction and ERT data was carried out for data collected on a ~2.3-km profile between the 200 Areas on the Hanford Site. While ERT and seismic refraction images have sensitivity to overlapping physical properties (porosity, moisture content, lithology), the resolution and physics used to acquire each of these datasets are different and therefore information can be different or mutually complementary. Performing a joint inversion provides a reasonable option for a coherent, coupled interpretation for mutually complementary datasets. Between the 200 Areas, there are few boreholes to interpret the geologic framework model, and these data sets were obtained to provide a first line of evidence toward identifying stratigraphic structure. The seismic refraction and ERT data were independently inverted during fiscal year 2022 and broadly showed a two-layer structure with a trough-like feature that is ~1 km wide and upwards of 150 m deep. The depth of the trough feature was greater in the ERT image compared to the seismic image, which indicated a maximum depth of approximately 110 m. The objective of the joint inversion described in this report was to invert the seismic refraction and ERT data together while constraining the ERT image to be structurally similar to the seismic refraction image. The approach was applied using the geophysical inverse modeling program E4D, which has the capability to invert first-arrival times from seismic refraction data and ERT resistances using a “cross-gradient” constraint. The application of cross-gradient constraints with different weights produces ERT models that show a high degree of similarity within the upper 100 m (above ~120 m elevation). None of the ERT models show an improved structural similarity to the seismic result; therefore, it is recommended that further attempts to jointly interpret these models focus on petrophysics and image resolution. Petrophysical measurements of core samples would improve knowledge of what drives the ERT response in this region and, along with downhole geophysical measurements, could be used to “ground truth” the surface-based geophysical results. Image resolution studies would provide insight into which regions of the inverted images are reliable and which regions are poorly constrained.
The Collaborative Seismic Earth Model: Generation 2
Geological interpretations, earthquake source inversions and ground motion modeling, among other applications, require models that jointly resolve crustal and mantle structure. With the second generation of the Collaborative Seismic Earth Model (CSEM2), we present a global multi-resolution tomographic Earth model that serves this purpose. The model evolves through successive regional- and global-scale refinements. While the first generation aggregated regional models, with this study, we ensure consistency between all individual submodels, resulting in a model that accurately explains wave propagation across scales. Recent regional tomographic models were incorporated, comprising continental-scale inversions for Asia and Africa, as well as regional inversions for the Western US, Central Andes, Iran, and Southeast Asia. Across all regional refinements, over 793,000 source-receiver pairs contributed. Moreover, the long-wavelength Earth model (LOWE) introduces large-scale structures outside of pre-existing local refinements. A full-waveform inversion for global anisotropic P-and S-wave speed structure over a total of 194 iterations with a minimum period of 50 s on a large data set of 1 hr of waveform data from 2,423 earthquakes and over 6 million source-receiver pairs ensures that regional updates in the crust and uppermost mantle translate into updates of deeper, global-scale structure. To test the performance of CSEM2, we evaluate waveform fits between observed and synthetic seismograms at 50 s for an independent data set on the global scale, and on the regional scale for lower periods. We accurately simulate waveforms within and across regional refinements, maintaining the original resolution of the submodels embedded in the global framework.
Block encoding of the three-dimensional heterogeneous Poisson equation with application to fracture flow
Quantum linear system (QLS) algorithms offer the potential to solve large-scale linear systems exponentially faster than classical methods. However, applying QLS algorithms to real-world problems remains challenging due to issues such as state preparation, data loading, and efficient information extraction. In this work, we study the feasibility of applying QLS algorithms to solve discretized three-dimensional (3D) heterogeneous Poisson equations, with specific examples relating to groundwater flow through geologic fracture networks. We explicitly construct a block encoding for the 3D heterogeneous Poisson matrix by leveraging the sparse local structure of the discretized operator. While classical solvers benefit from preconditioning, we show that block encoding the system matrix and preconditioner separately does not improve the effective condition number that dominates the QLS run-time. This differs from classical approaches where the preconditioner and the system matrix can often be implemented independently. Nevertheless, due to the structure of the problem in three dimensions, the quantum algorithm achieves a run-time of 𝑂(𝑁 2/3 polylog 𝑁 ⋅log (1/𝜖)), outperforming the best classical methods (with run times of 𝑂(𝑁log 𝑁 ⋅log (1/𝜖))) and offering exponential memory savings. These results highlight both the promise and limitations of QLS algorithms for practical scientific computing, and point to effective condition-number reduction as a key barrier in achieving quantum advantages.