Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Subsurface”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Adsorption Properties of Au−Ni Surface Alloys with a Nonstoichiometric Moiré Structure: A Density Functional Theory Study

Due to the large lattice mismatch between gold and nickel, gold–nickel surface alloys can form unique nonstoichiometric overlayer structures characterized by a moiré pattern and subsurface defects. For this work, we performed density functional theory (DFT) calculations to study the adsorption of molecular oxygen, atomic hydrogen, and atomic carbon on a gold–nickel(111) surface alloy with 0.46 monolayer gold randomly distributed in the surface layer. We observed six distinct adsorption structures for molecular oxygen characterized by intramolecular stretching frequencies of <700, 729, 795, 857, 929, and 1004 cm –1 , which describe well the experimentally observed high-resolution electron energy-loss spectra. Surface atomic hydrogen adsorption is associated with adsorbate–surface modes in the ∼1000 cm –1 range, while subsurface hydrogen can have features as low as ∼400 cm –1 . We observed a unique adsorption structure for atomic carbon inside the surface dislocation loop defect, which explains the experimentally observed low carbon-surface mode at ∼340 cm –1 . Our study sheds light on the unique adsorption properties of the gold–nickel surface alloys and helps with rationalizing vibrational frequency experimental studies for this system.

adsorption↗
DOE OSTI2025

Direct Deoxygenation of Phenol over Fe-Based Bimetallic Surfaces Using On-the-Fly Surrogate Models

We present an accelerated nudged elastic band (NEB) study of phenol direct deoxygenation (DDO) on Fe-based bimetallic surfaces using a recently developed Gaussian process regression (GPR) calculator. Our test calculations demonstrate that the GPR calculator achieves up to 3 times speedup compared to conventional density functional theory calculations while maintaining high accuracy, with energy barrier errors below 0.015 eV. Using GPR-NEB, we systematically examine the DDO mechanism on pure Fe(110) and surfaces modified with Co and Ni in both top and subsurface layers. Our results show that subsurface Co and Ni substitutions preserve favorable thermodynamics and kinetics for both C–O bond cleavage and C–H bond formation, comparable to those on the pure Fe(110) surface. In contrast, top-layer substitutions generally increase the C–O bond cleavage barrier, render the step endothermic, and result in significantly higher reverse reaction rates, making DDO unfavorable on these surfaces. This work demonstrates the effectiveness of GRR-accelerated transition state searches for complex surface reactions and provides insights into rational design of bimetallic catalysts for selective deoxygenation.

Aromatic compounds↗
DOE OSTI2023

Oxidation and Reduction of Polycrystalline Cerium Oxide Thin Films in Hydrogen

This study investigates the oxidation state of ceria thin films' surface and subsurface under 100 mTorr hydrogen using ambient pressure X-ray photoelectron spectroscopy. We examine the influence of the initial oxidation state and sample temperature (25-450 °C) on the interaction with hydrogen. Our findings reveal that the oxidation state during hydrogen interaction involves a complex interplay between oxidizing hydride formation, reducing thermal reduction, and reducing formation of hydroxyls followed by water desorption. In all studied conditions, the subsurface exhibits a higher degree of oxidation compared to the surface, with a more subtle difference for the reduced sample. The reduced samples are significantly hydroxylated and covered with molecular water at 25 °C. We also investigate the impact of water vapor impurities in hydrogen. We find that although 1 × 10 -6 Torr water vapor oxidizes ceria, it is probably not the primary driver behind the oxidation of reduced ceria in the presence of hydrogen.

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

Degeneracy in Intercalated Pb Phases under Buffer-Layer Graphene on SiC(0001) and Diffuse Moiré Spots in Surface Diffraction

First-principles density functional theory (DFT) is used to analyze the stability of Pb intercalated phases under buffer layer graphene on SiC(0001) as a function of the supercell size, Pb coverage, and degree of Pb ordering. By comparing the chemical potentials of such two-dimensional Pb structures, here we find that there is a family of structurally distinct thermodynamically preferred Pb subsurface configurations with minute stability differences. These differences are comparable to the thermal energies at about 450 °C, where the Pb intercalated phases are grown. High-resolution surface-diffraction experiments using Spot Profile Analysis Low-Energy Electron Diffraction (SPA-LEED) confirm this high degree of degeneracy of the Pb intercalated phases from broad, low-intensity moiré spots observed exclusively from intercalated Pb. The low intensity of the moiré spots implies the coexistence of structurally different subsurface Pb phases.

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

