Search NASA⌕ Search

SEARCH · Search NASA

Results for “shape estimation”

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 OSTI2026

ResStock Measure Documentation: Residential Two-Stage Geothermal Heat Pump (4.0 COP, 20.5 EER)

The goal of this work is to develop energy efficiency, demand flexibility, and other retrofit end-use load shapes (electricity, gas, propane, or fuel oil) that cover a majority of the high-impact, market-ready (or nearly market-ready) measures. "Measures" refers to retrofits that can be applied to buildings during modeling. An "end-use savings shape" is the difference in energy consumption between a baseline building and a building with an energy efficiency, demand flexibility, or other retrofit measure applied. It results in a time-series profile that is broken down by end use and fuel (electricity or on-site gas, propane, or fuel oil use) at each time step. ResStock is a highly granular, physics-based, bottom-up model that uses multiple data sources, statistical sampling methods, and advanced building energy simulations to estimate the annual subhourly energy consumption of the residential building stock across the United States. The baseline model intends to represent the U.S. residential building stock as it existed in 2018. Technical documentation for the inputs and assumptions in the baseline building stock model is available in Reyna et al. (2025). Calibration and validation of the baseline model results are available in the final technical report of the End-Use Load Profiles project (Wilson et al. 2022). This document focuses on a single end-use savings shape measure: Residential Two-Stage Geothermal Heat Pump (4.0 COP, 20.5 EER). This document builds on details established in the single-stage document (Maguire et al. 2025) to detail differences in the approach to modeling this higher efficiency, but more commonly deployed, type of geothermal heat pump. Specific EnergyPlus objects and product specific curves used are highlighted along with showing the results of this measure compared to the baseline and single-speed geothermal heat pumps. Two-speed geothermal heat pumps are able to save even more energy and on utility bills than single-speed products, albeit at the expense of a higher first cost.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

ResStock Measure Documentation: Residential Variable-Speed Geothermal Heat Pump (4.4 COP, 30.9 EER)

The goal of this work is to develop energy efficiency, demand flexibility, and other retrofit end-use load shapes (electricity, gas, propane, or fuel oil) that cover a majority of the high-impact, market-ready (or nearly market-ready) measures. "Measures" refers to retrofits that can be applied to buildings during modeling. An "end-use savings shape" is the difference in energy consumption between a baseline building and a building with an energy efficiency, demand flexibility, or other retrofit measure applied. It results in a time-series profile that is broken down by end use and fuel (electricity or on-site gas, propane, or fuel oil use) at each time step. ResStock (TM) is a highly granular, physics-based, bottom-up model that uses multiple data sources, statistical sampling methods, and advanced building energy simulations to estimate the annual subhourly energy consumption of the residential building stock across the United States. The baseline model intends to represent the U.S. residential building stock as it existed in 2018. Technical documentation for the inputs and assumptions in the baseline building stock model is available in Reyna et al. (2025). Calibration and validation of the baseline model results are available in the final technical report of the End-Use Load Profiles project (Wilson et al. 2022). This documentation focuses on a single end-use savings shape measure: Residential Variable-Speed Geothermal Heat Pump (GHP). This document provides the relevant new modeling information for variable-speed systems not previously covered in either the single-stage or two-stage documents. Variable-speed GHPs represent the most efficient option available for this technology: They provide the most savings, with up to 46% for the applicable portion of the housing stock, compared to 31% for less efficient single-stage GHPs. Additional results shown here detail how the savings change for sections of the housing stock with different base heating fuels and in different climate zones, and they show the savings potential by state for both heating and cooling. Utility bills and electric panel impacts are also shown and discussed.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Cosmological constraints from a joint DESI DR1 Full-Shape and DR2 BAO

