Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geologic uncertainty”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Modeling the distribution of the endangered Jemez Mountains salamander (Plethodon neomexicanus) in relation to geology, topography, and climate

The Jemez Mountains salamander (Plethodon neomexicanus; hereafter JMS) is an endangered salamander restricted to the Jemez Mountains in north-central New Mexico, United States. This strictly terrestrial and lungless species requires moist surface conditions for activities such as mating and foraging. Threats to its current habitat include fire suppression and ensuing severe fires, changes in forest composition, habitat fragmentation, and climate change. Forest composition changes resulting from reduced fire frequency and increased tree density suggest that its current aboveground habitat does not mirror its historically successful habitat regime. However, because of its limited habitat area and underground behavior, we hypothesized that geology and topography might play a significant role in the current distribution of the salamander. We modeled the distribution of the JMS using a machine learning algorithm to assess how geology, topography, and climate variables influence its distribution. The best habitat suitability model indicates that geology type and maximum winter temperature (November to March) were most important in predicting the distribution of the salamander (23.5% and 50.3% permutation importance, respectively). Minimum winter temperature was also an important variable (21.4%), suggesting this also plays a role in salamander habitat. Our habitat suitability map reveals low uncertainty in model predictions, and we found slight discrepancies between the designated critical habitat and the most suitable areas for the JMS. Because geological features are important to its distribution, we recommend that geological and topographical data are considered, both during survey design and in the description of localities of JMS records once detected.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Optimization of the Post-Operational Phase on Two Belgian Multi-Unit Nuclear Power Plants: the Case of the Non-Fissile Irradiated Core Items - 20156

The current legal framework in Belgium foresees the progressive phase out of nuclear power between October 2022 (Doel 3) and December 2025 (Doel 2). Upon its definitive shutdown, each unit of the Tihange and Doel sites will enter a Post-Operational Phase (POP) and be prepared for its Decontamination and Decommissioning (D and D). Prior to obtaining the D and D license, the Operator Electrabel is legally required to remove any non-fissile irradiated core items stored in the deactivation pools. The non-fissile irradiated core items consist essentially of control rods, poison rods and source thimbles as well as thimble plugs and foreign materials irradiated during operation: - Their significant content in highly radiant radionuclides (up to 6 TBq of Co-60 per kg of irradiated material) renders all existing operational waste management processes inadequate due to insufficient biological shielding; - Their high concentrations in long-lived radionuclides call for their disposal in a geological repository for which no final design nor waste acceptance criteria are expected prior to 2050. Uncertainties in the Belgian energy supply and security, however, require the Operator to be prepared for a partial nuclear phase out, where one or more units would benefit from lifetime extension while the remaining units would undergo decommissioning. The present paper aims at presenting how Electrabel, in partnership with Tractebel, addressed this challenge by maximizing the use of synergies within the respective sites as well as between both sites themselves, all the while accounting for site specificities. The most recent results and state of progress of the project will be detailed and the first lessons learned will be shared. The project has been split in multiple tasks and phased as follows: - An inventory phase aimed at mapping the contents, origin, composition and history of the non-fissile irradiated core items; - A pre-characterization phase based on neutron activation models; - A waste sorting phase aimed at separating waste forms for which an evacuation route exists from those for which such route does not exist; - A feasibility phase aimed at exploring all possible scenarios for the management of non-fissile irradiated core items and identifying the optimal feasible solution for each site; - A preparation phase (currently ongoing), developing further the optimal solution and ensuring that back-up solutions are available for any foreseeable change of context (licensing issue, modification in the nuclear phase-out program, etc.) and initiating early contacts with potential subcontractors for segmentation works and cask manufacturers, as well as the Belgian regulatory body and waste management agency. This phase also foresees the investigation of destructive and non-destructive radiological measurements to support the detailed characterization of the waste forms; - A realization phase (future work). (authors)

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

De-risking superhot EGS development through 3D play fairway analysis: Methodology development and application at Newberry Volcano, Oregon, USA