Dynamic Surface Restructuring in Cu(Au) Alloys Driven by Oxygen-Mediated Au Mobility

Alloying plays a crucial role in tuning the surface properties of metals, but the atomic-level mechanisms by which alloying elements influence surface structure dynamics under reactive conditions remain elusive. Using Cu(Au) in oxidizing environments as a model system, we reveal a dynamic oxygen-induced transformation of the topmost atomic layer into a periodically hill-and-valley morphology with reversible switching between undulated and flattened surface states. These interconversions are driven by the retreat of surface Au to the subsurface during oxygen adsorption and its resegregation to the surface upon oxygen desorption. This cyclical mobility establishes a feedback loop, allowing the surface to dynamically reconfigure in response to changes in the oxygen pressure. The results offer a broadly applicable framework for understanding atomic-scale surface restructuring in alloy systems, where differences in the chemical reactivity of alloying elements drive dynamic redistribution between surface and subsurface regions. As a result, this dynamic coupling has practical implications for designing corrosion-resistant coatings and metastable nanostructures with tunable catalytic properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Elucidating the Initial Steps in α-Uranium Hydriding Using First-Principles Calculations

Hydrogen embrittlement of uranium, which arises due to the formation of a structurally weak pyrophoric hydride, poses a major safety risk in material applications. Previous experiments have shown that hydriding begins on the top or near the surface (i.e., subsurface) of α-uranium. However, the fundamental molecular-level mechanism of this process remains unknown. Here, in this work, starting from pristine α-U bulk and surfaces, we present a systematic investigation of possible mechanisms for the formation of metal hydride. Specifically, we address this problem by examining the individual steps of hydrogen embrittlement, including surface adsorption, subsurface absorption, and the interlayer diffusion of atomic hydrogen. Furthermore, by examining these processes across different facets, we highlight the importance of both (1) hydrogen monolayer coverage and (2) applied tensile strain on hydriding kinetics. Taken together, by studying previously overlooked phenomena, this study provides foundational insights into the initial steps of this overall complex process. We anticipate that this work will guide near-term future development of multiscale kinetic models for uranium hydriding and subsequently identify potential strategies to mitigate this undesired process.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Hybrid Sorbents for 129 I Capture from Contaminated Groundwater

Radioiodine ( 129 I) poses a risk to the environment due to its long half-life, toxicity, and mobility. It is found at the U.S. Department of Energy Hanford Site due to legacy releases of nuclear wastes to the subsurface where 129 I is predominantly present as iodate (IO 3 – ). To date, a cost-effective and scalable cleanup technology for 129 I has not been identified, with hydraulic containment implemented as the remedial approach. Here, novel high-performing sorbents for 129 I remediation with the capacity to reduce 129 I concentrations to or below the US Environmental Protection Agency (EPA) drinking water standard and procedures to deploy them in an ex-situ pump and treat (P&T) system are introduced. Furthermore, this includes implementation of hybridized polyacrylonitrile (PAN) beads for ex-situ remediation of IO 3 – -contaminated groundwater for the first time. Iron (Fe) oxyhydroxide and bismuth (Bi) oxyhydroxide sorbents were deployed on silica substrates or encapsulated in porous PAN beads. In addition, Fe–, cerium (Ce)–, and Bi–oxyhydroxides were encapsulated with anion-exchange resins. The PAN–bismuth oxyhydroxide and PAN–ferrihydrite composites along with Fe- and Ce-based hybrid anion-exchange resins performed well in batch sorption experiments with distribution coefficients for IO 3 – of >1000 mL/g and rapid removal kinetics. Of the tested materials, the Ce-based hybrid anion-exchange resin was the most efficient for removal of IO 3 – from Hanford groundwater in a column system, with 50% breakthrough occurring at 324 pore volumes. The functional amine groups on the parent resin and amount of active sorbent in the resin can be customized to improve the iodine loading capacity. These results highlight the potential for IO 3 – remediation by hybrid sorbents and represent a benchmark for the implementation of commercially available materials to meet EPA standards for cleanup of 129 I in a large-scale P&T system.

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

Porous Colloidal Nanoparticles as Injectable Multimodal Contrast Agents for Enhanced Geophysical Sensing

