Rock Valley Fault Zone Absolute and Relative Earthquake Relocation to Improve Geologic Models
Seismological Society Annual Meeting Presentation on 4/21/2021 discussing the Rock Valley earthquake relocations.
SEARCH · Engineering Papers
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.
Seismological Society Annual Meeting Presentation on 4/21/2021 discussing the Rock Valley earthquake relocations.
Understanding CO2-brine-rock interactions in sandstone layers of the Conasauga Group, Bartow County, Georgia
Earth’s biosphere is thought to exert a substantial influence on regolith evolution and chemical weathering rates. However, ecosystems are also highly efficient at retaining and recycling nutrients. Thus, when the ecological demand for rock-derived nutrients (e.g., P, Ca, K) exceeds the rates of regolith supply, ecological retention and recycling strategies can minimize nutrient limitations. To evaluate the balance between nutrient recycling and new nutrient input, we combined a plant model that drives growth according to foliar P levels with a weathering model that includes regolith rejuvenation via erosion and export via chemical weathering according to water flow, regolith thickness, mineral dissolution rates, secondary minerals, and nutrient storage in organic and mineral phases. We find that plant growth is strongly dependent on the total regolith nutrient inventory, resulting in a strong correlation between plant productivity and erosion. Increased water export or decreased regolith thickness diminish the total inventory of nutrient corresponding to lower rates of recycling and lower plant growth. In contrast, purported biogenic drivers of weathering, such as enhanced mineral dissolution, only support higher growth rates at high erosion rates. At erosion rates typical of the global land surface, more rapid mineral dissolution combined with enhanced formation of secondary minerals, depletes the inventory of mineral P, resulting in no benefit for plant growth. We also find that the increased chemical weathering export does not scale directly with plant growth. For example, accelerated mineral weathering does increase chemical weathering export but not potential plant growth. Conversely, thicker regolith is associated with a small increase in weathering export, but a large increase in potential plant growth. Collectively, when plant growth is coupled to regolith weathering our calculations suggest that plant productivity is not directly correlated with silicate weathering fluxes, and that biotic drivers of silicate weathering may only be effective at high erosion rates not typical at the Earth’s surface.
As global demand for rare earth elements (REEs) increases, maintaining the production and supply chain is critical. Technologies capable of being used in the field and in situ in the subsurface for rapid REE detection and quantification facilitates the efficient mining of known resources and exploration of new and unconventional resources. Laser-induced breakdown spectroscopy (LIBS) is a promising technique for rapid elemental analysis both in the laboratory and in the field. Multiple articles have been published evaluating LIBS for detection and quantification of REEs; however, REEs in their natural deposits have not been adequately studied. In this work, detection and quantification of two REEs, La and Nd, have been studied in both synthetic and natural mineral matrices at concentrations relevant to REE extraction. Measurements were performed on REE-containing rock and coal samples (natural and synthetic) utilizing different LIBS instruments and techniques, specifically a commercial benchtop instrument, a custom benchtop instrument (single- and double-pulse modes), and a custom LIBS probe currently being developed for in situ, subsurface, borehole wall detection and quantification of REEs. Plasma expansion, emission intensity, detection limits, and double-pulse signal enhancement were studied. The limits of detection (LOD) were found to be 10/14 ppm for La and 15/25 ppm for Nd in simulated coal/rock matrices in single-pulse mode. Signal enhancement of 3.5 to 6-fold was obtained with double-pulse mode as compared to single-pulse operation.
Overview of Rock Valley Earthquake Relocation for the 2021 Nuclear Explosion Monitoring Program Review
Pronounced voltage hysteresis in Li-excess cathode materials is commonly thought to be associated with oxygen redox. However, these materials often possess overlapping oxygen and transition-metal redox, whose contributions to hysteresis between charge and discharge are challenging to distinguish. In this work, a two-step aqueous redox titration is developed with the aid of mass spectrometry (MS) to quantify oxidized lattice oxygen and Mn 3+/4+ redox in a representative Li-excess cation-disordered rock salt—Li 1.2 Mn 0.4 Ti 0.4 O 2 (LMTO). Two MS-countable gas molecules evolve from two separate titrant-analyte reactions, thereby allowing Mn and O redox capacities to be decoupled. The decoupled O and Mn redox coulombic efficiencies are close to 100% for the LMTO cathode, indicating high charge-compensation reversibility. As incremental Mn and O redox capacities are quantitatively decoupled, each redox voltage hysteresis is further evaluated. Overall, LMTO voltage hysteresis arises not only from an intrinsic charge-discharge voltage mismatch related to O redox, but also from asymmetric Mn-redox overvoltages. The results reveal that O and Mn redox both contribute substantially to voltage hysteresis. In conclusion, this work further shows the potential of designing new analytical workflows to experimentally quantify key properties, even in a disordered material having complex local coordination environments.
Abstract Fluorination has been identified as a key element for enabling the stable cycling of earth‐abundant manganese‐based disordered rock salt (DRX) cathodes. However, fluorination in the DRX bulk remains a challenge for scalable solid‐state synthesis. In this study, a tailored reaction pathway is proposed to synthesize a highly fluorinated DRX. It is demonstrated for the first time that the unconventional precursors, Li 6 MnO 4 , MnF 2 , and TiO 2 , can avoid the formation of Mn‐based intermediates (such as Li 2 (Mn,Ti)O 3, LiMnO 2 , and Mn 3 O 4 ), which, once formed, persist until the synthesis temperature reaches close to or above that required for fluorine volatility. Therefore, this method can form a highly fluorinated DRX with a composition of Li 1.23 Mn 0.40 Ti 0.37 O 2−y F y ( y = 0.29–0.34) at a low temperature (800 °C) relative to that required for conventional DRX solid‐state reactions (≥900 °C). Li 1.23 Mn 0.40 Ti 0.37 O 2−y F y ( y = 0.29–0.34) delivers a specific capacity above 300 mAh g −1 and a specific energy of 980 Wh kg −1 at 30 °C. Detailed characterization reveals that this DRX phase reversibly utilizes Mn 2+/3+ redox in the low‐voltage region and Mn 3+/4+ redox in the middle‐voltage range, whereas reversible oxygen redox is observed at high potentials.
With the wide adoption of Li‐ion batteries, Ni‐rich cathode is considered as one of the most promising candidates of cathodes due to its high energy density and low cost. However, stability decreased with increasing Ni content in the Ni‐rich cathode. To solve this bottleneck, many strategies, such as coating, doping, surface modification, and special morphologies, have been developed. Herein, we introduce a groundbreaking approach for enhancing Ni‐rich cathode through an innovative acid etching process that promotes Mn shell self‐assembly, inducing a rock‐salt phase on the surface. This method not only simplifies the Ni‐rich cathode modification process, but also significantly improves the structural stability and electrochemical performance of Ni‐rich cathode. Our findings demonstrate that developed single‐crystal Ni‐rich cathode shows 3–34 % better stability compared to both commercial modified Ni‐rich cathode and unmodified counterparts. The unique Mn shell effectively mitigates reversible phase shifts during cycling, contributing to a remarkable enhancement in cycling stability. Additionally, this novel fabrication technique paves the way for cost‐effective production of high‐performance cathode materials, offering substantial benefits for lithium‐ion battery technology. And this study proves the potential of this method in advancing the design and development of durable, high‐capacity cathode materials for next‐generation batteries.
Abstract The application of lithium rich disordered rock salts (DRX) as cathode materials has greatly expanded the materials space for high energy density cathodes. These materials are able to consistently achieve capacities higher than via a complex percolation based intercalation mechanism. Most current DRX materials face significant capacity fade when cycled over an extended period. One of the factors responsible for this could be deleterious side effect of interface reactions between the electrode and electrolyte at high voltage. We aim to focus on one aspect of this interaction, i.e the effect of exposure to electrolyte on the surface and its effect on the electrochemical properties of (LMNOF) powder. LMNOF was systematically treated in an electrolyte solution for varying periods of time at an elevated temperature. Treated samples exhibit surface layer modification (removal of F rich surface layer), leading to a 10 % improvement in the capacity retention behavior of LMNOF. Surface and bulk based measurements indicate an increase in disorder and gradual removal of F and Li. All of these processes mimic an aging process similar to cycling. This a priori formed interface allows for increased stability with respect to retention of capacity and oxygen loss.
Measuring hydro-mechanical properties of natural fractures is a prerequisite for optimizing hydraulic stimulation design and well placement. We completed experiments to characterize shear on natural fractures in schist, amphibolite, and rhyolite specimens drilled from EGS Collab Project’s field sites at the Sanford Underground Research Facility (SURF) in South Dakota. A triaxial direct shear method and coupled x-ray imaging were used to perform hydroshearing and mechanical shearing at the site’s in-situ stress conditions. This produced simultaneous measurements of fracture and matrix strength, permeability, stress-dependent aperture, dilation, and friction strength. Our results identified that only a subset of the natural fractures was weak enough for hydroshearing. Generally, hydroshearing increases fracture permeability by a factor of 10 or more and the enhancement is retainable over time. However, the shear slip does not always result in permeability enhancement. High content of phyllosilicates was found to associate with exceptionally weak fractures that also exhibited poor or even negative enhancement after stimulation. Combining our measurements with site data, we can predict that most observable fractures at the two EGS Collab sites do not meet the criteria for hydroshearing before tensile opening. In some cases, the visible fractures are low permeability and as strong as the adjacent rock. To induce hydroshearing before tensile opening, injection must target known weak and favorably oriented fractures with confirmed pre-existing permeability.
The Mars Curiosity rover carries an instrument, ChemCam, designed to measure the composition of surface rocks and soil using laser-induced breakdown spectroscopy (LIBS). The measured spectra from this instrument must be analyzed to identify the component elements in the target sample, as well as their relative proportions. This process, which we call disaggregation, is complicated by so-called matrix effects, which describe nonlinear changes in the relative heights of emission lines as an unknown function of composition due to atomic interactions within the LIBS plasma. In this work, we explore the use of the plasma physics code ATOMIC, developed at Los Alamos National Laboratory, for the disaggregation task. ATOMIC has recently been used to model LIBS spectra and can robustly reproduce matrix effects from first principles. The ability of ATOMIC to predict LIBS spectra presents an exciting opportunity to perform disaggregation in a manner not yet tried in the LIBS community, namely via Bayesian model calibration. However, using it directly to solve our inverse problem is computationally intractable due to the large parameter space and the computation time required to produce a single output. Therefore, we also explore the use of emulators as a fast solution for this analysis. We discuss a proof of concept Gaussian process emulator for disaggregating two-element compounds of sodium and copper. The training and test datasets were simulated with ATOMIC using a Latin hypercube design. After testing the performance of the emulator, we successfully recover the composition of 25 test spectra with Bayesian model calibration.
Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.
X-ray computed tomography (CT) images of rock specimens often contain artifacts which must be corrected before scientific analyses are performed. Here, we present a new workflow of automated image processing to utilize poor-quality X-ray CT scan images. The workflow runs on the open-source image analysis software and efficiently separates desired features from low-contrast scanned images. The new workflow is a two-step technique using contrast enhancement and automated feature segmentation to generate noise-free binary images. The results of binary images using the proposed workflow and using a conventional thresholding technique are analyzed to show the quality of the proposed method. The paper also presents a workflow of estimating the structural geometries of features in two and three dimensions. The results of the structural feature analyses and computational time were compared between the open-source (ImageJ) and commercial image analysis software (Bruker Computed Tomography Analyzer). The commercial software was more computationally efficient, but the task-specific macros in open-source software enabled the user-desired automation in image processing and data extraction of desired structural features of comparable quality.
Neptunium poses challenges for radioactive wastes management due to its long half-life and ingrowth from americium decay, as well as its high mobility in water in its oxidized forms. Np(V)-bearing aqueous solutions were contacted with polished and crushed gneiss and dolerite from the Yeniseisky site, Russia, which is under consideration as a location for an underground repository, both in oxic and anoxic conditions at 20 and 90 °C. Neptunium sorption on crushed gneiss and dolerite in oxic conditions at 20 °C was quite low, but increased substantially in anoxic conditions. Autoradiography on samples treated at 20 °C revealed that sorbed Np was associated mostly with dark colored iron-bearing minerals, and solvent extraction indicated partial reduction of Np(V) to Np(IV). Sorption of Np from aqueous solutions contacted with dolerite and gneiss at 90 °C was in all cases higher than at the lower temperature. Autoradiography of polished rock samples treated in anoxic conditions at 90 °C revealed that Fe-bearing minerals had 3–4 orders of magnitude more Np sorbed on them than the quartz and feldspar grains. X-ray photoelectron spectroscopy indicated this sorbed Np was present as NpO 2 . This study supports the concept that the geologic conditions present at the planned repository site will retard Np migration.
This paper presents a coupled thermal-hydraulic-mechanical (THM) model for simulating failure processes in geological media. The coupling between stress/strain and heat transfer describing how thermal expansion affects stress is based on the theory of thermo-elasticity. Additionally, the coupling between stress/strain and fluid flow describing how pore-pressure affects stress is governed by Biot's consolidation theory. Stress/strain-dependent porosity/permeability/capillary-pressure defines the influence of stress on fluid flow and thermal fields. The THM model is implemented in the solver module that links the multiphase, multicomponent fluid flow simulator TOUGH2 with the geomechanical simulator RFPA (Rock Failure Process Analysis), which is based on the concept that heterogeneity leads to non-linearity and progressive failure behavior. Capability of TOUGH-RFPA simulator was verified by two benchmark problems related to hydraulic-mechanical and thermal-mechanical response in a hollow cylinder. Then the applicability of TOUGH-RFPA was demonstrated for modeling cooling-induced damage around a very deep geothermal well with mesoscopic heterogeneous properties, which showed complex damage evolution and fracturing during well cooling and fluid pressurization.
Here, we formulate a time-dependent damage theory for rocks subjected to mechanical deformation and solid dissolution. The constitutive description is inspired by the transition state theory, which states that the rate of dissolution is a function of the reactive surface area measured through the crack density in the volume. We use a gradient-enhanced damage framework in which damage depends on the deformation of the material as well as on the amount of solid mass dissolved over time. The gradient-enhanced formulation is characterized by a three-field variational formulation with the solid displacement, nonlocal equivalent strain, and nonlocal rate of solid dissolution as the basic state variables. Traditionally, time-independent damage theories have only allowed damage to increase with increasing external load. In the proposed framework, the degree of damage may increase due to solid dissolution even when the external load is held fixed. In this way, solid dissolution is viewed as a process that is responsible for bringing about rate-dependent effects such as creep and stress-relaxation, which are two common features of geomaterial behavior.
Optimum solutions for geologic modeling and reservoir simulation in industries such as oil and gas recovery and carbon capture and storage require accurate characterization of reservoir properties, which are often heterogeneous. In this study, high-quality micro-computed tomography (CT) images (1.475-μm/pixel resolution) of a sandstone core acquired from the Bell Creek oil field, USA, were used to provide nondestructive analysis of pore- and core-scale heterogeneity across measurement scales of 94–566 μm. In addition to characterizing the as-received sample, the core sample was flooded with brine to evaluate the capacity of the core sample to receive injected fluids. The micro-CT images were systematically segmented into pore spaces and grains via machine learning (ML) steps including image preprocessing, label creation using a traditional ML method based on limited manual image annotation, and finally U-Net segmentation. The segmented image stacks were reconstructed into digital cubes of various scales of voxel lengths. The 3D porosity values were calculated for all the digital cubes, and the fractal dimensions of the cubes were estimated using a box-counting method. The results showed that smaller cubes had greater heterogeneity and that the porosity values could be accurately estimated by fractal dimension and voxel lengths using ML models. For the core sample with brine flooding, the ratio of pores filled by brine to the total pore space was related to the porosity and could also be accurately estimated by porosity, fractal dimension, and voxel lengths using ML models. In conclusion, the results of this study demonstrate that the concept of fractal dimension can be a useful vector to perform upscaling analysis of sandstone rock heterogeneity from the pore to core scale and that fractal dimensions can be used to estimate porosity values and pore space-filling capacity across those scales.
Volume-increasing replacement reactions can lead either to fracturing of the parent phase or to the formation of cohesive layers that passivate further reaction, but the factors that drive one outcome or the other are not understood. In this experimental study, we investigated the volume-increasing replacement of carbonate rocks (consisting of CaCO 3 and CaMg(CO 3 ) 2 ) of different porosity by witherite (BaCO 3 ). Samples were characterized using scanning electron microscopy, Raman spectroscopy, small-angle neutron/X-ray scattering, X-ray tomography, and scanning transmission electron microscopy. We observed the formation of a witherite reaction layer and witherite formation within pores and along grain boundaries. Despite this being a volume-increasing replacement reaction, newly formed witherite was porous, potentially allowing further replacement. Filled fractures were observed in the low-porosity carbonates, whereas witherite formed within pores in high-porosity carbonates. We conclude that fracturing of the parent phase versus passivation is contingent on the initial microstructure of the parent with an optimal degree of porosity required for fracturing.