Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geologic Formations”

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.

At least 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Characterizing the uppermost 100 m structure of the San Jacinto fault zone southeast of Anza, California, through joint analysis of geological, topographic, seismic and resistivity data

We present results from complementary geological, topographic, seismic and electrical resistivity surveys at the Sagebrush Flat (SGB) site along the Clark fault (CF) strand of the San Jacinto fault zone trifurcation area southeast of Anza, California. Joint interpretation of these data sets, each with unique spatiotemporal sensitivities, allow us to better characterize the shallow (<100 m) fault zone at this structurally complex site. Geological mapping at the surface shows the CF has three main subparallel strands within a <100 m zone with varying degrees of rock damage. These strands intersect units of banded gneiss and tonalite, and various sedimentary units. Near the surface, the weathered but more intact tonalite and gneiss to the southwest have relatively high VP. The low-lying flat sedimentary basins around the two southwestern-most CF strands and elevated damaged gneiss to the northeast have lowest VP <500 m s–1. The high relief of the northeast gneiss unit may in part be explained by its extensive damage and inferred increased relative rock uplift. Resistivity imaging shows the unconsolidated dry basin sediments (maximum >1300 Ohm.m) contrasted against the compacted fine-grained (potentially wet) materials within the CF core and the Bautista Formation (minimum <40 Ohm.m), which is slightly elevated above the flat basins. The inverse relationship between VP (increases) and resistivity (decreases) in the uppermost ~15 m can be characterized as log–log linear with slopes of –2.6 to –4. At depths >30 m, the velocity heterogeneity near the surface merges into larger-scale structures that are generally slower on the northeast side of the CF core compared to the southwest side (as much as ~40 per cent reduction in average VP). A previous study revealed a 20–37 per cent variability in peak ground velocities across the SGB site from local earthquakes. The upper end of that range is associated with the near-surface unconsolidated sedimentary basins and northeast damaged gneiss unit. Preliminary analysis of time-dependent topography mostly shows effects of changing vegetation and anthropogenic activity.

Geochemistry & Geophysics↗
DOE OSTI2023

Effects of anisotropic shale creep on the stress and permeability evolution of a geological nuclear waste repository

To ensure long-term safety and performance, geological nuclear waste repositories require low-permeability barriers such as bentonite buffers and/or shale host rock. Shale is not only known for its low permeability but also for its trend to undergo time-dependent deformation (i.e., creep), which could heal damage, but the effects of shale creep on the long-term performance of nuclear waste repositories have not been clearly understood. In particular, the anisotropic nature of shale (i.e., bedding) could have a significant effect on its creep behavior, and consequently, on the long-term performance of nuclear waste repositories. In this research, numerical simulations were carried out with the objective of showing the effects of anisotropic shale creep on the stress and permeability evolution of a generic geological nuclear waste repository in shale. The TOUGH-FLAC simulator, a thermo-hydromechanically (THM) coupled numerical code, was used for the simulations. To achieve the objective, comparisons were performed between the results of anisotropic shale creep simulations and those of different simulation cases, namely, no creep (i.e., elastic), isotropic creep, and long-term creep shale cases. Results of the comparisons show that the elastic and isotropic creep shale cases respectively led to the overestimation and underestimation of stress and permeability in the repository, whereas the long-term creep shale case, which accumulated greater creep in later periods than in earlier periods, helped to keep large shear and tensile stresses from developing while maintaining compressive spherical stress, resulting in consistently low permeability levels. These results indicate that performance assessments with elastic and isotropic creep formation models will provide the upper and lower bound estimates of stress and permeability, while more reasonable estimates will be given by an anisotropic creep formation model, and that shale with long-term creep characteristics will be beneficial in many aspects of the safety and performance of nuclear waste repositories.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Influences of Mg Ions and Al Ions on Deposition Behavior of Supersaturated Silicic Acid - 20261