Injecting fluids into underground geologic structures is crucial for the development of long-term strategies for managing captured carbon and facilitating sustainable energy extraction operations. Here, we have previously reported that the injection of metal–organic frameworks (MOFs) into the subsurface can enhance seismic monitoring tools to track fluids and map complex structures, reduce risk, and verify containment in carbon storage reservoirs because of their absorption capacity of low-frequency seismic waves. Here, we demonstrate that water-based Cr/Zn/Zr MOF colloidal suspensions (nanofluids) are multimodal geophysical contrast agents that enhance near-wellbore logging tools. Based on experimental fluid-only measurements, MIL-101(Cr), ZIF-8, and UiO-66 nanofluids have distinct complex conductivity and/or low-field nuclear magnetic resonance (NMR) signatures that are relevant to field-deployed technologies, implying the potential to enhance near-wellbore monitoring of CO 2 injection and associated processes with downhole logging tools. Small- and wide-angle X-ray scattering characterization of ~0.5 wt % MIL-101(Cr) suspensions confirmed phase stability and provided insight into the fractal nature of colloidal nanoparticles. Finally, low-field (2 MHz) NMR measurements of MIL-101(Cr) nanofluid injection into a prototypical Berea sandstone demonstrate how paramagnetic high-surface area MOFs may dominate the relaxation times of hydrogen-bearing fluids in porous geologic matrices, enhancing the mapping of near-surface and near-wellbore transport pathways and advancing sustainable subsurface energy technologies.

NMR↗
DOE OSTI2023

High-Sensitivity Low-Energy Ion Spectroscopy with Sub-Nanometer Depth Resolution Reveals Oxidation Resistance of MoS 2 Increases with Film Density and Shear-Induced Nanostructural Modifications of the Surface

For decades, density has been attributed as a critical aspect of the structure of sputter-deposited nanocrystalline molybdenum disulfide (MoS 2 ) coatings impacting oxidation resistance and wear resistance. Despite its importance, there are few examples in the literature that explicitly investigate the relationship between the density and oxidation behaviors of MoS 2 coatings. Aging and oxidation are primary considerations for the use of MoS 2 coatings in aerospace applications as they inevitably experience prolonged storage in water and oxygen-rich environments prior to use. Oxidation that is either limited to the first few nanometers or through the bulk of the coating can result in seizure due to high initial coefficients of friction or component failure from excessive wear. High-sensitivity low-energy ion spectroscopy (HS-LEIS) and Rutherford backscattering spectrometry (RBS) are both used to understand the extent of oxidation throughout the first ~10 nanometers of the surface of pure sputtered nanocrystalline MoS 2 coatings after high-temperature aging and how it is impacted by the density of coatings as measured by RBS. Results show that low-density coatings (ρ = 3.55 g/cm 3 ) exhibit a more columnar microstructure and voiding, which act as pathways for oxidative species to penetrate and interact with edge sites, causing severe surface and subsurface oxidation. Furthermore, HS-LEIS of surfaces sheared prior to oxidation reveals that the oxidation resistance of low-density MoS 2 coatings can be significantly improved by shear-induced reorientation of the surface microstructure to a basal orientation and elimination of pathways for oxygen into the bulk through compaction of surface and subsurface voids.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Accelerated Cu2O Reduction by Single Pt Atoms at the Metal-Oxide Interface§

The reducibility of metal oxides, when serving as the catalyst support or the active sites themselves, plays an important role in heterogeneous catalytic reactions. Here we present an integrated experimental and theoretical study that reveals how the addition of small amounts of atomically dispersed Pt at the metal/oxide interface dramatically enhances the reducibility of a Cu2O thin film by H2. X-ray photoelectron spectroscopy (XPS) and temperature-programmed desorption (TPD) results reveal that upon oxidation, a PtCu single-atom alloy (SAA) surface is covered by a thin Cu2O film and is therefore unable to dissociate H2. Despite this, in situ studies using ambient pressure (AP) XPS reveal that the presence of a small amount of Pt under the oxide layer can, at the single-atom limit, promote the reduction of Cu2O by H2 at room temperature. We build two density functional theory based surface models to better understand these experimental findings: A Cu2O/Cu(111)-like surface oxide layer, known as the ‘29’ oxide, in which Pt is alloyed into the Cu(111) surface as well as a PtCu SAA. Our calculations suggest that the increased activity is due to the presence of atomically dispersed Pt under the surface oxide layer, which weakens the Cu-O bonds in its immediate vicinity, thus making the interface between subsurface Pt and the surface oxide a nucleation site for the formation of metallic Cu. This initial step in the reduction process results in the presence of surface Pt atoms surrounded by metallic Cu patches, and the Pt atoms become active in H2 dissociation, which consequently accelerates the reduction of the oxide layer. This work demonstrates how isolated Pt atoms at the metal/oxide interface of a Cu-based catalyst accelerate the reduction of the oxide and, therefore, helps maintain the active, reduced state of the catalyst under reaction conditions, as is required for the synthesis of methanol from CO2 hydrogenation.

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