The DErisking Exploration for multiple geothermal Plays in magmatic ENvironments (DEEPEN) project aims to de-risk exploration of geothermal plays in magmatic systems, with a focus on superhot and supercritical geothermal systems. This article presents advancements in geoscientific assessment of superhot enhanced geothermal system (EGS) resources as part of the DEEPEN project, particularly in exploration and identification of these systems. While there exists a general consensus on key components of conventional hydrothermal systems, there has not been a prior consensus on such components for EGS, let alone superhot EGS. The DEEPEN project identified these key components for superhot EGS plays, and used them to modify the traditional geothermal play fairway analysis (PFA) approach to identify favorable areas in superhot EGS plays in 3D. This was done through modeling both favorability and uncertainty separately for each key component of a superhot EGS play, incorporating diverse 3D geoscientific datasets, including geologic features, models produced from direct observations, single inversions, and joint inversions. The PFA is applied to Newberry Volcano as a form of validation. Finally, the results appear to validate the methodology, aligning with expectations derived from conceptual modeling, and highlighting the area targeted for EGS stimulation well NWG 55-29 as favorable, and suggesting potential additional areas worthy of exploration at Newberry.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Multi‐Layered Systems for Permanent Geologic Storage of CO 2 at the Gigatonne Scale

Abstract The effectiveness of Carbon Capture and Storage (CCS) as an imperative decarbonization technology relies on the sealing capacity of a fine‐grained caprock to permanently store CO 2 deep underground. Uncertainties in assessing the caprock sealing capacity increase with the spatial and temporal scales and may delay CCS deployment at the gigatonne scale. We have developed a computationally efficient transport model to capture the dynamics of basin‐wide upward CO 2 migration in a multi‐layered setting over geological time scales. We find that massive capillary breakthrough and viscous flow of CO 2 , even through pervasively fractured caprocks, are unlikely to occur and compromise the storage security. Potential leakage from the injection reservoir is hampered by repetitive layering of overlying caprocks. This finding agrees with geologic intuition and should be understandable by the public, contributing to the development of climate policies around this technology with increased confidence that CO 2 will be indefinitely contained in the subsurface.

58 GEOSCIENCES↗
DOE OSTI2024

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

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

Stochastic Modeling Workflow to Generate Representative Geologic Variability in Training Dataset for SMART Initiative

The poster discusses the modeling workflow to generate ensemble of geologic realizations of the Illinois Basin Decatur Project (IBDP) site, based on available site characterization data and inherent uncertainty of those data, for use by project collaborators in DOE SMART Initiative (Phase 2) to build their forward modeling, history matching, and optimization workflows. This poster is summarized from the technical report for the SMART project submitted to U.S. DOE earlier this year.

Ganesh, Priya Ravi↗
DOE OSTI2025

Overview of the SCEC/USGS Community Stress Drop Validation Study Using the 2019 Ridgecrest Earthquake Sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress-drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude-spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress-drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two-week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. Here, in this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite-fault modeling, ground-motion-based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between-method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress-drop values.

58 GEOSCIENCES↗
DOE OSTI2024

Probabilistic Programming for Transportable Source Characterization and Uncertainty Quantification of the North Korean Nuclear Tests 2006–2017

Here, we introduce a transportable technique to determine the yield and depth of burial (DOB) from seismic source spectra of underground nuclear explosions. We demonstrate this technique on the six declared North Korean nuclear tests. This approach derives source spectra in absolute units from regional phase (Pg) amplitudes by correcting the observations for geometric spreading, attenuation, and site amplification. We couple the source spectra and explosion source models with a probabilistic programming framework that integrates deep learning techniques and Bayesian modeling. This approach permits the exchange of information across various data categories to quantify both the data and model uncertainty. This technique stands out as an innovative use of broad‐area propagation models, making it transportable across various geologic settings. This method proves to be effective in scenarios with diverse and/or limited observational data, even when the source depth is unknown. We present new independent estimates of absolute yield and DOB that are consistent with the prior assessments, underscoring the potential of this method in enhancing transportable nuclear explosion monitoring capabilities.

58 GEOSCIENCES↗
NASA NTRS1995

Evidence for Regional Basin Formation in Early Post-Tessera Venus History: Geology of the Lavinia Planitia Area (V55)

On Venus, global topography shows the presence of highs and lows including regional highly deformed plateaus (tesserae), broad rifted volcanic rises, linear lows flanking uplands, and more equidimensional lowlands (e.g. Lavinia and Atalanta planitiae) Each of these terrain types on Venus has relatively distinctive characteristics, but origins are uncertain in terms of mode of formation, time of formation, and potential evolutionary links. There is a high level of uncertainty about the formation and evolution of lowlands on Venus. We have undertaken the mapping of a specific lowlands region of Venus to address several of these major questions. Using geologic mapping we have tried to establish: What is the sequence of events in the formation and evolution of large-scale equidimensional basins on Venus? When do the compressional features typical of basin interiors occur? What is the total volume of lava that occurs in the basins and is this similar to other non-basin areas? How much subsidence and downwarping has occurred after the last major plains units? WE have undertaken an analysis of the geology of the V55 Lavinia Planitia quadrangle in order to address many of these issues and we report on the results here.