Silicic acid can change into various species such as deposit, colloidal, and soluble form under certain conditions: for example, gradual pH changes around geological disposal repositories. In particular, silicic acid deposits narrow groundwater flow paths, thereby restricting the flow of groundwater. Such effects may contribute to a delay in radionuclide migration. Previous studies have reported that silicic acid supersaturated in sodium and calcium ions can narrow the flow paths. This study examined the effect of magnesium (Mg) and aluminum (Al) ions on the deposition of supersaturated silicic acid on the solid phase. The deposition behavior of supersaturated silicic acid on the solid phase was investigated in a batch system while changing the initial supersaturation concentration, the coexisting ion concentration, and the amount of added solid phase. In the presence of Mg ions, silicic acid precipitation increased with increasing Mg ion concentration due to the influence of the electrolyte, and crosslinking reactions. The coexistence of Al ions significantly suppressed silicic acid deposition due to the formation of a covering on the solid surface during Al ion hydrolysis under conditions of high Al concentration (1.4 mM). On the basis of these experimental results, the apparent deposition rate constant of silicic acid was estimated to be in the range from 10{sup -11} to 10{sup -10} m/s. The Damkoehler number (the dimensionless number expressed by the ratio of the characteristic time of mass transfer and the chemical reaction constant) was sufficient to cause the clogging of the host rock by deposition. This suggests that supersaturated silicic acid deposition, regardless of the presence of Mg ions (and Al ions), contributes to improving barrier performance against radionuclide migration in the host rock. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Calculation of the Vertical Strata Load of Utility Tunnel Crossing Ground Fissure Zone

A ground fissure is a geological disaster in which the vertical dislocation of strata causes surface rupture. Ground fissures can cause extreme harm to the surface and underground buildings. Ground fissure activity can result in different settlement on the two sides of the strata, which will generate additional stress (pressure) that differs from the stress of the general stratum on underground structures across the ground fissure zone. It is essential to assess the effective stress of strata in the design of underground engineering structures across a ground fissure zone. The Xi’an ground fissure through a utility tunnel was focus of the research, and a physical model and data for oblique crossing of the 45° ground fissure were analyzed. A model of the utility tunnel structure was established, including the surrounding soil load as an active ground fissure environment. This model was used to calculate the vertical formation pressure of the overlying soil on the utility tunnel. A method to calculate the overlying load on the utility tunnel caused by ground fissure activity was proposed and compared with the calculation based on the A. Marston principle. The results showed that the ground fissure load calculation method based on the strata-holding effect can effectively calculate the earth pressure of the surrounding soil layer of the utility tunnel in the cross-ground fissure section. The results of this work provide guidance and reference value for the design of a utility tunnel in an area with the potential for a ground fissure.

Zhang, Dan↗
DOE OSTI2021

Early Holocene cold snaps and their expression in the moraine record of the eastern European Alps

Glaciers preserve climate variations in their geological and geomorphological records, which makes them prime candidates for climate reconstructions. Investigating the glacier–climate system over the past millennia is particularly relevant first because the amplitude and frequency of natural climate variability during the Holocene provides the climatic context against which modern, human-induced climate change must be assessed. Second, the transition from the last glacial to the current interglacial promises important insights into the climate system during warming, which is of particular interest with respect to ongoing climate change. Evidence of stable ice margin positions that record cooling during the past 12 kyr are preserved in two glaciated valleys of the Silvretta Massif in the eastern European Alps, the Jamtal (JAM) and the Laraintal (LAR). We mapped and dated moraines in these catchments including historical ridges using beryllium-10 surface exposure dating (10Be SED) techniques and correlate resulting moraine formation intervals with climate proxy records to evaluate the spatial and temporal scale of these cold phases. The new geochronologies indicate the formation of moraines during the early Holocene (EH), ca. 11.0 ± 0.7 ka (n = 19). Boulder ages along historical moraines (n = 6) suggest at least two glacier advances during the Little Ice Age (LIA; ca. 1250–1850 CE) around 1300 CE and in the second half of the 18th century. An earlier advance to the same position may have occurred around 500 CE. The Jamtal and Laraintal moraine chronologies provide evidence that millennial-scale EH warming was superimposed by centennial-scale cooling. The timing of EH moraine formation coincides with brief temperature drops identified in local and regional paleoproxy records, most prominently with the Preboreal Oscillation (PBO) and is consistent with moraine deposition in other catchments in the European Alps and in the Arctic region. This consistency points to cooling beyond the local scale and therefore a regional or even hemispheric climate driver. Freshwater input sourced from the Laurentide Ice Sheet (LIS), which changed circulation patterns in the North Atlantic, is a plausible explanation for EH cooling and moraine formation in the Nordic region and in Europe.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