We present a cosmological analysis combining full-shape (FS) clustering measurements from the Dark Energy Spectroscopic Instrument (DESI) DR1 with baryon acoustic oscillation (BAO) measurements from DESI DR2. To achieve a robust combination that accounts for the correlation between the two data releases, we employ the ShapeFit compression method and estimate the joint covariance using EZmocks. This compressed approach inherently mitigates the prior volume effects that have previously dominated Bayesian constraints from DESI data with minimal external priors. Consequently, we obtain — for the first time within a Bayesian framework — reliable DESI-only constraints on extensions to ΛCDM using only a Big Bang Nucleosynthesis prior on the baryon density and a wide prior on the spectral index. In flat ΛCDM, we find Ω m = 0.3035 ± 0.0085, h = 0.6876 ± 0.0059, and σ 8 = 0.822 ± 0.034. For the w 0 w a CDM dynamical dark energy model, we measure w 0 = -0.49 ± 0.25 and w a = -1.52 ± 0.77, improving constraints by ∼ 30% relative to the analogous DR1 measurement and reducing the discrepancy with ΛCDM to 1.4σ when compared to BAO only analyses. We also report competitive limits on the sum of neutrino masses and spatial curvature. This work demonstrates that the ShapeFit compression provides a prior-robust and computationally efficient pathway to constrain beyond-ΛCDM physics with large-scale structure.

baryon acoustic oscillations↗
DOE OSTI2025

A novel approach for large-scale wind energy potential assessment

Increasing wind energy generation is central to grid decarbonization, yet methods to estimate wind energy potential are not standardized, leading to inconsistencies and even skewed results. This study aims to improve the fidelity of wind energy potential estimates through an approach that integrates geospatial analysis and machine learning (i.e., Gaussian process regression). We demonstrate this approach to assess the spatial distribution of wind energy capacity potential in the Contiguous United States (CONUS). We find that the capacity-based power density ranges from 1.70 MW/km2 (25th percentile) to 3.88 MW/km2 (75th percentile) for existing wind farms in the CONUS. The value is lower in agricultural areas (2.73 ± 0.02 MW/km2, mean ± 95 % confidence interval) and higher in other land cover types (3.30 ± 0.03 MW/km2). Notably, advancements in turbine manufacturing could reduce power density in areas with lower wind speeds by adopting low specific-power turbines, but improve power density in areas with higher wind speeds (>8.35 m/s at 120m above the ground), highlighting opportunities for repowering existing wind farms. Wind energy potential is shaped by wind resource quality and is regionally characterized by land cover and physical conditions, revealing significant capacity potential in the Great Plains and Upper Texas. The results indicate that areas previously identified as hot spots using existing approaches (e.g., the west of the Rocky Mountains) may have a limited capacity potential due to low wind resource quality. Improvements in methodology and capacity potential estimates in this study could serve as a new basis for future energy systems analysis and planning.

Dai, Tao↗
DOE OSTI2025

Elucidating the Competitive Hydrodeoxygenation of Lignin-Derived Oxygenates over Bulk MoO 3 Catalyst through Kinetic Analysis

Ambient pressure hydrodeoxygenation (HDO) of lignin-derived oxygenates over molybdenum oxide-based catalysts is an effective strategy to produce chemicals that can be directly integrated into our existing petrochemical infrastructure. Complexities pertaining to the simultaneous kinetic and mechanistic analysis of HDO have limited research endeavors to single-compound systems. Although valuable insight into the catalytic reaction has been gained through this approach, it provides limited understanding of the competitive adsorption behavior manifest in a realistic multioxygenate reaction environment. To address this shortcoming, simultaneous gas-phase acetone and anisole HDO was performed at 330 °C and ≤1 bar H 2 partial pressure over bulk MoO 3 . Propene, propane, and benzene were the HDO products formed, showing a similar product distribution to the single-compound system. Selectivity to propene and propane was ∼14 times higher than benzene, even at three times higher anisole partial pressure compared to acetone. A negative anisole HDO (−0.97 ± 0.22) rate order with varying acetone partial pressure suggested a strong inhibition effect on anisole HDO by acetone. Conversely, with increasing anisole partial pressure, a rate order of −0.07 ± 0.12 was observed for acetone HDO, implying a weak impact of anisole cofeed on acetone HDO. A kinetic-driven approach was taken to estimate the relative adsorption constants of the oxygenates. Acetone exhibited a 6.4 times higher adsorption propensity on the HDO active site than anisole. Relative adsorption constants for phenolics increased with increasing basicity of the oxygenate but decreased for aliphatic molecules, suggesting a volcano-shaped relationship. The results suggest the possibility of an optimal electron density around the molecule’s oxygen atom to maximize the molecule’s adsorption strength.

