SEARCH · Engineering Papers

Results for “rock”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Geomechanical characterization of Rock Valley carbonates

In Rock Valley, Nevada, the USA, past earthquake sequences and shallow faulting in Paleozoic carbonates remain poorly understood. Carbonates are typically more ductile rocks that do not experience large stress drops or fracture coalescence, which contradicts previous observations in the region. As part of the Source Physics Experiment, an effort has been made to experimentally characterize the petrophysical and geomechanical properties of carbonates from Rock Valley. Well core and outcrop samples were used to determine the difference between shallow-buried Tertiary limestones and deeply buried Paleozoic limestones and dolostones. The carbonates possessed porosities between 3 and 9%, with V P and V S in the Tertiary carbonates around 2700 and 1800 m/s, respectively, and 6000 and 3100 m/s, respectively, in the Paleozoic carbonates. Microstructural characterization revealed that the Paleozoic carbonates contained significantly more pre-existing damage and alteration than the Tertiary carbonates, though the deformation varies from localized to diffuse. Unconfined Brazilian tests and triaxial tests with confining pressures between 0 and 50 MPa showed that, although samples all experienced failure, the dolostones typically failed at greater stresses and possessed greater E/ν ratios than the limestones. Furthermore, the Hoek–Brown failure criterion was used to construct failure envelopes with the compressive and tensile failure tests. Velocity, density, and porosity measurements were used to construct a hypothetical velocity-depth profile and compare the results with the theoretical petrophysical measurements. Using brittleness analysis, the likely failure conditions were determined to be low-porosity dolomitic rock for fault nucleation in the Paleozoic basement.

Brittleness↗

DOE OSTI2024

Validating the rhenium proxy for rock organic carbon oxidation using weathering profiles

Chemical weathering over geological timescales acts as a source or sink of atmospheric carbon dioxide (CO 2 ), while influencing long-term redox cycling and atmospheric oxygen (O 2 ) at Earth's surface. There is a growing recognition that the oxidative weathering of rock organic carbon (OC petro ) can release more CO 2 than is locally drawn down by silicate weathering, and may vary due to changes in erosion and climate. The element rhenium (Re) has emerged as a proxy to track the oxidative weathering of OC petro , yet uncertainties in its application remain namely that we lack a systematic assessment of the comparative mobility of Re and OC petro during sedimentary rock weathering. Here we measure Re and OC petro loss across gradients in rock weathering at 9 global sites, spanning a range of initial OC petro values from ~0.2 % to >10 %. We use titanium to account for volume changes during weathering and assess Re and OC petro loss alongside major elements that reflect silicate (Na, Mg), carbonate (Ca, Mg) and sulfide (S) weathering. Across the dataset, Re loss is correlated with OC petro loss but not with loss of any other major element. Further, across the weathering profiles, the average molar ratio of OC petro to Re loss was 0.84 ± 0.15, with 8 out of 9 sites having a ratio >0.74. At one site (Marcellus Shale), the average ratio was lower at 0.58 ± 0.11. The excess loss of Re matches expectations that, typically, between ~0 and 20 % of the Re liberated by sedimentary rock weathering derives from silicate or sulfide phases, while some OC petro may be physically or chemically protected from weathering. Overall, our measurements provide validation for the Re proxy of OC petro oxidation and allow future work to further improve our knowledge of regional and global-scale rates of this important source of CO 2 in the geochemical carbon cycle.

58 GEOSCIENCES↗

DOE OSTI2021

Mineralogy-dependent sulfide oxidation via polysulfide and thiosulfate pathways during weathering of mixed-sulfide bearing mine waste rock