Adsorption of REEs to Kaolinite via Ion Exchange and Surface Complexation as a Function of Water Chemistry

Rare earth elements (REEs) are critical components of modern technology behind renewable energy, transportation, and electronics but have a limited current supply. A substantial portion of global REE production relies on ion adsorption deposits. A high abundance of kaolinite in REE enrichment zones within these deposits suggests that kaolinite controls the subsurface migration of REEs. This study aimed to improve the current understanding of REE binding to kaolinite under varying water chemistry conditions. We conducted batch experiments with kaolinite (KGa-2) and three REEs (Nd, Dy, and Yb) at varying pH, electrolyte concentration, dissolved inorganic carbon (DIC), low molecular weight organic acids (citric and oxalic acids), and total REE concentration conditions. Increasing electrolyte concentration inhibits REE adsorption at pH < 7, suggesting that ion exchange contributes to adsorption at these pH values. DIC affects adsorption above pH 7–8 by forming strong aqueous complexes with heavy REEs. Citric acid decreases REE adsorption via aqueous complexation of REEs at pH > 5 but does not affect adsorption at pH < 5. The surface complexation model captures the main adsorption trends with two mechanisms: ion exchange on basal planes at pH < ∼6 and inner-sphere surface complexation to edge sites at pH > ∼6. Equilibrium constants for surface complexation increase in the order of Yb > Dy > Nd, indicating a higher strength of adsorption for heavy REEs. This study demonstrates how water chemistry conditions control the adsorption mechanisms that may determine the mobility of REEs in subsurface environments rich in kaolinite.

58 GEOSCIENCES↗
DOE OSTI2025

Spatial Distribution and Clustering of Glycosaminoglycans in Electrospun Gelatin-Based Scaffolds

The extracellular matrix (ECM) is comprised of components like collagen, elastin, and glycosaminoglycans (GAGs). Electrospun fibrous scaffolds are designed to replicate the form and composition of the native ECM, often requiring blending of various ECM component mimics to enhance cellular responses. However, the spatial distribution of blended components within these fibers remains unclear. This study investigates the spatial distribution of chondroitin sulfate-C (CSC) in electrospun gelatin-based scaffolds. scanning electron microscopy (SEM), attenuated reflectance-Fourier transform infrared (ATR-FTIR) spectroscopy, X-ray photoelectron spectroscopy (XPS), and Time-of-flight Secondary Ion Mass Spectrometry (ToF-SIMS) were applied for surface and subsurface chemical characterization of the fibrous scaffolds. SEM confirmed a fibrous morphology, while ATR-FTIR and XPS analyses indicated the presence of CSC through the identification of sulfate groups. ToF-SIMS imaging, alongside K-means clustering and Ripley’s K function, revealed a nonuniform CSC distribution with higher concentrations at the top layer of the scaffold. This study demonstrates that CSC presentation at the fiber surface varies with depth and differs from bulk incorporation while reveals nanoscale clustering and spatial heterogeneity at both the surface and subsurface of electrospun gelatin fibers. These findings define an underexplored design consideration with potential to influence cell–scaffold interactions.

Animal derived food↗
DOE OSTI2026

Integration of the Biot–Gassmann Fluid Substitution Method and Machine Learning-Based Velocity–Stress Relationship for Estimating In Situ Stresses