acid sites↗
DOE OSTI2024

Multiscale porosity characterization in additively manufactured polymer nanocomposites using micro-computed tomography

Extrusion-based additive manufacturing (AM) of polymer composites exhibits complex thermally driven phenomena that introduce severe discontinuities in the internal structure across length scales, especially voids or porosity. This study utilizes a high-throughput porosity characterization technique to analyze large datasets from numerous micro-computed tomography (mu CT) scans to capture the influence of AM print parameters on the size, shape, and location of porosity across multiple print layers (up to a few cm) on fused granular fabrication (FGF) printers. The materials investigated include nanocomposite formulations based on commercially relevant nylon-12 and polyether ketone ketone (PEKK) materials comprising nano- or micro- sized fillers. The estimated global porosity follows an inverse linear correlation against the bulk density of the printed samples. Increasing the extrusion multiplier (EM) and the nozzle temperature while decreasing the print speeds reduces the global porosity. Outlier analyses (local porosity morphology) show that faster print speeds and higher extrusion rates result in long, slender inter-layer voids, while lower nozzle temperatures lead to large, symmetrical, inter-bead voids (at the bead junction). Lack of active chamber temperature increases inter-layer and intra-bead voids with a two-fold increase in global porosity. Overall, the micro filler-reinforced composites exhibit higher global porosity than nanofiller-reinforced composites, which is attributed to the increased mismatch in the thermal expansion coefficient between the filler and the polymers used in the study.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Heliostat sizing methodology for concentrating solar thermal industrial process heat projects

This study presents a method to obtain a heliostat size that minimizes the levelized cost of heat (LCOH) of a heliostat-based concentrating solar thermal system for applications of solar heating for industrial processes at operating temperatures from 565 to 1550°C. The method extends prior work by embedding a routine for system design that obtains near-optimal subsystem sizes, increasing the fidelity of drive cost functions, and adding an optical performance model to supplement the previously developed cost models, which we update to reflect current pricing trends. An illustrative business case is developed for Daggett, California, targeting specified annual thermal energy outputs of 50 to 400 GWh th . Optical performance is modeled using verified estimates from the literature. A surrogate heliostat cost model, derived from commercial heliostat designs and scaled for production volume, installation, and operations and maintenance costs, is used to develop cost functions. Results show that heliostat size strongly affects the LCOH, producing a characteristic U-shaped trend with a robust near-optimal window of 7-20 m 2 ; the heliostat size producing the lowest project cost in our study grows slightly as the project size increases, and is reduced as the operating temperature increases. The findings in this study are consistent with the general trend of smaller heliostats being deployed at existing projects for high-temperature industrial process heat and reflect the significant reduction in power electronics and other per-heliostat costs. The methodology we propose is general and can be tailored to revised cost curves as the technology continues to evolve.

14 SOLAR ENERGY↗
DOE OSTI2024

Coercivity of (Fe 0.7 Co 0.3 ) 2 B Nanowire and Its Bonded Magnet