In situ deformation of antigorite-olivine two-phase mixtures: Implications for dynamics and seismic anisotropy in the mantle wedge

Water released from hydrous minerals in subducting slabs reacts with the overlying plate, resulting in widespread serpentinization in the mantle wedge. Deformation of serpentinized peridotites has been invoked to explain forearc seismic anisotropy, yet studies of the mechanical properties and deformation behaviors of serpentine-bearing multiphase aggregates remain limited. Here we deformed olivine-antigorite mixtures containing 70, 50, and 20 vol.% of antigorite at 2.5 – 7.6 GPa, 673 K and strain rates of ∼10 –5 –10 –4 s –1 . Elasto-viscoplastic self-consistent simulations, constrained by synchrotron X-ray diffraction (XRD) data, were used to estimate lattice strain, stress–strain partitioning, crystallographic preferred orientations (CPO) development, and aggregate strength. Selected run products were also analyzed by electron backscatter diffraction for comparison with the CPO results obtained from XRD experiments. We found olivine transitions from A- or B-type to C-type when antigorite fraction drops to 20 vol.%, coinciding with a microstructural change from interconnected weak layers to a load-bearing framework (LBF). An additional run on a sample Atg50/Ol50 with preexisting microstructures suggested the formation of LBF was promoted by these microstructures, although the preexisting antigorite CPO has been overprinted at 20.8 % strain and could be erased completely by subsequent deformation in nature. Estimated viscosity of the two-phase mixtures suggests that low-degree serpentinization (≤20 %) in the mantle wedge may increase the strength of olivine-rich peridotite and hinder slab-mantle decoupling, whereas high-degree serpentinization (≥50–70 %) weakens the peridotite and favors decoupling if sufficient viscosity contrast (>10) develops. Seismic anisotropy shows a nonlinear dependence on antigorite fraction: antigorite CPO governs the anisotropy of the mixtures with ≥50 vol.% antigorite, whereas olivine CPO dominates at low fractions (∼20 vol.%). The presence of pre-existing microstructures reduces seismic anisotropy of the deformed mixtures, however the persistence of pre-existing CPO in actively subducting slabs remains uncertain, making their significance over geological timescales questionable.

Crystallographic preferred orientation↗
DOE OSTI2022

QA/QC of the East River, Colorado, discharge and geochemical time series datasets (Almont, BCC, and Pump House) to be used for modeling of hydrogeochemical balance

The following datasets were QA/QC-ed (Quality Assurance/Quality Control): 1. Brush Creek Confluence (BCC) discharge data (from Helen Malenda, USGS, Colorado School of Mines), which were calculated using the pressure transducer data and rating curves. The original 15 min time series data were presented as mean daily discharge. 2. Almont discharge data from United States Geological Survey (USGS). The original data were in 15 min time intervals, and were averaged to mean daily discharge time series. 3. Pump House discharge data as mean daily discharge (downloaded from the SFA portal). 4. BCC and Pump House chemistry data from SFA data portal and/or original spreadsheets provided by Roelof Versteeg. The following challenging QA/QC problems of the datasets were resolved: Missing data with the duration of gaps up to >1 month; Duplicated dates; Anomalies and outliers of discharge and concentrations; Time stamps of measurements of the discharge and concentrations are not aligned (hydrogeochemical balance calculations require the timestamps to be aligned). All QA/QC-ed datasets are given as csv files. The csv files were prepared using the xts files with multiple worksheets, which are also included in the data packages. Figures of the QA/QC-ed datasets are given in the jpeg and pdf formats. The QA/QC-ed datasets have been used to quantify discharge and chemical concentrations in river water in order to understand riverine exports of water and dissolved constituents in the East River watershed. These datasets served as a basis in the presentation given by P. Fox et al. at the 2021 Goldschmidt Conference.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Nano-Engineered Catalyst Supported on Ceramic Hollow Fibers for the Utilization of CO 2 in Dry Reforming to Produce Syngas