Tracking the S oxidation pathway in sulfide-bearing mine waste-rock piles is complicated by variations in water content, O 2 and Fe 3+ concentrations, microbial diversity, mineralogy, the occurrence of a variety of S species associated with incomplete oxidation, and nonlinear coupling between physicochemical processes. Synchrotron-based S K-edge X-ray absorption near edge structure (XANES) spectroscopy facilitates the identification and quantification of S species (i.e., S 2– , $S_{2}O_{2}^{2–}$, S 8 , $S_{2}O_{3}^{2–}$, and $S_{2}O_{4}^{2–}$) in a weathered mine waste-rock pile containing pyrite, sphalerite, galena, pyrrhotite, and chalcopyrite. Mineralogy-dependent polysulfide and thiosulfate pathways both affect the S oxidation within the waste-rock pile. Currently, the polysulfide pathway, with S 8 as the intermediate S species, generates high concentrations of dissolved metals (from the rapid oxidation of monosulfides including sphalerite, galena, pyrrhotite, and chalcopyrite). As these monosulfides are depleted, the thiosulfate pathway of pyrite oxidation, including intermediate oxidation products (i.e., S 8 and $S_{2}O_{3}^{2–}$), will become the dominant oxidation pathway, resulting in the potential for generation of additional acidic drainage. This information highlights the importance of identifying the sulfide-mineral oxidation reaction pathways when attempting to estimate acid generation potential, and when developing remediation strategies for the storage and management of sulfide-bearing waste rock.

58 GEOSCIENCES↗

DOE OSTI2020

Determination of permeability for hydrocarbon release due to excavation-induced stress redistribution in rock salt

Due to a stress redistribution after the excavation of an underground tunnel for radioactive waste disposal, an Ed/DZ (excavation disturbed/damaged zone) will be generated in the near field of the opening, resulting in significant changes in the hydraulic and mechanical properties of the rock mass in the zone. Initially more or less randomly distributed hydrocarbons at grain boundaries in rock salt, which sometimes can only be observed with ultraviolet light, can then be mobilised and migrate at a potentially significant rate towards the excavation. Within the international cooperative project DECOVALEX 2019, the migration mechanism of such fluid inclusions in rock salt is being studied intensively. Specifically, a multi-scale modelling strategy has been developed. A macroscale coupled hydro-mechanical modelling of an underground excavation was performed to determine hydraulic and time-dependent deviatoric stress conditions, by taking into account the rock salt creep behaviour. Under the obtained macro-scale constraints, micro-scale modelling of a pathway dilation along halite grain boundaries was performed using different model strategies: a) coupled hydromechanical modelling with a consideration of hydraulic pressure-induced dilatant deformation, b) nonlinear dynamic model taking account of fluid migration, stress-dependent grain boundary wetting and shear-induced dilatancy of salt, and c) phase-field modelling of flow pathway propagation. The permeability increase resulting from the pathway dilation is estimated to be as high as two orders of magnitude. Based on the permeability determined, a series of pressure buildups measured from a borehole with a high hydrocarbon release rate, a total of 430 build-ups within a monitoring time of 938 days, can be simulated with a macro-scale compressible flow model accounting for different zones around the opening.

58 GEOSCIENCES↗

DOE OSTI2022

Characterization of Acoustic Emissions From Analogue Rocks Using Sparse Regression‐DMDc

Abstract Moisture loss in rock is known to generate acoustic emissions (AE). Phenomena that result in AE during drying are related to the movement of fluids through the pores and induced‐cracks that arise from differential mineral shrinkage, especially in clay‐bearing rock. AE from the movement of fluids occurs from the reconfiguration of fluid interfaces during drying, while AE from mineral shrinkage involves the debonding within or between minerals. Here, analogue rock samples were used to examine the differences in the AE signatures when one or both AE source‐types are present. An unsupervised sparse regression model, Dynamic Mode Decomposition with control, that extends Dynamic Mode Decomposition is used to characterize the AE signals recorded during the drying of porous analogue rock samples fabricated with ordinary Portland cement, with and without clay. This method can effectively and accurately reconstruct acoustic signals emitted from samples that only experience moisture loss without cracking. However, the method struggles to reconstruct signals from samples with intricate crack networks that formed during drying because AE generating mechanisms can emit contemporaneously, and the resulting waves propagate through drying‐induced cracks that can lead to multiple internal reflections. Thus, the differential reconstruction accuracy of time series generated by different underlying physical processes provides a robust filter for reducing large data catalogs. In general, both dynamics and sparse initiating events are learned directly from data and this method exposes a data hierarchy based on the complexity of the intrinsic dynamics.

58 GEOSCIENCES↗

DOE OSTI2022

Structural Characterization of a Toppling Rock Slab From Array-Based Ambient Vibration Measurements and Numerical Modal Analysis