(Fe 0.7 Co 0.3 ) 2 B are potential permanent magnets material due to its large saturation magnetization and high Curie temperature. However, it has moderate magnetocrystalline anisotropy (MCA) and low coercivity. One way to improve its coercivity is to combine the contributions from magnetocrystalline- and magnetic-shape anisotropy by preparing (Fe 0.7 Co 0.3 ) 2 B nanowires. We study the effects of size, morphology, and surface defects on the hard magnetic properties of nanowires using micromagnetic simulation. The hard magnetic properties of (Fe 0.7 Co 0.3 ) 2 B nanowire-bonded magnets are estimated, including the role of inter-wire magnetostatic interaction. By considering the existence of local reductions in MCA energy of up to 30% on the surface layer of nanowires, the anisotropic bonded magnet with a 65% vol. of (Fe 0.7 Co 0.3 ) 2 B nanowires would have typical remanence, B r = 7.6–8.4 kG, coercivity, H ci = 9.6–9.9 kOe, and maximum energy product, (BH) m = 14–17.8 MGOe. Developing effective technology for synthesizing nanowires and fabricating corresponding bonded magnets is promising for manufacturing practical magnets based on the magnetic phase with a relatively low or moderate MCA, such as (Fe 0.7 Co 0.3 ) 2 B.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dimorphos’s Material Properties and Estimates of Crater Size from the DART Impact

On 2022 September 26, the Double Asteroid Redirection Test (DART) spacecraft intentionally collided with Dimorphos, the moon of the binary asteroid system 65803 Didymos. This collision provided the first full-scale test of a kinetic impactor for planetary defense. Images from DART’s DRACO camera revealed Dimorphos to be an oblate spheroid covered in boulders of varying sizes and shapes. Very little was known about Dimorphos prior to DART’s impact, including its shape, structure, and material properties. Approach observations and those following the DART impact have provided crucial knowledge that narrows the parameter space relevant to modeling the impact into Dimorphos. Here we present the results of a suite of hydrocode simulations of the DART impact on Dimorphos. Despite remaining uncertainties, initial models of DART’s kinetic impact provide important information about the results of DART (e.g., potential crater size and morphology, ejecta mass) and the properties of Dimorphos. Simulations here suggest that Dimorphos has near-surface strength ranging from a few Pascals to tens of kPa, which corresponds to crater sizes of ~40–60 m. Simulated crater sizes provide a crucial comparison metric for the European Space Agency Hera mission when it arrives at the Didymos system. Hera’s measurement of crater size in combination with measurement of Dimorphos’s mass will allow us to assess our simulations and provide the information needed to make the DART impact experiment both the first test of a planetary defense mitigation mission and the first full-scale planetary defense simulation validation exercise.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Search for new physics in final states with semivisible jets or anomalous signatures using the ATLAS detector

A search is presented for hadronic signatures of beyond the Standard Model (BSM) physics, with an emphasis on signatures of a strongly coupled hidden dark sector accessed via resonant production of a 𝑍′ mediator. The ATLAS experiment dataset collected at the Large Hadron Collider from 2015 to 2018 is used, consisting of proton-proton collisions at $\sqrt{𝑠}$ = 13 TeV and corresponding to an integrated luminosity of 140 fb −1 . The 𝑍′ mediator is considered to decay to two dark quarks, which each hadronize and decay to showers containing both dark and Standard Model particles, producing a topology of interacting and noninteracting particles within a jet known as “semivisible.” Machine learning methods are used to select these dark showers and reject the dominant background of mismeasured multijet events, including an anomaly detection approach to preserve broad sensitivity to a variety of BSM topologies. A resonance search is performed by fitting the transverse mass spectrum based on a functional form background estimation. No significant excess over the expected background is observed. Results are presented as limits on the production cross section of semivisible jet signals, parametrized by the fraction of invisible particles in the decay and the 𝑍′ mass, and by quantifying the significance of any generic Gaussian-shaped mass peak in the anomaly region.

particle dark matter↗
DOE OSTI2025

Financial Analysis of the High Flow Experiment conducted at the Glen Canyon Dam during Water Year 2023