The objective of this project was to develop a novel catalytic reactor containing nano-engineered catalysts for the utilization of CO 2 (captured from coal-fired power plants and other CO 2 emitting sources) in dry methane reforming (DMR) (CO 2 + CH 4 → 2 H 2 + 2 CO) to produce synthesis gas (syngas). The technology aims to reduce CO2 emissions by developing beneficial uses for CO 2 from coal-fired power plants. It also offers an alternative to mitigate CO 2 emissions in areas where geologic storage may not be an optimal solution and/or utilization could significantly offset the costs of carbon capture and sequestration. The nano-engineered Ni-based catalyst was prepared by atomic layer deposition (ALD). The Ni particles were as small as ~2-4 nm. The nano-engineered catalyst showed CH 4 conversion >95%, H 2 /CO ratio in the range of 0.7-1.0, and CH 4 reforming rate as high as 2,500 L/h/gNi at 850 ºC and pressure of 15-25 psia. The Ni-based ALD catalyst also showed good stability in DMR reaction during a 200-h continuous operation at 850 °C. This is due to strong bonding between the nanoparticles and substrates since the Ni nanoparticles were chemically bonded to the substrate during the ALD process. The high thermal stability maintains the high dispersion of Ni nanoparticles, which can inhibit coke formation because their step edges are small enough to limit carbon nucleation and growth. Technoeconomic analysis (TEA) indicates the levelized cost of syngas (LCOS) is $172/ton with our technology, which is lower than the equivalent (molar) cost of hydrogen produced by steam methane reforming (SMR) or autothermal reforming (ATR). The major operating cost is natural gas feed and fuel, and the levelized cost is highly sensitive to the price of natural gas and relatively insensitive to the CAPEX. Revenues from syngas could have a significant impact on the net cost of electricity (COE), depending on the cost of natural gas and the selling price of syngas, estimated at $36 per MWh if the syngas were sold at $195 per ton. Following DOE NETL’s guidance, a lifecycle analysis (LCA) was conducted to compare with SMR. The functional unit for the basis of comparison was defined as 1kg carbon monoxide in the product stream. The global warming potential (GWP) of our process was found to be 40% lower than the state-of-the-art SMR process. The sensitivity analysis confirms the emissions are most sensitive to the natural gas fuel requirements to deliver heat to the process.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Effects of He Ion Irradiation in Simple Borosilicate Glasses for HLW Immobilisation - 20324

High-level radioactive waste (HLW) consists of fission products and minor actinides. Borosilicate glass is the host material of choice in most countries to safely immobilize HLW, followed by eventual geological disposal. HLW glasses will experience both self-heating and self-irradiation effects due to the decay of fission products and minor actinides present therein. During the first 500 years of geological disposal, beta-gamma decay will be the main source of radiation due to short half-lives of fission products, as compared to very long-lived actinides which will be responsible for irradiation over much longer timescales. High-energy alpha particles and alpha recoil nuclei can potentially have significant impact (stored energy, volume change, accumulated dose, chemical durability) on the performance of the wasteform over longer timescales. Thus, to ensure safe and economic disposal over geological periods (10{sup 4} - 10{sup 6} years), the assessment of performance and properties in terms of radiation effects becomes indispensable. This work aims to simulate the effects of alpha decay in radioactive waste glasses. Two international glasses, namely UK 'MW' lithium-sodium borosilicate and Indian sodium-barium borosilicate, which are being used as HLW immobilization matrices by the UK and India respectively, were prepared by a melt-quench method and annealed below their glass transition temperatures. Both glasses were irradiated with 650 KeV He{sup 2+} ions with doses ranging from 10{sup 15} to 10{sup 17} ions/cm{sup 2}. Structural analysis was carried out using electron paramagnetic resonance (EPR) and Raman spectroscopies. First derivative continuous wave EPR spectra were recorded for all pristine and irradiated glass monoliths at room temperature using an X-band (frequency 9.6 GHz) EPR Bruker spectrometer. Results from EPR analyses showed the formation of irradiation induced paramagnetic centers at g ∼ 2.00 in all irradiated glasses. Raman spectra of pristine and irradiated samples were acquired, and showed notable changes in glass structure post-irradiation, with the emergence of a broad Raman band at 1560 cm{sup -1} upon irradiation with high fluences. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