Head, J. W.↗
DOE OSTI2022

Investigating the Physical Drivers for the Increasing Tropical Cyclone Rainfall Hazard in the United States

In this study, we investigate both the changes of tropical cyclone (TC) rainfall hazard in the United States under climate change and the relative importance of the factors that cause the changes. We find that under the SSP5 8.5 scenario, the 100-year TC rainfall level can increase by up to 320% along the U.S. coastline by the end of this century. The influence of TC rainfall-producing ability increase is more significant than the influence of TC frequency increase on the increase of the 100-year TC rainfall level (up to 180% vs. 60% increase). Among the different physical drivers for the increase in storm rainfall-producing ability, the increase of TC intensity is the leading factor, followed by changes in TC duration and atmospheric temperature. The projected increase of TC rainfall hazard is robust against the uncertainty in the TC frequency projection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Temperature Uncertainty Modeling with Proxy Structural Data as Geostatistical Constraints for Well Siting: An Example Applied to Granite Springs Valley, NV, USA

Utilizing existing temperature and structural information around Granite Springs Valley, Nevada, we build 3D stochastic temperature models with the aim of evaluating the 3D uncertainty of temperature and choosing between candidate exploration well locations . The data used to support the modeling are measured temperatures and structural proxies from 3D geologic modeling, the latter considered "secondary" data. Two stochastic geostatistical techniques are explored for incorporating the structural proxies: cosimulation and local varying mean. With both the cosimulation and local varying mean methods, many equally likely temperature models (i.e., realizations) are produced, from which temperature probability profiles are calculated at candidate well locations. To aid in choosing between the candidate locations, two quantities summarize the temperature probabilities: Vprior and entropy. Vprior quantifies the likelihood for economic temperatures at each candidate location, whereas entropy identifies where new information has the most potential to reduce uncertainty. In general, the cosimulation realizations have smoother spatial structure, and extrapolate high temperatures at candidate locations that are located along the direction of the longest spatial correlation, which are down dip from existing temperature logs. The smooth realizations result in tight temperature probability profiles that are easier to interpret, but they have unrealistic temperature reversals in some locations because the cosimulation technique does not enforce a conductive geothermal gradient as a baseline (i.e., linearly increasing temperature with depth). The local varying mean results produce realizations with more realistic geothermal gradients, with temperatures increasing downward since a depth-temperature relationship is included. However, because they have much noisier spatial nature compared to cosimulation, it is harder to interpret the temperature probability profiles. The different local varying mean results allow the geologist to determine which proxy (e.g., dilation versus distance to fault termination) should be used given the specific geothermal system. In general, Vprior from local varying mean results identify locations that are close to high values for the structural proxies: areas with highe r probabilities for higher temperatures. The entropy results identify where uncertainty is greatest and therefore new drilling information could be most useful. Though these techniques provide useful information, even when applied to areas of sparse data, our comp arison of these two techniques demonstrates the need for new geothermal geostatistics techniques that combine the advantages of these two methods and that are tailored to the spatial uncertainty issues inherent in geothermal exploration.

3D temperature modeling↗
DOE OSTI2025

"Hidden" hydrothermal technical potential & technoeconomics: Revealing permeability & fluids with more data

Historical hydrothermal estimates have largely relied on temperature or heat flow estimates ignoring the need for natural flowing fluids. More accurate hydrothermal estimates require some indication of permeability and fluids that naturally exist in the subsurface. This paper describes a novel approach that includes proxies of permeability and fluids in hydrothermal estimates by leveraging the relatively data-rich Great Basin. Specifically, nameplate capacities (megawatts) of operating geothermal plants, negative (0 megawatt) locations and 48 geophysical and geologic features are used to used in eXtreme Gradient Boosting (XGBoost) regression to make hydrothermal capacity predictions. Additionally, this work inputs the XGBoost-based hydrothermal predictions into the Renewable Energy Potential (reV) model to quantify technical capacity, its uncertainty and techno-economics. Compared to historical hydrothermal estimates, these predictions adhere to the 37 operating geothermal plants and negative locations. We present a method for subsampling the negative sites to bring the labels into balance that uses the geologic domain knowledge to proportionally represent negatives. Overall, the distributions of the hydrothermal technical capacity and the site levelized cost of energy are respectively much tighter, lower and more accurate than the previous estimates for the Great Basin, as they include geological and geophysical surrogates for permeability and fluids. Percentile (50th and 90th, median and high estimate, respectively) models provide bookends for these metrics.