The Glen Canyon Dam (GCD) is a Colorado River Storage Project (CRSP) power resource that is a component of the Salt Lake City Area Integrated Projects (SLCA/IP). The 2016 record of decision (ROD) for the GCD long-term experimental and management plan (LTEMP) final Environmental Impact Statement (EIS) specified criteria for GCD monthly water releases, daily and hourly operating limits, and experimental releases. This report examines the financial implications of the high flow experiment (HFE) conducted at GCD during the spring of Water Year (WY) 2023 as required by the LTEMP HFE Protocol. This report is part of a series of reports that describe the financial costs of LTEMP experimental releases since the 2016 ROD was adopted in January 2017. Previous reports analyzed the impact of several past HFEs and Bug Flow Experiments. This report focuses on the HFE conducted in April 2023. For this experimental release, financial costs of approximately $1.33 million were incurred because the HFE required sustained water releases exceeding the power plant’s maximum turbine flow rate. In addition, during the experiment, operators were not allowed to shape GCD power production, either to follow Firm Electric Service (FES) customer day-ahead energy deliveries or to respond to market prices. This study identifies the main factors contributing to the HFE costs and examines the interdependencies among these factors. It applies an integrated set of tools to estimate Western Area Power Administration (WAPA) financial impacts by simulating GCD under two types of cases; namely, (1) a “With Experiment” case that mimics the operations that actually occurred and (2) a “Without Experiment” case that simulates operations under the assumption that the HFE did not occur. The “With Experiment” case mimics operations during the HFE and the entire month the HFE occurred. It complies with LTEMP hourly and daily operating criteria. The “Without Experiment” case assumes that the HFE did not occur. The monthly water release volume is assumed to be identical under both cases. The Colorado River Storage Project Python-based model (CRiSPPy) model was the main modeling tool used to simulate the dispatch of the GCD hydropower plant and associated water releases from Lake Powell. In the modeling process, the research team used extensive data sets and historical information on SLCA/IP power plant characteristics, hydrologic conditions, and WAPA’s power purchases and sales prices. In addition to estimating the financial impact of the HFE, the team used the CRiSPPy model to gain insights into the interplay among ROD operating criteria, exceptions made to criteria to accommodate the HFE, and WAPA operating practices.

13 HYDRO ENERGY↗
DOE OSTI2026

Heliostat Sizing Methodology for Solar Heat for Industrial Processes

This study presents a method to obtain a heliostat size that minimizes the levelized cost of a heliostat-based concentrating solar thermal system for industrial process heat (IPH) applications at operating temperatures from 565 to 1550 degrees Celsius. The method extends prior work by embedding a routine for system design that obtains near-optimal subsystem sizes, increasing the fidelity of drive cost functions, and adding an optical performance model. An illustrative business case is developed for Daggett, California, targeting specified annual thermal energy outputs of 50 to 400 GWhth. Optical performance is modeled using verified estimates from the literature. A surrogate heliostat cost model, derived from commercial heliostat designs and scaled for production volume, installation, and operations and maintenance costs, is used to develop cost functions. Results show that heliostat size strongly affects the levelized cost of heat (LCOH), producing a characteristic U-shaped trend with a robust near-optimal window of 8 - 12 m2; the heliostat size producing the lowest project cost in our study grows slightly as the project size increases, and is reduced as the operating temperature increases. The findings in this study are consistent with the general trend of smaller heliostats under deployment at existing projects for high-temperature industrial process heat and reflect the significant reduction in power electronics and other per-heliostat costs. The methodology we propose is general and can be tailored to revised cost curves as the technology continues to evolve.

14 SOLAR ENERGY↗
DOE OSTI2026

Brightest Cluster Galaxy ellipticity as proxy for halo shape: Orientation bias, assembly bias, and potential selection effects in SZ-selected clusters