Recent advancements have shown that in situ stresses can be reliably estimated through an integrated machine/deep learning (ML/DL)-based framework, which relies on models trained and validated using true triaxial ultrasonic velocity (TUV) experimental data that involve measurements of ultrasonic velocity in saturated rocks under varying stress configurations. However, when the goal is to interpret lower frequency measurements, it may be more appropriate to run experiments on dry rocks and then obtain Biot–Gassmann-derived equivalent saturated velocities (low-frequency approximation) and employ these quantities for training ML/DL models to predict in situ stress. Whether the dispersion effect of frequency on the velocity–stress relationship substantially impacts in situ stress prediction is an important and unresolved question. This work presents an enhancement of ML/DL-based workflow by training and implementing ML/DL models using equivalent saturated acoustic velocities (low-frequency) obtained by applying Biot–Gassmann fluid substitution on the ultrasonic velocities of dry cores. The models were trained on TUV data sets derived from three subsurface cores extracted from the geothermal well 16B(78)-32 at the Utah FORGE site. Each core was subjected to 75 unique stress configurations for velocity measurement in the dry state. The ML/DL trained on the TUV data set with equivalent saturated velocities demonstrated promising performance to predict in situ stress in subsurface geological rocks using velocity–stress relationships with R 2 of 0.86, 0.971, and 0.975 and root mean squared error (RMSE) of 2.59, 1.92, and 1.80 for validation/testing phases of vertical, minimum horizontal, and maximum horizontal stress models, respectively. Additionally, interpretation and explanation by Shapley additive explanations (SHAP) analysis further improved scientific validation and model reliability for estimating in situ stresses.

colloids↗
DOE OSTI2021

Close Observation of Hydraulic Fracturing at EGS Collab Experiment 1: Fracture Trajectory, Microseismic Interpretations, and the Role of Natural Fractures

Abstract Despite the wide application of hydraulic fracturing in enhanced geothermal system (EGS) development and unconventional hydrocarbon production, our understanding of the many factors affecting the propagation of hydraulic fractures has relied on circumstantial evidence. This is partly due to the scarcity of direct observations in the subsurface. The EGS Collab project attempts to address these issues in the context of EGS research by performing intermediate‐scale (∼10 m) hydraulic stimulation experiments in a thoroughly characterized and heavily instrumented underground testbed. This paper analyzes the data collected from the first suite of hydraulic fracturing tests in this testbed, consisting of seven stimulation episodes. High‐quality microseismic data delineate five planar features very clearly. Combining fracture‐wellbore intersection observations from distributed temperature sensing (DTS) and visual observations from an open‐hole well, as well as prior in‐situ stress measurements, we conclude with high certainty that the four larger planes were hydraulic fractures. The growth of the hydraulic fractures was temporarily halted by a prominent, open natural fracture in the testbed but they eventually crossed and slightly reoriented under continued stimulation. Mineral‐filled (i.e., healed) natural fractures, though prevalent in the testbed, did not have perceivable effects on hydraulic fracture propagation. The high‐quality, mutually corroborating data sets allowed conclusions to be drawn with high confidence and attests to the advantage of intermediate‐scale experiments in subsurface research.

58 GEOSCIENCES↗
DOE OSTI2021

Validation of the Community Land Model Version 5 over the Contiguous United States (CONUS) using in situ and remote sensing data sets

The Community Land Model (CLM) is an effective tool to simulate the biophysical and biogeochemical processes and their interactions with the atmosphere. Although CLM Version 5 (CLM5) constitutes various updates in these processes, its performance in simulating energy, water and carbon cycles over the Contiguous United States (CONUS) at scales which land surface changes and hydrometeorological and hydroclimatological applications are more locally relevant is yet to be assessed. In this study, we conducted three simulations at 0.125? during 1979-2018 over the CONUS using different configurations of CLM, namely CLM5-biogeochemistry (CLM5BGC), CLM4.5BGC, and CLM5-satellite phenology (CLM5SP). We validated and compared their simulations against multiple remote-sensed and in-situ datasets. Overall, the parametric and structural updates (e.g., carbon cost for nitrogen uptake, variable soil thickness, dry surface layer) in CLM5 improve its ability in capturing terrestrial biogeochemical dynamics. The low evapotranspiration in CLM5BGC is associated with biases in simulating vegetation phenological characteristics rather than soil water limitations. The mismatch between CLM5BGC-simulated peak leaf area index and reference data can be attributed to CLM5BGC's inability in simulating phenology of trees and grasses. The differences between CLM-simulated irrigation and reference estimates can be attributed to differences between processes represented in models and in reality, and uncertainties in input and validation datasets. Evaluation against observations at small catchments suggest that hydrologic parameters needed to be calibrated to improve simulations of runoff, especially subsurface runoff. Additional efforts are needed to incorporate spatially-distributed plant phenology and physiology parameters and regional-specific agricultural management practices (e.g., planting, harvest).