13 HYDRO ENERGY↗
NASA NTRS1998

Plate Motions, Regional Deformation, and Time-Variation of Plate Motions

The significant results obtained with support of this grant include the following: (1) Using VLBI data in combination with other geodetical, geophysical, and geological data to bound the present rotation of the Colorado Plateau, and to evaluate to its implications for the kinematics and seismogenic potential of the western half of the conterminous U.S. (2) Determining realistic estimates of uncertainties for VLBI data and then applying the data and uncertainties to obtain an upper bound on the integral of deformation within the "stable interior" of the North American and other plates and thus to place an upper bound on the seismogenic potential within these regions. (3) Combining VLBI data with other geodetic, geophysical, and geologic data to estimate the motion of coastal California in a frame of reference attached to the Sierra Nevada-Great Valley microplate. This analysis has provided new insights into the kinematic boundary conditions that may control or at least strongly influence the locations of asperities that rupture in great earthquakes along the San Andreas transform system. (4) Determining a global tectonic model from VLBI geodetic data that combines the estimation of plate angular velocities with individual site linear velocities where tectonically appropriate. and (5) Investigation of the some of the outstanding problems defined by the work leading to global plate motion model NUVEL-1. These problems, such as the motion between the Pacific and North American plates and between west Africa and east Africa, are focused on regions where the seismogenic potential may be greater than implied by published plate tectonic models.

Gordon, R. G.↗
DOE OSTI2024

Upscaling Methods Applied to a Fine-Scale Reservoir Model

This study was conducted as part of the Southwest Regional Partnership on Carbon Sequestration (SWP) project to evaluate how upscaling fine-scale simulation models to coarse-scale simulation models impacted the results. The focus was on the Farnsworth Unit (FWU) and its Morrow' B' Sandstone reservoir, specifically the west half of the field. Due to data limitations and the geologic characteristics of the surrounding area, the upscaling was limited to the west half of the FWU rather than a broader basinscale model. The primary aim was to explore how upscaling impacts numerical simulation models, particularly regarding CO 2 -enhanced oil recovery (EOR) and storage capacity predictions. Upscaling was necessary to reduce computational demands when transitioning from high-resolution geological models to coarser grids, as large-scale simulations with finer grids can be computationally prohibitive. This study expands on previous work by the SWP to understand how additional upscaling, applied to already fine-scale numerical simulation models, affects reservoir performance simulations (Ampomah, Balch, & Grigg, 2015). This is key to understanding how loss of resolution can affect coarsescale model results that may be used for large sensitivity analyses, uncertainty quantifications, and training data for machine learning applications.

02 PETROLEUM↗
DOE OSTI2021

Interferometric SAR modelling of near surface data to improve geological model in the Surat Basin, Australia

Here we present a study on geophysical inverse modelling for subsurface structural properties of an unconventional hydrocarbon site that was monitored previously by Interferometric Synthetic Aperture Radar (InSAR) technology for surface deformation. A static three-dimensional geomodel along with extracted property maps replicates the depth of each underlying stratigraphic unit and structural feature with the density of each geological layer. We examine the hypothesis that integration of elastic properties of each formation layer with InSAR observations in a stratified elastic medium will lead to a viscoelastic geophysical inverse problem that can be solved to estimate fractional volume change at the reservoir level. Moreover, we examine synthetic scenarios in which the elastic properties of the formations are perturbed before determining the resulting impact on the rate of surface deformation. The results show that although the slope of underlying formations, their density and depth can define the extent and pattern of a deformation signal, their properties have a marginal impact on volumetric change compared to the dense network of shallow depth Coal Seam Gas (CSG) mining wells. Besides, it is also demonstrated that the inversion of InSAR deformation maps can resolve the uncertainties associated with low-resolution seismic interpretation as well as filling the data gaps within seismic acquisitions. A significant contribution of this investigation to the geological basin modelling involves a) introducing a remote and non-invasive technology such as InSAR to improve geophysical mapping of subsurface structures such as faults in areas with sparse or no reflective seismic information, and b) applying a multi-layer viscoelastic geophysical source model for an unconventional hydrocarbon reservoir such as CSG.

58 GEOSCIENCES↗
NASA NTRS2014

Relationship between LIBS Ablation and Pit Volume for Geologic Samples: Applications for in situ Absolute Geochronology