A simulation study of carbon storage with active reservoir management

As part of the Integrated Midcontinent Stacked Carbon Storage Hub (IMSCS-Hub) project led by Battelle Memorial Institute, a study was conducted to determine the feasibility of storing carbon dioxide (CO 2 ) in the stacked saline rock formations of the Sleepy Hollow Field (SHF), located in Red Willow County, southern Nebraska. A series of CO 2 injection simulation scenarios, with and without active reservoir management (ARM; brine extraction), were evaluated to investigate the feasibility of storing 50+ million tonnes (Mt) of CO 2 . The results indicated CO 2 injection combined with ARM may enable permanent storage of 50+ Mt of CO 2 . The area of review (AOR), the area in which underground sources of drinking water (USDWs) might be endangered during CO 2 injection, was assessed for the simulation scenarios. In comparison to a case without ARM, brine extraction resulted in a much smaller AOR, covering an area of 42 square miles (108.8 km 2 ), roughly one-fourth the size of an AOR resulting from CO 2 injection without ARM. Our findings, presented in this paper, indicate that CO 2 injection with ARM can improve the CO 2 storage capacity of a geologic storage complex up to 100% while also reduce the rate of pressure buildup in the subsurface, resulting in a 75% reduction in AOR. This may help in lowering carbon capture and storage project costs, risks, and effort needed to meet monitoring requirements for a storage project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Probabilistic risk-based Area of Review (AoR) determination for a deep-saline carbon storage site

Regulatory oversight of a geologic carbon sequestration (GCS) project relies on iterative estimations, throughout the project lifetime, of the area where increased risks to underground sources of drinking water (USDWs) may occur due to injection of CO2. This area, referred to as Area of Review (AoR), is typically delineated by predicting the migration of fluid between the reservoir and the lowermost USDW via an open wellbore using predictions from physics-based reservoir simulators. The inherent uncertainty in input parameters used in reservoir modeling therefore affects the accuracy of determining the AoR for a project. Furthermore, the standard analytical approaches for calculating a critical pressure to delineate the risk area yield an infinite AoR for cases where the injection reservoir is overpressured relative to the USDW. A methodology is presented in this paper to better characterize the risk to USDWs while accounting for the uncertainty in reservoir modeling, with an application to a permitted GCS project with an overpressured injection formation, FutureGen 2.0. The methodology is demonstrated using the National Risk Assessment Partnership’s open-source integrated assessment model (NRAP-Open-IAM) to develop a probabilistic estimate of impact risk to USDW quality. CO2 and pressure predictions from the reservoir modeling conducted using the STOMP-CO2 simulator for the FutureGen 2.0 site are used in a NRAP-Open-IAM model with reservoir, wellbore, and aquifer components to: (1) assess the extent of potential leakage into the USDW for the predicted reservoir pressure conditions; (2) evaluate the extent of potential impact using “no-net-degradation” thresholds; and (3) account for uncertainty in reservoir permeabilities.

03 NATURAL GAS↗
DOE OSTI2020

Viscosity of the Mercurian Magma Ocean: Implications of Sulfur-Free and Sulfur-Bearing Magma Oceans for Differentiation and Crustal Petrogenesis