Accurate assessments of the internal structure and boundary conditions of unstable rock slopes are imperative for evaluating landslide hazard scenarios. However, instability characterization at depth remains challenging and is often limited by costly or invasive subsurface investigations. Here, we develop a new approach coupling array-based ambient vibration modal analysis and numerical modeling to improve structural characterization of rock slope instabilities at depth. We used ambient noise cross-correlation on 4 hr of seismic data recorded by an array of 30 nodal geophones at a 500-m-long toppling rock slab in Utah, USA to identify modal frequencies between 0.8 and 3.5 Hz and derive modal displacements. We show that transverse and longitudinal bending modes span the length of the instability, indicating an interconnected slab. Further, statistical comparison of field results with outputs from >1,000 finite element models with varying boundary conditions showed that the instability depth varies between 40–70 and 10–20 m in the middle and lateral regions, respectively. Our approach yields new information on the structural conditions of rock cliff and column instabilities at depth, which is not easily obtained by other means but is imperative for change detection monitoring and improved hazard assessments.

58 GEOSCIENCES↗

DOE OSTI2026

Stimulated geologic hydrogen: from mechanistic control to engineered rock transformation

Geologic hydrogen (GeoH 2 ) generated from subsurface iron-rich rock–water reaction (i.e., serpentinization) is emerging as a promising candidate for the next primary energy source. Yet, accelerating GeoH 2 production from geological to human timescales via enhanced serpentinization remains a formidable scientific and technical challenge. Here, in this Review, we decipher the mechanistic control and explore strategies for accelerating in situ, engineered iron-rich rock transformation into carbon-free GeoH 2 by orders of magnitude. Serpentinization rate is hindered by low porosity and permeability of source rocks, suboptimal temperatures, unfavorable water chemistry, inefficient Fe 2+ -to-Fe 3+ conversion, thermodynamic constraints, and low reactive surface area. While closed-system experiments provide valuable mechanistic insights, open-system conditions with fluid circulation are more crucial for economically viable GeoH 2 production. We assess stimulation techniques from enhanced hydrocarbon and geothermal recovery as tools to be adopted or adapted for increasing reactive surface areas for stimulated GeoH2 production. We estimate that 7.40 × 10 5 to 1.73 × 10 6 million metric tons (Mt) of hydrogen could be engineered over 20 to 50 years from about 10% iron-rich rocks within 10 km depth of continental crust. Enabling GeoH 2 as a viable energy source requires not only advancing scientific frontiers but also forming a global GeoH 2 research network and innovation ecosystem to address the critical scientific, technical, societal, economic, and policy challenges.

08 HYDROGEN↗

DOE OSTI2020

Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies

The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress in modeling and experimental approaches towards the characterization of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in deep clay/shale/argillaceous rock. International collaboration activities such as heater tests and postmortem analysis of samples recovered from these have elucidated key information regarding changes in the engineered barrier system (EBS) material exposed to years of thermal loads. Chemical and structural analyses of sampled bentonite material from such tests has as well as experiments conducted on these are key to the characterization of thermal effects affecting bentonite clay barrier performance and the extent of sacrificial zones in the EBS during the thermal period. Thermal, hydrologic, and chemical data collected from heater tests and laboratory experiments has been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches (e.g., continuum vs. discrete) to tackle issues related to flow and transport at various scales of the host-rock and EBS design concept. Consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the EBS design. This is particularly important for backfilled repository concepts where temperature plays a key role in the EBS behavior and long-term performance. This report describes multiple R&D efforts on disposal in argillaceous geologic media through development and application of coupled THMC process models, experimental studies on clay/metal/cement barrier and host-rock (argillite) material interactions, molecular dynamic (MD) simulations of water transport during (swelling) clay dehydration, first-principles studies of metaschoepite (UO 2 corrosion product) stability, and advances in thermodynamic plus surface complexation database development. Drift-scale URL experiments provides key data for testing hydrological-chemical (HC) model involving strong couplings of fluid mixing and barrier material chemical interactions. The THM modeling focuses on heater test experiments in argillite rock and gas migration in bentonite as part of international collaboration activities at underground research laboratories (URLs). In addition, field testing at an URL involves in situ analysis of fault slip behavior and fault permeability. Pore-scale modeling of gas bubble migration is also being investigated within the gas migration modeling effort. Interaction experiments on bentonite samples from heater test under ambient and elevated temperatures permit the evaluation of ion exchange, phase stability, and mineral transformation changes that could impact clay swelling. Advances in the development, testing, and implementation of a spent nuclear fuel (SNF) degradation model coupled with canister corrosion focus on the effects of hydrogen gas generation and its integration with Geologic Disposal Safety Assessment (GDSA). GDSA integration activities includes evaluation of groundwater chemistries in shale formations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Modular, Crushed-Rock Thermal Energy Storage Pilot Design (Final Report)