The orientation of triaxial galaxy clusters with respect to the line-of-sight is expected to be one of the prime sources of scatter and potential bias in optical observables (e.g., richness and weak-lensing signal) of galaxy clusters. In this work, we use the observed shape of the central Brightest Cluster Galaxy (BCG) as proxy for the orientation along the line-of-sight for clusters selected via the Sunyaev-Zel'dovich (SZ) effect from the South Pole Telescope (SPT) and Atacama Cosmology Telescope (ACT) surveys, matched to optically selected clusters from the Dark Energy Survey Year 3 (DES). We construct two samples of clusters that are designed to be identical in SZ mass estimate and redshift but with the roundest vs. the most elliptical BCGs, which we expect to correspond to BCGs (and clusters) with major axes aligned along the line-of-sight vs. in the plane of the sky, respectively. We find that the optical richness of round-BCG clusters is $\sim 10$% larger than that of elliptical-BCG clusters, in agreement with the expectation from projection effects and presenting the first such detection in data. The density profiles, however, are not in agreement with the expectation from projection effects: the 1-halo term (below $6~h^{-1}\rm{Mpc}$) of both the weak-lensing and galaxy density profiles are the same for the subsamples, contrary to previous studies based on X-ray selected clusters. In the 2-halo regime (above $6~h^{-1}\rm{Mpc}$), we find a significant excess of the elliptical-BCG cluster profiles compared to the round-BCG cluster profiles, which is the opposite of the expectation from numerical simulations. We hypothesize that the intrinsic shape of the BCG reflects not just the orientation angle, but also intrinsic properties of the cluster which can affect both the SZ signal and the amplitude of the 2-halo term.

Srinivasan, Radhakrishnan [SUNY, Stony Brook] (ORC↗
DOE OSTI2026

Upper ocean mixing, surface heat fluxes, and heat content variability in the upper 150 m during Hurricane Laura (2020)

Observations (from thermistor sensors and a hurricane glider) and model-based temperature were used to answer questions including the following: 1) What physical and environmental conditions explain the rapid intensification of Hurricane Laura? 2) How did the pre-existing warm mixed layer at Stone mooring (StM) influence the degree of cooling in the mixed layer? 3) How did vertical mixing, surface heat fluxes, and advective processes collectively shape the mixed layer heat evolution at StM? Observations measured by thermistor sensors showed that Hurricane Laura induced a cooling of -1.2oC in the mixed layer when Hurricane Laura’s wind speed increased to 47 m/s on 26 August. This contrasts with the -1.04oC mixed layer temperature change estimated from model-based temperature. At StM, the presence of a 31oC warm mixed layer and elevated heat content (60-80kJ/cm2) effectively preconditioned the upper ocean ahead of Hurricane Laura’s passage. After Hurricane Laura passed, turbulent mixing associated with the storm transported the heat anomalies downward into subsurface layers. In addition, this study shows that a surface heat flux of 5.04 kJ/cm2 supported the intensity of Hurricane Laura as it traversed over the StM. Using a mixed layer heat balance model, this study shows that entrainment flux and surface flux mostly contributed to the observed change in the mixed layer temperature.

Tsai, Senam↗
DOE OSTI2024

Self-diffusion of liquid deuterium hydride and liquid tritium

Here, we present a quasi-elastic neutron scattering study of liquid deuterium hydride carried out using the Disk Chopper Spectrometer at the National Institute of Standards and Technology. Under saturated vapor pressure, the self-diffusion constant of deuterium hydride obeys an Arrhenius law D = D 0 exp(-E A /k B T), where the prefactor D 0 is given by D 0 = 9.5 ± 1.2 Å 2 /ps and the activation energy is given by E A = 58 ± 2 K. We apply the quantum law of corresponding states to the known diffusion constants of the hydrogen isotopologues. From this application, we estimate that D 0 ≈ 9.1 Å 2 /ps and E A ≈ 75 K in liquid tritium. Young’s theory of quantum-mechanical effects in van der Waals fluids is shown to apply to the diffusion constants of the liquid hydrogens. Our results underscore the importance of nuclear quantum effects in shaping the properties and behavior of the hydrogen isotopologues.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Phylogenetic and ecological drivers of the avian lung mycobiome and its potentially pathogenic component