Cheng, Yanyan↗
DOE OSTI2022

On an Extensive Late Hydrologic Event in Gale Crater as Indicated by Water‐Rich Fracture Halos

Abstract We analyze spatially pervasive, light‐toned “halos” associated with fractures in a sedimentary unit (Stimson) of Gale crater, Mars, and report a similar network of halos discovered in a separate geologic group (Bradbury). Through a dedicated active neutron measurement campaign, we provide independent confirmation of the water‐rich nature of these features. Together with mineralogical and geochemical data, these features are consistent with abundant hydrated amorphous silica (opal‐A). We suggest that the mineral and amorphous assemblages are indicative of formation under low‐temperature and predominantly low‐pH conditions (passive silica enrichment) with minor contribution of silica (active silica enrichment) from adjacent units. We show that there is significant amorphous silica in the array of sedimentary rocks in Gale crater, allowing them to play a role in an active silica enrichment phase of halo formation. We suggest that the involved alteration event was short lived and our finding of vast halo networks in a distant, older unit implies a more vast network of hydrologic subsurface conduits than previously known. This relatively recent subsurface hydrologic system was present long after the transition from a warm and wet to a cold and dry Martian environment, extending the habitability conditions on Mars to an epoch that is generally considered not favorable for life on the surface. Finally, our bulk H quantification of these features, which ranges from ∼3–6 wt% H 2 O‐equivalent‐H, suggests that the amorphous material in halos hosts ample supplies of readily released water, making them a considerable resource at the otherwise dry Martian equator.

Gabriel, Travis S. J.↗
DOE OSTI2021

Seismic Ambient Noise Analyses Reveal Changing Temperature and Water Signals to 10s of Meters Depth in the Critical Zone

The critical zone sustains terrestrial life, but we have few tools to explore it efficiently beyond the first few meters of the subsurface. Using analyses of high-frequency ambient seismic noise from densely spaced seismometers deployed in the forested Shale Hills subcatchment of the Susquehanna Shale Hills Critical Zone Observatory (SSHCZO), we show that temporal changes in seismic velocities at depths from ~1 m to tens of m can be detected. Furthermore, these changes are driven by variations at the land surface. The Moving-Window Cross-Spectral (MWCS) method was employed to measure seismic-velocity changes in coda waves at hourly resolution in 10 different frequency bands. We observed a diurnal signal, a seasonal signal, and a meteorological-event-based signal. These signals were compared to time-series measurements of precipitation, well water levels, soil moisture, soil temperature, air temperature, latent heat flux, and air pressure in the heavily instrumented catchment. Most of the velocity changes can be explained by variations in temperature that result in thermoelastic strains that propagate to depth. But some double minima in seismic velocity time-series observed after large rain events were attributed in part to the effects of water infiltration. These results show that high-frequency ambient noise data may in some locations be used to detect changes in the critical zone from ~1 to ~100 m or greater depth with hourly resolution. But interpretation of such data requires multiple environmental data sets to deconvolve the complex interrelationships among thermoelastic and hydrological effects in the subsurface critical zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Agreement of Analytical and Simulation-Based Estimates of the Required Land Depth in Climate Models

Previous analytical and simulation-based analyses suggest that deeper land surface models are needed to realistically simulate the terrestrial thermal state in climate models, with implications for land-atmosphere interactions. Analytical approaches mainly focused on the subsurface propagation of harmonics such as the annual temperature signal, and a direct comparison with climate-change model output has been elusive. Here, this study addresses the propagation of a harmonic pulse fitted to represent the timescale and amplitude of anthropogenic warming. Its comparison to land model simulations with stepwise increased bottom boundary depth leads to an agreement between the simulation-based and analytical frameworks for long-term climate trends. Any depth increase gradually decreases the relative error in the subsurface thermodynamics, and a minimum depth of 170 m is recommended to simulate the ground climate adequately. The approach provides an accurate estimate of the required land-model depth for climate-change simulations and assesses the relative bias in insufficiently deep land models.

54 ENVIRONMENTAL SCIENCES↗