The goal of this project was to design a next-step pilot to advance near-term energy storage integrated with a fossil plant to provide a facility capable of being viable and effective in a market with growing penetration of variable renewable energy (VRE). Thermal energy storage (TES) represents an ideal technology for this purpose. The completed effort included a feasibility study to prepare for the Phase II pre-front end engineering design (pre-FEED) for implementing a crushed-rock TES system integrated with a natural gas combined cycle (NGCC) plant. The crushed rock storage technology, which is being developed by Brenmiller Energy, is a modular TES system termed bGen™, which can accommodate both thermal and electrical inputs and output steam, hot water, or hot air. For this application, the estimated efficiency is 80% thermal to thermal. For the feasibility study, the Brenmiller technology was designed to operate on a slipstream from NYPA’s Eugene W. Zeltmann Power Project (Zeltmann) NGCC plant in Astoria, New York. The projected size of the system will be up to 4 MWe with at least 4 hours of storage duration, or 16 MWh-e total. The study also included a techno-economic evaluation of a 200 MWh commercial-scale demonstration. Prior to this project, EPRI had reviewed Brenmiller’s technology, which is being built to demonstrate bGen™ at 1.7 MWe on a solar plant (Rotem) and has been designed for an NGCC facility in Italy, assessing it at technology readiness level (TRL) 5. Brenmiller is also conducting a separate 1-MWth pilot with NYPA that pairs a bGen™ module with a microturbine for a combined-heat-and-power (CHP) application to improve efficiency and provide flexibility. The next-step pilot being designed as part of this project would represent a 5-fold increase in scale, versus Rotem, and would show the technology’s ability to provide effective and economical energy storage, bringing the technology to TRL 6. This pilot would represent the next-to-last demonstration scale before the technology could be commercial ready at GWh-e scales in the 2030 timeframe. The main objective of the work completed by the Electric Power Research Institute, Inc. (EPRI), Brenmiller Energy (Brenmiller), New York Power Authority (NYPA), and United E&C (formerly AECOM) was to perform a Phase I feasibility study on the integration of a crushed-rock thermal energy storage (TES) with a fossil plant. Under this project, the EPRI-led team successfully completed a feasibility study to prepare for the potential future Phase II pre-front end engineering design (pre-FEED) to implement a crushed-rock TES system integrated with a natural gas combined cycle (NGCC) plant. Specific deliverables under this project included the Technology Maturation Plan, Conceptual Study, Techno-Economic Assessment, Technology Gap Assessment, Project Plan for Phase II (submitted as Phase II Renewal Application), Commercialization Plan, and the Final Report (this document). This Final Report includes a compilation of the various summary reports that were prepared during the 12-month schedule of the Phase I project execution under award DE-FE0032017.

03 NATURAL GAS↗

DOE OSTI2023

The role of stress and fluid saturation on the acoustic response of fractured rock

Standard rock physics models are formulated to describe the behavior of porous sedimentary reservoirs, with clean sandstones being the archetypal system; however, many situations demand geophysical monitoring of rocks with significantly different structures, such as low porosity, fractured reservoirs. Conventional models also suggest that these “stiff” reservoirs can be challenging to monitor seismically due to small fluid substitution effects, but the presence of fractures leads to stress dependence which may be leveraged for remote monitoring purposes. Using samples from the Duperow Formation (dolostone) obtained from the Danielson test well in Kevin Dome, MT, we conducted ultrasonic and multi-scale structural (profilometry, synchrotron micro-tomography, pressure sensitive film) measurements on naturally fractured core in order to characterize the effects of fluid substitution and effective stress on the acoustic response of fractured reservoir rock with a focus in particular on the textural and seismic characteristics of natural fractures. We find that changes in effective stress can yield changes in velocity of up to 20% and changes in attenuation up to 200%. Measured fluid substitution effects are resolvable, but stress effects dominate. These measurements provide insight into the physical processes controlling acoustic response of fractured rocks in general and can also be used to inform monitoring efforts in fractured reservoirs.