Vertebrate lungs contain diverse microbial communities, but little is known about the drivers of community composition or consequences for health. Microbiome assembly by processes such as dispersal, coevolution, and host-switching can be probed with comparative surveys; however, few studies exist for lung microbiomes, particularly for the fungal component, the mycobiome. Distinguishing among fungal taxa that are generalist or specialist symbionts, potential pathogens, or incidentally inhaled spores is urgent because of potential for emerging diseases. Here, we characterize the avian lung mycobiome and test the relative influences of environment, phylogeny, and functional traits. We used metabarcoding and culturing from 195 lung samples representing 32 bird species across 20 families. We identified 526 fungal taxa as estimated by distinct sequence types (zOTUs) including many opportunistic pathogens. These were predominantly from the phylum Ascomycota (79%) followed by Basidiomycota (16%) and Mucoromycota (5%). Yeast and yeast-like taxa (Malassezia, Filobasidium, Saccharomyces, Meyerozyma, and Aureobasidium) and filamentous fungi (Cladosporium, Alternaria, Neurospora, Fusarium, and Aspergillus) were abundant. Lung mycobiomes were strongly shaped by environmental exposure, and further modulated by host identity, traits, and phylogenetic affinities. Our results implicate migratory bird species as potential vectors for long-distance dispersal of opportunistically pathogenic fungi.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Strain-Based Assessment of Shale Caprock during Cyclic Underground Hydrogen Storage

Successful large-scale underground hydrogen storage (UHS) in depleted gas reservoirs depends on the integrity of the overlying caprock to prevent hydrogen loss during cyclic injection and depletion. Prior studies on crushed Marcellus shale, a potential caprock, indicate that cyclic hydrogen injection and depletion induces microstructural changes, increasing porosity and permeability. However, the extent of these changes in intact shale remains unclear. This study presents a strain-based experimental approach to quantify volumetric strain evolution in intact Marcellus shale matrix under unconstrained stress conditions. A quadrant-shaped shale sample without visible fractures underwent eight hydrostatic pore-pressure cycles (injection to 1500 psi and depletion to 500 psi in 250 psi steps). Linear strain gauges measured strain in three orthogonal directions. Results indicate progressive plastic strain accumulation, leading to an ∼12% increase in matrix porosity after eight cycles, with an estimated 19% increase after 30 cycles. This porosity increase follows a logarithmic trend, suggesting a diminishing effect in later cycles. Additionally, permeability and diffusive mass flux are projected to rise by ∼70% over 30 cycles, enhancing hydrogen migration risk. The shale matrix also exhibited mechanical stiffening over successive cycles, limiting large-scale deformation but not preventing porosity enhancement. A new parameter, α, was introduced to characterize shale sensitivity to cyclic loading, aiding UHS caprock assessments. These findings underscore the necessity of incorporating cyclic loading effects in UHS site selection and operational strategies to ensure long-term storage integrity.

08 HYDROGEN↗
DOE OSTI2024

Terrestrial photosynthesis inferred from plant carbonyl sulfide uptake

Terrestrial photosynthesis, or gross primary production (GPP), is the largest carbon flux in the biosphere, but its global magnitude and spatiotemporal dynamics remain uncertain. The global annual mean GPP is historically thought to be around 120 PgC yr −1 , which is about 30–50 PgC yr −1 lower than GPP inferred from the oxygen-18 ( 18 O) isotope and soil respiration. This disparity is a source of uncertainty in predicting climate–carbon cycle feedbacks. Here, in this study, we infer GPP from carbonyl sulfide, an innovative tracer for CO 2 diffusion from ambient air to leaf chloroplasts through stomata and mesophyll layers. We demonstrate that explicitly representing mesophyll diffusion is important for accurately quantifying the spatiotemporal dynamics of carbonyl sulfide uptake by plants. From the estimate of carbonyl sulfide uptake by plants, we infer a global contemporary GPP of 157 (±8.5) PgC yr −1 , which is consistent with estimates from 18 O (150–175 PgC yr −1 ) and soil respiration ($149^{+29}_{-23}$) PgC yr −1 ), but with an improved confidence level. Our global GPP is higher than satellite optical observation-driven estimates (120–140 PgC yr –1 ) that are used for Earth system model benchmarking. This difference predominantly occurs in the pan-tropical rainforests and is corroborated by ground measurements, suggesting a more productive tropics than satellite-based GPP products indicated. As GPP is a primary determinant of terrestrial carbon sinks and may shape climate trajectories, our findings lay a physiological foundation on which the understanding and prediction of carbon–climate feedbacks can be advanced.

Lai, Jiameng↗