Mercury is believed to have had a magma ocean in its early evolution whose solidification dynamics contributed to the planet’s present day surface composition. The surface is characterized by two main geological provinces; the Northern Smooth Plains (NSP),and the Heavily Cratered Terrain (HCT) and Intercratered Terrian (IcT). These provinces are defined by the differences in morphology, composition, and age. The NSP are located in the northern hemisphere and have low Mg/Si, Ca/Si, and S/Si ratios and high Al/Si ratios, and are compositionally analogous to a magnesian alkali basalt. The HCT/IcT cover the majority of the planet and are older than the NSP, and are characterized by high Mg/Si, Ca/Si, and S/Si ratios and low Al/Si ratios. The dichotomy in surface composition motivates the study into the structure of the mantle as a product of magma ocean solidification. Two important and interrelated parameters that determine the efficiency of crystal fractionation during magma ocean solidification are convection vigor and viscosity of the magma ocean. Crystal suspension in a cooling magma ocean could be efficient or inefficient which suggests two end member possibilities for magma ocean solidification. The first is the formation of a compositionally stratified mantle reflecting a cumulate sequence produced along a fractional liquid line of descent. The second is the formation of a well-mixed, batch crystallized mantle by high viscosity and/or vigorously convecting magma ocean. Determining which of these end members, or where on the spectrum between them is most relevant to Mercury’s magma ocean can be experimentally constrained through high P-T falling-sphere viscometry experiments.

Mouser, M. D.↗
DOE OSTI2020

A spatial-statistical investigation of surface expressions associated with cyclic steaming in the Midway-Sunset Oil Field, California

In the Midway Sunset Oil Field in Central California, operators inject steam into the shallow diatomite formation to enhance heavy oil recovery through imbibition, wettability alteration, and viscosity reduction, among other mechanisms. The injected steam, however, does not always remain in the reservoir or return through the wells. In two zones in the study area, the steam comes out at the surface, creating sinkholes, seeps, and steam outlets. These phenomena, called “surface expressions,” pose safety and environmental hazards. Even though these surface expressions are a widespread problem in Central California, they are not well documented and understood. Possible causes of the surface expressions include: high injection pressure, structurally controlled flow patterns, leakage of steam through old improperly abandoned wells, high injection volumes, or flow along naturally occurring faults, among other possible factors. This work examines attributes of the zones with surface expressions in order to determine factors that may contribute to their occurrence. Spatial statistical analysis using logistic regression, random forests, and classification trees is used to explore the relationship between the surface expressions and geological and production-related attributes. The results point to a significant spatial correlation between the surface expressions and two predictors: concentration of plugged wells and geologic seal thickness. The results guide follow-up studies to further investigate the role of well abandonment and seal thickness in the occurrence of surface expressions.

02 PETROLEUM↗
DOE OSTI2021

Density Functional Tight-Binding Simulations Reveal the Presence of Surface Defects on the Quartz (101)–Water Interface

Understanding the structure and reactivity of quartz–water interfaces is critical for numerous applications in the geological, environmental, and biological sciences. However, disagreements on the atomic-level structure of the interfaces between experiments and simulations are hampering our ability to predict the surface reactivity. Here, we used density functional tight-binding (DFTB)-based molecular dynamics to simulate a series of quartz (101) surfaces having different types and densities of surface defects in water and compared them with the structures determined by X-ray reflectivity measurements. The DFTB simulations are able to reproduce previous classical and quantum mechanical predictions of the pristine quartz (101)–water interface that disagree with experimental observations. To remedy this situation, a set of defective quartz surfaces having various surface silicon (Si) vacancies were built as indicated by recent experimental studies. We found that the rotation of surface [SiO 4 ] tetrahedra near Si vacancies can lead to outward displacements of Si atoms similar to those observed in the experiments. The presence of additional surface Si vacancies caused inward relaxations of terminal oxygens through the formation of hydrogen bonds. Finally, the overall results indicate that the quartz (101)–water interface may include a mixture of geminal (≡Si–(OH) 2 )- and vicinal (≡Si–OH)-type silanol groups together with the presence of surface Si vacancies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Thermal Characterization and Exhumation of Northwest San Juan Basin Area, NM