58 GEOSCIENCES↗

DOE OSTI2025

Helium Release During Fracture and Granular Fragmentation of Rocks

Geogenic Helium-4 ( 4 He) in-situ increases locally in regions of large deformation generated naturally or anthropogenically. This gas release by deformation is a potential geochemical precursor signal for subsurface deformation. To evaluate the applicability of 4 He degassing for correlating deformation in different lithologies, we conducted high force crush tests, up to 97,800 N axial load, to assess the total 4 He released during fragmentation of the rocks. We observed that the highest 4 He released occurred in the sedimentary rocks and that release correlated strongly with lithologic age and U/Th content. Microstructural changes of the pre- and post-test rocks indicate that the degree of grain size reduction relates directly to the total 4 He released during crushing. The range of in-place 4 He was calculated based XRF measurements of uranium and thorium in each lithology, with the results indicating that the majority of the trapped 4 He was not released. However, the 4 He released by deformation depended upon how the each rock deformed during deformation and the degree of grain size reduction. We postulate that 4 He precursor signals can be used to understand subsurface deformation only if geomechanical and geochemical conditions for 4 He enrichment in a lithology are met.

Deformation signals↗

DOE OSTI2021

Borehole breakout modeling in arkose and granite rocks

The existence of a deep borehole in the Earth’s crust disturbs the local stresses and creates a stress concentration that may result in breakout and damage to the borehole. Maintaining wellbore integrity mitigates environmental impacts such as groundwater contamination, gas leakage to the atmosphere, and fluid spills and seepage at the surface. In this paper, the stability of deep boreholes (5 km) is examined by laboratory experiments and numerical models in the context of nuclear waste disposal in Israel. Two rock types in southern Israel are considered: the crystalline basement (granite) and the Zenifim Formation (arkose). A series of room-temperature triaxial rock deformation experiments were conducted at different confining pressures. This mechanical characterization was then used to parameterize the elastic properties and damage behavior of the rocks. This facilitated modeling the stability of the deep boreholes by two different formulations of damage rheology: a dynamic-oriented formulation used to model deformation immediately after the creation of the open hole and a quasi-static formulation used to model longer stress corrosion regime. The calibrated modeling results indicate greater stability with Zenifim arkose than the crystalline granite for deep borehole conditions despite the granite having a greater triaxial compressive strength. Dissipation associated with dilation and porous compaction in the arkose during deformation plays a significant stabilizing role in the borehole compared to crystalline rocks. These results suggest that common strength-based borehole stability assessment may lead to inaccurate predictions. Three-dimensional modeling of bottom-hole stress conditions and the effects of transient borehole geometry show conventional two-dimensional analysis may not be conservative when predicting borehole damage.

58 GEOSCIENCES↗

DOE OSTI2020

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗

DOE OSTI2020

The formation and evolution of the Moon’s crust inferred from the Sm-Nd isotopic systematics of highlands rocks

Ages determined for magnesian and ferroan anorthosite crustal rock suites overlap, suggesting they formed contemporaneously about 4.3–4.5 Ga. A notable exception is the Sm-Nd age previously determined on Mg-suite gabbronorite 67667 which is at least 100 Ma younger than the youngest ferroan anorthosite. New chronologic measurements of 67667 presented here yield concordant Sm-Nd and Rb-Sr mineral isochron ages of 4349 ± 31 Ma and 4368 ± 67 Ma, suggesting the sample is older than previous estimates. Furthermore, a whole rock Sm-Nd isochron of Mg-suite rocks from the Apollo 14, 15, 16, and 17 landing sites yields an age of 4348 ± 25 Ma, indicating that Mg-suite magmatism was widespread and roughly contemporaneous on the lunar nearside. Here, analysis of Sm-Nd internal isochron ages confirms that Mg-suite magmatism was restricted to a period between about 4.33 and 4.35 Ga at the Apollo 14, 15, 16, and 17 landing sites and was synchronous with magmatism at the Apollo 16 site associated with the ferroan anorthosite suite between 4.35 and 4.37 Ga. Magnesian- and ferroan anorthosite suite rocks with ages younger than ~4.33 Ga appear to have experienced slow cooling in the deep lunar interior, so that the ages record when the samples cooled below the closure temperature of the Sm-Nd isotopic system and not the time they crystallized.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Temporal relationships among lunar crustal rocks