In planetary sciences, in situ absolute geochronology is a scientific and engineering challenge. Currently, the age of the Martian surface can only be determined by crater density counting. However this method has significant uncertainties and needs to be calibrated with absolute ages. We are developing an instrument to acquire in situ absolute geochronology based on the K‐Ar method. The protocol is based on the laser ablation of a rock by hundreds of laser pulses. Laser Induced Breakdown Spectroscopy (LIBS) gives the potassium content of the ablated material and a mass spectrometer (quadrupole or ion trap) measures the quantity of 40Ar released. In order to accurately measure the quantity of released 40Ar in cases where Ar is an atmospheric constituent (e.g., Mars), the sample is first put into a chamber under high vacuum. The 40Arquantity, the concentration of K and the estimation of the ablated mass are the parameters needed to give the age of the rocks. The main uncertainties with this method are directly linked to the measures of the mass (typically some μg) and of the concentration of K by LIBS (up to 10%). Because the ablated mass is small compared to the mass of the sample, and because material is redeposited onto the sample after ablation, it is not possible to directly measure the ablated mass. Our current protocol measures the ablated volume and estimates the sample density to calculate ablated mass. The precision and accuracy of this method may be improved by using knowledge of the sample's geologic properties to predict its response to laser ablation, i.e., understanding whether natural samples have a predictable relationship between laser energy deposited and resultant ablation volume. In contrast to most previous studies of laser ablation, theoretical equations are not highly applicable. The reasons are numerous, but the most important are: a) geologic rocks are complex, polymineralic materials; b) the conditions of ablation are unusual (for example, variable vacuum pressure), and c) the ablation is made with hundreds of successive laser pulses. In this work, we aim to understand the effects that occur on LIBS spectra when a homogeneous rock or a mineral is ablated under high vacuum. Understanding these effects is important to define best practices for LIBS measurements and may lead to improved measurement (or possibly prediction) of the ablated volume. We will describe our laboratory approach and first results, and discuss its utility for situ absolute geochronology campaigns.

Devismes, D.↗
DOE OSTI2022

Influence of sampling frequency and estimation method on phosphorus load uncertainty in the Western Lake Erie Basin, Ohio, USA

Accurate estimates of nutrient loads are necessary to identify critical source areas and quantify the impact of management practices on pollutant export. Previous studies have investigated nutrient load estimate uncertainty, but they often focus on nutrient loads estimated using an interpolation method for large-scale watersheds with short-term datasets. The study objective was to quantify uncertainty in soluble reactive phosphorus (SRP), total phosphorus (TP), and suspended solids (SS) load estimates from two small (<10 3 km 2 ) agricultural watersheds in the western Lake Erie Basin resulting from different sampling frequencies. Each watershed had high temporal resolution datasets of discharge (15 min) and nutrient concentration (1 to 3 samples per day) collected over a 30-year period (1990–2020). Firstly, SRP, TP, and SS loads were calculated using the high temporal resolution datasets, which was assumed as “true loads”. Secondly, the high temporal concentration data were decomposed to semiweekly, weekly, biweekly, and monthly sampling and annual loads were estimated using four common load estimation methods to assess the effect of sampling frequency and load estimation method on load estimate error. Across the four different methods, the composite method had the lowest relative root mean square and absolute bias, but the rectangular interpolation method was the most precise. Furthermore, even with semiweekly sampling, the composite method resulted in an unacceptable level of precision (average imprecision = 39 %), while the interpolation method resulted in an unacceptable bias (average absolute bias = 16 %). Because neither method could provide acceptable accuracy and precision at the lowest decrease in sampling (e.t. semiweekly sampling), continued daily sampling is recommended in these watersheds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Spatial controls of methane uptake in upland soils across climatic and geological regions in Greenland

Abstract In the Arctic, the spatiotemporal variation of net methane uptake in upland soils depends on unresolved interactive controls between edaphic and microbial factors not yet included in current models, underpinning the uncertainty of upscaling the Arctic methane budget. Here we show that upland soils in Greenland are consistent methane sinks (−1.83 ± 0.19 nmol methane g −1 dw d −1 ) across a N-S (64–83 °N) pedoclimatic transect. We demonstrate that methane oxidizers abundance, soil pH, and available soil copper are important controls on the spatial variation in methane oxidation. We revised a soil biogeochemical model with a high-resolution land classification and meteorological data for Greenland and tested it against our methane uptake measurements. The model simulated well the magnitudes of observed methane uptake but not the spatial variation across all sites. This work provides novel insights into the controls of methane uptake, which are critical for the accuracy of methane budgets.

Environmental Sciences & Ecology↗