The San Juan Basin (SJB), located in southwestern Colorado northwestern New Mexico, containing Paleozoic through middle Cenozoic strata, formed as a partitioned basin approximately 80 Ma in response to the Laramide orogeny. The SJB is a commercially mature, petroleum-producing basin, and currently being explored for CCUS and geothermal resources. Modern heat flow within the SJB is spatially variable with higher geothermal gradients in its northern and eastern portions. The temporal history of the basin thermal history is essential for understanding capacity for carbon storage and geothermal exploration. We present thermochronometric analyses of two surface and two subsurface samples to constrain spatial-temporal thermal evolution of the northwestern San Juan Basin. Apatite (U-Th)/He thermochronometric analyses were conducted on Cliff House and Kirtland Formation outcrops and subsurface samples of the Ojo Alamo and Pictured Cliffs Formations from a DOE funded, CCUS project pilot-well. In addition, vitrinite reflectance data and 1D basin modelling constrain possible time-temperature pathways of the samples in an area where no previous studies have constrained the uplift/exhumation of the Hogback Monocline using thermochronometry. Modeling suggests that elevated subsurface temperatures developed simultaneously with regional late Oligiocene volcanism in the San Juan Volcanic field and timing of maximum burial in the SJB. Late Miocene to Pliocene cooling/exhumation through the apatite (U-Th)/He partial retention zone support stratigraphic evidence of the Colorado Plateau uplift, possibly regional epeirogenic uplift driven by mantle processes.

02 PETROLEUM↗
DOE OSTI2020

SECARB Semi-Annual Field Review Reports - Early Test

The Southern States Energy Board (SSEB) managed the overall Southeast Regional Carbon Sequestration Partnership (SECARB) Phase III project and The Bureau of Economic Geology (BEG), at the University of Texas at Austin, managed all activities associated with the Early Test field site near Cranfield, Mississippi (the Early Test or Cranfield Project). BEG was assisted in the management of site activities by Sandia Technologies (Sandia). BEG and Sandia implemented reporting procedures that included monthly and quarterly technical progress reports, as well as periodic reporting on key field activities. Subsequent sections of this report catalog the field activities and are organized in a semiannual chronology. The Cranfield Project was located approximately 12 miles east of Natchez, Mississippi. Denbury Onshore, LLC, began operating a commercial CO 2 flood of the field (using the subsurface injection of CO 2 for EOR) in June 2008. Carbon dioxide from the Jackson Dome, a natural source near Jackson, Mississippi, was delivered to the Cranfield oilfield via pipeline. Beginning in October 2008, the SECARB Early Test team characterized the surface and subsurface of the Cranfield site. During the field test CO 2 was injected into the lower Tuscaloosa Formation, a regionally extensive saline formation with the potential to hold millions of tons of CO 2 emissions. By January 2015, the team injected more than 11 million metric tons of CO 2 and monitored a total stored CO 2 mass of 5,326,975 metric tons. SECARB’s study operations occurred in four integrated research program areas within Cranfield field: (1) the High Volume Injection Test area (HiVIT); (2) the Detailed Area of Study (DAS); (3) the Geomechanical area; and (4) the near-surface observatory, also called the “P-site.” Carbon dioxide injection activities were conducted at the HiVIT and the DAS. Figure 1.1-1 provides a depiction of the DAS.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2024

DECOVALEX-2023: Task F2 Salt Final Report