Temporal relationships among the three most common suites of lunar crustal rocks have been investigated by obtaining new high precision ages on Felsic/Alkali-suite Quartz monzodiorite Clast B from breccia 15405 and Magnesian-suite norite 78235/6/8/55/56 and comparing them to previously dated ferroan anorthosite sample 60025. The weighted average age of 4337.19 ± 0.49 Ma of 15405 Clast B is defined by zircon U-Pb and Pb-Pb ages as well as mineral isochron Sm-Nd and Nd-Nd ages. It is identical to the weighted average age for Apollo 17 norite 78235/6/8/55/56 of 4334.1 ± 3.5 Ma which is defined by Pb-Pb ages measured on baddeleyites in this investigation and less precise Pb-Pb and Sm-Nd ages reported in the literature. Both ages are ∼ 25 Ma younger than the weighted average of Sm-Nd and Pb-Pb ages reported in the literature on ferroan anorthosite 60025 of 4359.3 ± 2.3 Ma. The fact that ages of all three samples are defined by multiple U-Pb, Pb-Pb, Sm-Nd, and 142 Nd- 143 Nd chronometers provide confidence that they record the igneous crystallization history of the samples and do not represent disturbances or mixing lines with no temporal significance. Here, the extent to which these three ages represent broader scale magmatism is difficult to evaluate. Nevertheless, the age defined for 15405 Clast B, 78235/6/8/55/56, and 60025 are contemporaneous with the peak of ages observed in detrital zircons from the Apollo 12, 14, 15, and 17 landing sites (4340 ± 20 Ma), a Mg-suite Sm-Nd whole rock isochron defined by samples from Apollo 14, 15, 16, and 17 landing sites (4348 ± 25 Ma), and a Ferroan Anorthosite-suite Sm-Nd whole rock isochron defined by samples from the Apollo 15 and 16 landing sites (4354 ± 29 Ma). This implies that Ferroan Anorthosite-suite magmatism is temporally distinct and earlier than magmatism associated with the Mg-suite and the Felsic/Alkali-suite, as predicted by the lunar magma ocean model of lunar differentiation. The short 35 ± 10 Ma interval between primary ferroan anorthosite magmatism and secondary magmatism suggests that the lunar crust formed over a limited period of time. Although heat from decay of long-lived isotopes, large impacts, tidal heating associated with interactions between the Earth and Moon, and density driven overturn of the magma ocean have all been invoked to explain production of ancient secondary crustal magmatism, only tidal heating and cumulate overturn are consistent with the apparent short duration of secondary crustal magmatism and the great depth of crystallization implied for some Mg-suite samples. The initial ε 143 Nd values derived from the 15405 Clast B and 78238 Mg-suite norite isochrons, as well as a Mg-suite whole rock isochron are −0.23 ± 0.11, −0.27 ± 0.74, and −0.25 ± 0.09, respectively. They are identical within uncertainty indicating that Mg-suite and Felsic/Alkali-suite magmas were derived from materials that had the same time averaged Sm/Nd ratios since the formation of the solar system. This, combined with the contemporaneous nature of 15405 Clast B and 78235/6/8/55/56 Mg-suite norite, is consistent with evolution of both samples, and likely both magma suites, from a common source through closed system fractional crystallization or partial melting processes.

Chronology↗

DOE OSTI2025

Apatite geochemistry as a tool for understanding the petrogenesis of layered mafic-ultramafic rocks in the Bushveld Complex, South Africa

The sources of the magmas that formed the Rustenburg Layered Suite of the Bushveld Complex in South Africa remain debated, despite decades of research. Vertical and lateral variation in bulk rock and mineral separate Sr-Nd isotopic compositions, which generally indicate enriched sources, demonstrate that the layered sequence was formed by the emplacement of multiple batches of magma, crucially resulting in episodes of PGE-Cr-V mineralisation. The Lu-Hf isotope compositions of zircon are, however, at odds with the bulk rock Sr-Nd isotopic heterogeneity as they show near homogeneous compositions throughout the layered sequence (εHf (2.06 Ga) =−8). This lack of variation in Hf isotope composition has been attributed to deep, continental lithospheric mantle-related and/or crustal contamination of plume-derived Bushveld magmas. In this study, we analysed the major, trace element and Sr-Nd isotope geochemistry of apatite in the Rustenburg Layered Suite. Apatite occurs as an intercumulus mineral in the lowermost regions and a cumulus mineral in the uppermost regions of the layered sequence and can therefore be used to test existing models for the isotopic disequilibrium between bulk rock Sr-Nd and zircon Hf isotopic compositions. Apatite is largely chlorapatite in the lowermost regions and fluorapatite in the uppermost regions of the layered sequence. The Merensky Reef is unusual in that it comprises both chlorapatite and fluorapatite. Apatite throughout the layered sequence is generally unzoned and shows no evidence of late-stage alteration. Trace element data show that apatite is enriched in L/HREE, with common negative Eu-Sr anomalies. These trace element signatures are consistent with a magmatic origin for the apatite grains, with prior, or concurrent, plagioclase crystallization from the same melt. Variability in in situ Sr and Nd isotope compositions of apatite is recorded throughout the layered sequence with εNd (2.06 Ga) compositions varying between −2.5 and − 10.2 and initial 87 Sr/ 86 Sr compositions varying between 0.7079 and 0.7103 (for the Marikana dikes only). The variability in Sr-Nd isotope compositions of apatite is consistent with the bulk rock (and mineral separate) variation in Sr-Nd isotope compositions, suggesting apatite preserves primary magmatic compositions in the Rustenburg Layered Suite.

Apatite↗

DOE OSTI2024

Modeling nuclear waste disposal in crystalline rocks at the Forsmark and Olkiluoto repository sites – Evaluation of potential thermal–mechanical damage to repository excavations

We conduct coupled thermo-hydro-mechanical modeling of a KBS-3V repository design in crystalline rocks, using data and conditions from the Forsmark in Olkiluoto repository sites in Sweden and Finland. The study focuses on repository performance related to the impact of thermal and hydraulic evolution on the potential for thermal–mechanical damage to underground repository excavations. For the designs and conditions considered at the Forsmark and Olkiluoto repository sites, the simulations show a peak temperature well under the adopted performance target of a 100°C maximum temperature, whereas there is still a high potential for thermal–mechanical damage to the KBS-3V waste deposition holes. The thermal–mechanical damage is much more likely if rock permeability is so low that it delays saturation and swelling of bentonite-clay-based backfill beyond the time for the thermal–mechanical peak, which occurs 50 to 100 years after nuclear waste deposition. We also found that sidewalls of the KBS-3V emplacement tunnels are vulnerable to tensile fracturing due to the combined effect of thermal stressing and backfill swelling. The study highlights a strong interaction between bentonite-based backfill and host rock through capillary suction along with induced rock desaturation. A careful design and selection of the bentonite-clay-based backfill materials for KBS-3V tunnels and deposition holes can facilitate a timely saturation and backfill swelling that in turn can minimize thermal–mechanical damage.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Mobility of Radionuclides in Fractured Carbonate Rocks: Lessons from a Field-Scale Transport Experiment

Current research on radionuclide disposal is mostly conducted in granite, clay, saltstone, or volcanic tuff formations. These rock types are not always available to host a geological repository in every nuclear waste-generating country, but carbonate rocks may serve as a potential alternative. To assess their feasibility, a forced gradient cross-borehole tracer experiment was conducted in a saturated fractured chalk formation. The mobility of stable Sr and Cs (as analogs for their radioactive counterparts), Ce (an actinide analog), Re (a Tc analog), bentonite particles, and fluorescent dye tracers through the flow path was analyzed. The migration of each of these radionuclide analogs (RAs) was shown to be dependent upon their chemical speciation in solution, their interactions with bentonite, and their sorption potential to the chalk rock matrix. The brackish groundwater resulted in flocculation and immobilization of most particulate RAs. Nevertheless, the high permeability of the fracture system allowed for fast overall transport times of all aqueous RAs investigated. This study suggests that the geochemical properties of carbonate rocks may provide suitable conditions for certain types of radionuclide storage (in particular, brackish, high-porosity, and low-permeability chalks). Nevertheless, careful consideration should be given to high-permeability fracture networks that may result in high radionuclide mobility.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