The subject of Task F of DECOVALEX-2023 concerns performance assessment modelling of radioactive waste disposal in deep mined repositories. The primary objectives of Task F are to build confidence in the models, methods, and software used for performance assessment (PA) of deep geologic nuclear waste repositories, and/or to bring to the fore additional research and development needed to improve PA methodologies. In Task F2- (salt), these objectives have been accomplished through staged development and comparison of the models and methods used by participating teams in their PA frameworks. Coupled-process submodels and deterministic simulations of the entire PA model for a reference scenario for waste disposal in domal salt have been conducted. The task specification has been updated continuously since the initiation of the project to reflect the staged development of the conceptual repository model and performance metrics. Thermal, hydrological, mechanical, and chemical properties of individual components of the engineered and natural system were chosen for relevance by participating teams. The salt reference case system was characterized using data and measurements collected at relevant underground research laboratories (URLs), field sites, and simulation results from teams with specialized modelling capability. Participating teams made a wide range of model assumptions from compartmentalized networks to full 3D models of the salt formation. No single contributed model includes full-fidelity representation of all the features, events, and processes (FEPs) detailed in the task specification, but almost all features and processes are represented in at least one model. Despite differences in the modelling strategies developed by participating teams, all models indicate that salt compaction and radionuclide diffusion are key processes in the repository, and for the FEPs and model scenario considered, little of the disposed radionuclides will migrate beyond the repository seal over the 100,000 year simulations. In general, the model output quantities have the largest differences over the short term and near the waste. The models tend to be more similar further from waste and at later time. Disparities between the models are believed to be due to differing simplifications from the task specification, some of which are chosen simplifications to reduce complexity, and some are restrictions imposed by the modelling tools. A second round of this task has been accepted for DECOVALEX-2027 in conjunction with Task F1 on crystalline PA modelling. The future round includes waste package heating, improved modelling of salt creep closure, additional comparisons of coupled-process sub-models, and the impact of repository engineering design on radionuclide migration in the repository. Participants will also propose and finalize a set of uncertain inputs for the reference case simulations, propagate these uncertainties in a set of realizations, and conduct sensitivity analyses on the simulation results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Data from: 'Abiotic influences on continuous conifer forest structure across a subalpine watershed'

This package archives the core data used for analysis and inference in 'Abiotic influences on continuous conifer forest structure across a subalpine watershed' (Worsham et al., 2025). All data were collected in the East River, Washington Gulch, Slate River, and Coal Creek watersheds of Colorado. In the paper, we quantified the relative influence of climate, topographic, edaphic, and geologic factors on conifer stand structure and composition, and their functional relationships, at the watershed scale. We used waveform LiDAR data to derive spatially continuous stand structure metrics. We fused these with a species-level classification map to estimate tree species abundance. We applied generalized additive and generalized boosted models to evaluate the covariability of structural and compositional metrics with abiotic variables. The package contains the essential products required for reproducing our analysis and the tables and figures reported in the publication. The products comprise four classes: (1) geospatial data, (2) tabular data used for inferential analysis, (3) tabular data describing analytical results and performance statistics, and (4) a data user guide. (1) includes discretized waveform LiDAR data, locations and attributes of individual tree crowns, sampling locations and domain boundaries, a canopy height model, and raster files of estimated forest structural and compositional metrics at 100 m grid scale. (2) includes all response and explanatory variable values applied in inferential models. Response variables include conifer forest stand density, basal area, 95th percentile height, quadratic mean diameter, and others. Explanatory variables include climatic water deficit, actual evapotranspiration, elevation, heat load, soil available water content, and others. (3) includes results of training and testing several individual tree detection (ITD) algorithms, as well as inferential modeling results. (4) is a PDF user guide for this data package, including detailed descriptions and data dictionaries for all files. The data package root contains 17 assets: 8 compressed tape archive (.tar.gz) files, 5 comma-separated values (.csv) files, 3 Geographic Tagged Image File Format (GeoTIFF) (.tif) files, and 1 Portable Document Format (.pdf) file. The compressed .tar.gz archives contain ESRI shapefiles (.shp) .tif, compressed LASer (.laz), and .csv files. The archives must first be decompressed using the widely distributed command-line software utility TAR. All other files, including constituent files within the .tar.gz archives, can be opened in the open-source R statistical computing environment. Alternatively, .csv files may also be read in any simple text editor software or Microsoft Excel. Geospatial files including .shp and .tif files can also be opened in GIS software, such as QGIS (open-source) or ESRI ArcGIS (proprietary). The .pdf Data User Guide can be read with Adobe Acrobat Reader or other compatible readers.

2018 NEON and 2025 CHESS Campaigns↗