Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model initialization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

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

Cavitation model of the inflationary stage of Big Bang

In this study, we propose a model for the initial stage of the development of the universe analogous to cavitation in a liquid in a negative pressure field. It is assumed that at the stage of inflation, multiple breaks of the metric occur with the formation of areas of physical vacuum in which the generation of matter occurs. The proposed model explains the large-scale isotropy of the universe without ultrafast inflationary expansion and the emergence of a large-scale cellular (cluster) structure, as a result of the development of cavitation ruptures of a false vacuum. It is shown that the cavitation model can be considered on par with (or as an alternative to) the generally accepted inflationary multiverse model of the Big Bang.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Potential Signature of Population III Pair-instability Supernova Ejecta in the BLR Gas of the Most Distant Quasar at z = 7.54*

The search for Population III stars has fascinated and eluded astrophysicists for decades. One promising place for capturing evidence of their presence must be high-redshift objects; signatures should be recorded in their characteristic chemical abundances. We deduce the Fe and Mg abundances of the broadline region (BLR) from the intensities of ultraviolet Mg ii and Fe ii emission lines in the near-infrared spectrum of UKIDSS Large Area Survey (ULAS) J1342+0928 at z = 7.54, by advancing our novel flux-to-abundance conversion method developed for quasars up to z ~ 3. We find that the BLR of this quasar is extremely enriched, by a factor of 20 relative to the solar Fe abundance, together with a very low Mg/Fe abundance ratio: [Fe/H] = +1.36 ± 0.19 and [Mg/Fe] =-1.11 ± 0.12, only 700 million years after the Big Bang. We conclude that such an unusual abundance feature cannot be explained by the standard view of chemical evolution that considers only the contributions from canonical supernovae. While there remains uncertainty in the high-mass end of the Population III initial mass function, here we propose that the larger amount of iron in ULAS J1342+0928 was supplied by a pair-instability supernova (PISN) caused by the explosion of a massive Population III star in the high-mass end of the possible range of 150–300 M ⊙ . Chemical evolution models based on initial PISN enrichment well explain the trend in [Mg/Fe]-z all the way from z < 3 to z = 7.54. We predict that stars with very low [Mg/Fe] at all metallicities are hidden in the galaxy, and they will be efficiently discovered by ongoing new-generation photometric surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Actinide Production in the Neutron-rich Ejecta of a Neutron Star Merger

The rapid neutron-capture ("r-") process is responsible for synthesizing many of the heavy elements observed in both the solar system and Galactic metal-poor halo stars. Simulations of r-process nucleosynthesis can reproduce abundances derived from observations with varying success, but so far they fail to account for the observed overenhancement of actinides, present in about 30% of r-process-enhanced stars. In this work, we investigate actinide production in the dynamical ejecta of a neutron star merger (NSM) and explore whether varying levels of neutron-richness can reproduce the actinide boost. We also investigate the sensitivity of actinide production on nuclear physics properties: fission distribution, β-decay, and mass model. For most cases, the actinides are overproduced in our models if the initial conditions are sufficiently neutron-rich for fission cycling. Here, we find that actinide production can be so robust in the dynamical ejecta that an additional lanthanide-rich, actinide-poor component is necessary in order to match observations of actinide-boost stars. We present a simple actinide-dilution model that folds in estimated contributions from two nucleosynthetic sites within a merger event. Our study suggests that while the dynamical ejecta of an NSM are likely production sites for the formation of actinides, a significant contribution from another site or sites (e.g., the NSM accretion disk wind) is required to explain abundances of r-process-enhanced, metal-poor stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Literature Review of Documented Persistent Secondary Uranium Sources at DOE LM Sites

The U.S. Department of Energy Office of Legacy Management (LM) has the responsibility to protect human health and the environment at 101 legacy uranium sites across the nation by implementing long-term surveillance and maintenance requirements. Included in these responsibilities is to provide for the safe and environmentally sound disposal, long-term stabilization, and control of uranium mill tailings in a manner that minimizes or eliminates radiation health hazards to the public at 21 Uranium Mill Tailings Radiation Control Act (UMTRCA) Title I sites; 6 UMTRCA Title II sites; and 9 Comprehensive Environmental Response, Compensation, and Liability Act/Resource Conservation and Recovery Act sites. These sites are former uranium mills (processing sites), or sites with wastes from processing sites encapsulated in engineered disposal cells (disposal sites), or sites where components of atomic weapons were manufactured. During site characterization, groundwater contamination was often discovered and evaluated or modeled for future movement with the inclusion of uranium sorption in the groundwater plume. Much of this original work did identify uranium processing wastes as the original contaminant source. However, more recent work has identified the existence and influence of persistent secondary uranium sources that may prevent uranium in groundwater from attenuating as quickly as initial evaluations or modeling indicated. By understanding these processes, better decisions can be made to choose the appropriate remediation and closure strategies, identify long-term surveillance and maintenance activities, and more accurately estimate life-cycle costs. The purpose of this report is to summarize the literature available on documented persistence secondary uranium sources at LM sites as part of the Applied Studies and Technology Persistent Secondary Contaminant Sources Project.

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

Geometry-aware training of factorized layers in tensor Tucker format

Reducing parameter redundancies in neural network architectures is crucial for achieving feasible computational and memory requirements during train and inference of large networks. Given its easy implementation and flexibility, one promising approach is layer factorization, which reshapes weight tensors into a matrix format and parameterizes it as the product of two rank-r matrices. However, this family of approaches often requires an initial full-model warm-up phase, prior knowledge of a feasible rank, and it is sensitive to parameter initialization.In this work, we introduce a novel approach to train the factors of a Tucker decomposition of the weight tensors. Our training proposal proves to be optimal in locally approximating the original unfactorized dynamics and stable for the initialization. Furthermore, the rank of each mode is dynamically updated during training.We provide a theoretical analysis of the algorithm, showing convergence, approximation and local descent guarantees. The method's performance is further illustrated through a variety of experiments, showing remarkable training compression rates and comparable or even better performance than the full baseline and alternative layer factorization strategies.

Zangrando, Emanuele [Gran Sasso Science Institute ↗
DOE OSTI2020

Asteroseismic modelling of solar-type stars: a deeper look at the treatment of initial helium abundance

ABSTRACT Detailed understanding of stellar physics is essential towards a robust determination of stellar properties (e.g. radius, mass, and age). Among the vital input physics used in the modelling of solar-type stars which remain poorly constrained, is the initial helium abundance. To this end, when constructing stellar model grids, the initial helium abundance is estimated either (i) by using the semi-empirical helium-to-heavy element enrichment ratio, (ΔY/ΔZ), anchored to the standard big bang nucleosynthesis value, or (ii) by setting the initial helium abundance as a free variable. Adopting 35 low-mass, solar-type stars with multiyear Kepler photometry from the asteroseismic ‘LEGACY’ sample, we explore the systematic uncertainties on the inferred stellar parameters (i.e. radius, mass, and age) arising from the treatment of the initial helium abundance in stellar model grids. The stellar masses and radii derived from grids with free initial helium abundance are lower compared to those from grids based on a fixed ΔY/ΔZ ratio. We find the systematic uncertainties on mean density, radius, mass, and age arising from grids which employ a fixed value of ΔY/ΔZ and those with free initial helium abundance to be ∼ 0.9 per cent, ∼ 2 per cent, ∼ 5 per cent, and ∼ 29 per cent, respectively. We report that the systematic uncertainties on the inferred masses and radii arising from the treatment of initial helium abundance in stellar grids lie within the expected accuracy limits of ESA’s PLATO, although this is not the case for the age.

Nsamba, Benard↗
DOE OSTI2025

Demonstrating Hierarchical System Development With the Common Community Physics Package Single‐Column Model: A Case Study Over the Southern Great Plains

This study demonstrates a specific application of the hierarchical system development (HSD) approach to investigate, analyze, and attribute model issues within the Unified Forecast System (UFS), with a focus on process isolation. By evaluating a non‐precipitating, shallow cumulus case at the Atmospheric Radiation Measurement Southern Great Plains site in the UFS global forecast against the observation, the investigation identifies a warmer and deeper daytime convective planetary boundary layer (PBL) and misrepresented nocturnal PBL transition. Hypothesis testing, which employs the Common Community Physics Package (CCPP) single‐column model (SCM) and uses the same physics as the UFS global model, confirms that these issues are attributed to the model physics and initialization. Specifically, misrepresented PBL processes are linked to problematic surface condition and a lack of cloud formation, which may stem from deficiencies in PBL and cloud microphysics parameterizations and their interactions. The UFS initial condition contributes to an earlier, excessively collapsed daytime convective boundary layer and a lack of decoupling between the stable boundary layer and residual layer late in the afternoon. This work introduces an avenue for the community to engage with the application of HSD, along with the CCPP and CCPP SCM, to understand the interplay of model physics, disentangle the roles of model components, as well as facilitate model and forecast improvement.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Evidence of Hexadecapole Deformation in Uranium-238 at the Relativistic Heavy Ion Collider

State-of-the-art hydrodynamic simulations of the quark-gluon plasma are unable to reproduce the elliptic flow of particles observed at the BNL Relativistic Heavy Ion Collider (RHIC) in relativistic 238 U + 238 U collisions when they rely on information obtained from low-energy experiments for the implementation of deformation in the colliding 238 U ions. Here, we show that this is due to an inappropriate treatment of well-deformed nuclei in the modeling of the initial conditions of the quark-gluon plasma. Past studies have identified the deformation of the nuclear surface with that of the nuclear volume, though these are different concepts. In particular, a volume quadrupole moment can be generated by both a surface hexadecapole and a surface quadrupole moment. This feature was so far neglected in the modeling of heavy-ion collisions, and is particularly relevant for nuclei like 238 U, which is both quadrupole deformed and hexadecapole deformed. With rigorous input from Skyrme density functional calculations, we show that correcting for such effects in the implementation of nuclear deformations in hydrodynamic simulations restores agreement with BNL RHIC data. This brings consistency to the results of nuclear experiments across energy scales, and demonstrates the impact of the hexadecapole deformation of 238 U on high-energy collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Milestone 1.2.10: Steady-state H 2 “roll over” point data for aluminum alloys 1100 and 6061

Extended (> 50 years) dry storage is being evaluated by the U.S. Department of Energy (DOE) for the disposition of ~ 18 metric tons of aluminum-clad spent nuclear fuel (ASNF). Transition of the current ASNF inventory into dry storage—using the standard DOE canister—necessitates a rigorous, predictive understanding of the long-term physical and chemical factors that may influence the integrity of the proposed storage canister, including radiolytic molecular hydrogen (H 2 ) generation. Current model predictions employ initial radiolytic yields of H 2 , the values of which change as the cladding’s H 2 -precursor inventory is depleted and H 2 itself becomes progressively more involved in radiolytic and surface dissociation processes. Consequently, the absorbed radiation dose that this steady-state H 2 yield corresponds to is essential for the evaluation and improvement of model predictions. Here, we report our findings on the long-term generation of H 2 from the gamma irradiation (≤ 36 MGy) of corroded AA1100 and AA6061 coupons in helium environments at ambient temperature and ~ 50% RH. Our findings show that AA1100 systems reached steady-state by ~ 36 MGy, while higher doses were necessary for AA6061 systems. This discrepancy was attributed to the AA6061 coupons developing a thicker corrosion layer that led to the trapping of H 2 and its precursors, and potentially additional chemistries, ultimately delaying the depletion of H 2 precursors and the system’s “roll over” point. Further, current model predictions—based on previous AA1100 data—do not show steady-state attainment until above 120 MGy, which is not the case for the AA1100 data collected here. Consequently, the new alloy dependent data presented here are important for the continued improvement of predictive computer models for evaluating the feasibility of extended storage of ASNF in helium backfilled canisters.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

THERMAL STRESS ANALYSIS OF A SPENT NUCLEAR FUEL CANISTER

The potential for chloride induced stress corrosion cracking (CISCC) in spent nuclear fuel dry storage canisters is a current topic of research by the US Department of Energy Spent Fuel and Waste Science and Technology program. One of the important prerequisites for CISCC is a tensile residual stress state. This study utilizes computational models to provide an initial analysis of thermal stresses in a generic spent nuclear fuel canister. A STAR-CCM+ thermal fluid model provides a temperature profile analysis of the canister based on four different ambient temperature conditions. An ANSYS APDL finite element model incorporates the temperature profiles to analyze the thermal stresses in the canister. The calculated thermal stress magnitudes are in the range of 10 MPa to 80 MPa, which could be significant for crack propagation through the canister wall. The maximum thermal stress is not a concern for the structural failure of the canister, but it is high enough that it is a significant contribution to the total stress state of the canister wall, which also includes residual stress from fabrication, internal pressure from helium cover gas, and potential transient mechanical loads, such as earthquakes. This paper describes the initial finite element thermal stress analysis that was completed in 2020, reports the initial findings, and identifies areas where model refinement is still needed to complete this work in the future.

Jensen, Ben J.↗
DOE OSTI2021

A Spatially Explicit Individual-Based Modeling Approach to Evaluate the Cumulative Effects of Wind Energy Development on the Greater Sage-Grouse: Sensitivity Analysis and Validation of Key Parameters

To address the need for an approach to evaluate cumulative ecological impacts of wind energy development, particularly those on critical wildlife habitats and species, the authors developed a proof-of-concept, landscape-based, spatially explicit individual-based modeling (IBM) framework for populations of the greater sage-grouse (Centrocercus urophasianus) in Albany County, Wyoming, based on published and other available information on the life history of the species. This sage-grouse IBM estimates the spatial and temporal movements of sage-grouse based on modeled effects of habitat characteristics and developments on sage-grouse condition, survivorship, and reproduction. After completing the initial version of the model, we conducted a series of sensitivity analyses of key model parameters and validation of lek occurrences and spatial distributions to better understand performance of the model relative to observed values of important life history characteristics. The parameter sensitivity analysis was conducted as a series of one-factor-at-a-time tests in which the values of selected parameters were individually varied to evaluate the corresponding effects on model output. We examined the model’s response to changes in parameter values to gain insight into the overall robustness of the model, identify key driver parameters, and initiate qualitative validation. Five functions selected for sensitivity analyses were (1) the rate of change in body condition in relationship to habitat suitability index (HSI) of occupied habitat (condition change rate), (2) the effect of competition between individuals with overlapping home ranges expressed as a percentage reduction in HSI (competition factor), (3) value of condition at which survivorship probability is equal to zero (zero survivorship condition), (4) the relationship between female body condition and clutch size, and (5) the relationship between female body condition and nest success. We varied the parameter values to be tested for the sensitivity analysis around the original value used in the model (i.e., tested values were higher and lower than the model’s original value) to examine the effect of this variation on the predicted output of the model. Model outputs evaluated in the sensitivity analysis included (1) predicted age-class distribution (as measured by the percentage of immature individuals within the population); (2) predicted sex ratio; (3) predicted age-sex class distribution (percent of population in each of six age-sex classes, i.e., male and female juveniles, yearlings, and adults); (4) lifespan; (5) population size; and (6) spatial distribution. Of the five parameters analyzed, condition change rate, survivorship, clutch size, and nest success did not substantially affect any of the evaluated model outputs. For the predicted age-class distribution, the percentage of immature individuals predicted by the model for these four parameters ranged from 61.9% to 65.0%, which was reasonable with respect to the reference values of 51.4 to 57.8% in published studies. The ratio of females to males ranged from 1.09 to 1.14, which was comparable to the reference values of 1.2 to 3.0 presented in an existing study. The predicted mean lifespan ranged from 1.62 to 1.74 years, which appears to be reasonable with respect to the reference values of 0.9 to 1.1 years for sharp-tail grouse and greater prairie chickens presented by published studies. The predicted population size ranged from 1,877 to 3,140. This prediction would be reasonable with respect to our estimate of 5,000 yearlings and adults, which was based on the USFWS estimate for the state of Wyoming scaled to the number of leks in the county. The predicted spatial distribution was not substantially affected by the parameter values examined. In our sensitivity analyses, the only tested parameter that had a noticeable effect on model output was competition factor. Over the range of values tested for this parameter, predicted population size exhibited the largest range (from 1,833 to 4,737) of all parameters tested. As the competition factor increased, the predicted population size decreased. This reflects the effect competition had on the number of individuals that could co occur in high-quality habitat patches. Despite this apparently significant effect, the range of values appears to be reasonable with respect to our estimate of 5,000 yearlings and adults, which was based on the USFWS estimate for the state of Wyoming scaled to the number of leks in the county.

17 WIND ENERGY↗
DOE OSTI2023

Evaluating skill in predicting the Interdecadal Pacific Oscillation in initialized decadal climate prediction hindcasts in E3SMv1 and CESM1 using two different initialization methods and a small set of start years

Abstract It is a daunting challenge to conduct initialized hindcasts with enough ensemble members and associated start years to form a drifted climatology from which to compute the anomalies necessary to quantify the skill of the hindcasts when compared to observations. This limits the ability to experiment with case studies and other applications where only a few initial years are needed. Here we run a set of hindcasts with CESM1 and E3SMv1 using two different initialization methods for a limited set of start years and use the respective uninitialized free-running historical simulations to form the model climatologies. Since the drifts from the observed initial states in the hindcasts toward the uninitialized model state are large and rapid, after a few years the drifted initialized models approach the uninitialized model climatological errors. Therefore, hindcasts from the limited start years can use the uninitialized climatology to represent the drifted model states after about lead year 3, providing a means to compute forecast anomalies in the absence of a large hindcast sample. There is comparable skill for predicting spatial patterns of multi-year Pacific sea surface temperature anomalies in the domain of the Interdecadal Pacific Oscillation using this method compared to the conventional methodology with a large hindcast data set, though there is a model dependence to the drifts in the two initialization methods.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Automated Model Updating Effect on the Linear and Nonlinear Dynamic Responses of Historical Masonry Structures

Finite element models supply the simulation of the problems to obtain the structural behaviors under different loading conditions. The updating processes based on changing the uncertain parameters such as material properties, boundary conditions, meshing size etc. are made to create actual finite element models in order to represent the real conditions of the problems. This procedure can be made using two methods, namely manual model updating and global/local automatic model updating. The global/local automatic model updating procedures are very popular owing to the applicability to all kind of engineering structures, minimizing the differences nearly to zero, availability of damage localization and structural health monitoring. In this paper a historical masonry armory building selected as an application and its finite element analyses, experimental measurements and automated model updating processes are carried out. The finite element model of the structures was developed with ANSYS software and first three natural frequencies are attained. The experimental dynamic characteristics such as natural frequencies, damping ratios and mode shapes are obtained using ambient vibration method. The differences between experimental and numerical frequencies are minimized up to the range of 0–1%, by using local updating procedure. The linear and nonlinear seismic analyses of the building are conducted on initial and updated models. At the end of the analyses, maximum displacements, principal stresses and strains are given with detail. In addition, crack distributions diagrams are attained to validate the stresses accumulation points. It can be seen that there is an increasing trend in the displacement, stress and strain values after finite element model updating both linear and nonlinear analysis. In addition, linear and nonlinear analysis results are compared with each other to evaluate the nonlinear analysis effect. It is concluded that the displacement and strain values are increased after nonlinear analysis. But, there is not any agreement in stress values.

42 ENGINEERING↗
DOE OSTI2022

Monte Carlo event generator for initial conditions of conserved charges in nuclear geometry

At top collider energies where baryon stopping is negligible, the initial state of heavy-ion collisions is overall charge neutral and predominantly composed of gluons. Nevertheless, there can also be significant local fluctuations of the baryon number, strangeness, and electric charge densities about zero, perturbatively corresponding to the production of quark/antiquark pairs. These previously ignored local charge fluctuations can permit the study of charge diffusion in the quark-gluon plasma (QGP), even at top collider energies. In this paper we present a new model denoted ICCING (initial conserved charges in nuclear geometry) which can reconstruct the initial conditions of conserved charges in the QGP by sampling a (g→q¯q) splitting probability over the initial energy density. We find that the new charge distributions generally differ from the bulk energy density; in particular, the strangeness distribution is significantly more eccentric than standard bulk observables and appears to be associated with the geometry of hot spots in the initial state. The new information provided by these conserved charges opens the door to studying a wealth of new charge- and flavor-dependent correlations in the initial state and ultimately the charge transport parameters of the QGP.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Combustion Characteristics of a Tier II Gasoline Certification Fuel and its Surrogate with Iso-butanol: experiments and detailed numerical modeling

This presentation will describe results of an experimental and numerical study of combustion of a Tier II gasoline certification fuel and its surrogate blended with iso-butanol as a potential biofuel additive. The experiments were carried out under conditions where gas transport originated entirely by fuel evaporation with external convection (both forced and natural) were minimized. The resulting near spherically symmetric droplet flames were amenable to detailed numerical modeling. The droplet initial diameters were in the narrow range of 0.55 mm to 0.62 mm while and burning is performed under standard atmosphere. The combustion characteristics such as evolution of droplet, flame and soot shell diameters for indolene/iso-butanol and surrogate/iso-butanol mixtures are compared. The presence of iso-butanol is experimentally found to reduce the formation of soot as the iso-butanol concentration was increased. Simulations are also reported using a detailed kinetic chemistry comprised of 324 species and 17608 reactions using the CRECK POLIMI mechanism. The detailed numerical model included unsteady gas and liquid transport, soot chemistry and radiative transport. Droplet burn rates were well predicted while flame diameters showed more variability depending on how the flame was defined in the numerical simulation (peak gas temperature and OH concentration). Simulations also show the influence of iso-butanol concentration on flame temperature, greenhouse gas emissions and particulate emissions.

Guo, Songtao↗
DOE OSTI2020

Observations for Model Intercomparison Project (Obs4MIPs): status for CMIP6

The Observations for Model Intercomparison Project (Obs4MIPs) was initiated in 2010 to facilitate the use of observations in climate model evaluation and research, with a particular target being the Coupled Model Intercomparison Project (CMIP), a major initiative of the World Climate Research Programme (WCRP). To this end, Obs4MIPs (1) targets observed variables that can be compared to CMIP model variables; (2) utilizes dataset formatting specifications and metadata requirements closely aligned with CMIP model output; (3) provides brief technical documentation for each dataset, designed for nonexperts and tailored towards relevance for model evaluation, including information on uncertainty, dataset merits, and limitations; and (4) disseminates the data through the Earth System Grid Federation (ESGF) platforms, making the observations searchable and accessible via the same portals as the model output. Taken together, these characteristics of the organization and structure of obs4MIPs should entice a more diverse community of researchers to engage in the comparison of model output with observations and to contribute to a more comprehensive evaluation of the climate models. At present, the number of obs4MIPs datasets has grown to about 80; many are undergoing updates, with another 20 or so in preparation, and more than 100 are proposed and under consideration. A partial list of current global satellite-based datasets includes humidity and temperature profiles; a wide range of cloud and aerosol observations; ocean surface wind, temperature, height, and sea ice fraction; surface and top-of-atmosphere longwave and shortwave radiation; and ozone (O 3 ), methane (CH 4 ), and carbon dioxide (CO 2 ) products. A partial list of proposed products expected to be useful in analyzing CMIP6 results includes the following: alternative products for the above quantities, additional products for ocean surface flux and chlorophyll products, a number of vegetation products (e.g., FAPAR, LAI, burned area fraction), ice sheet mass and height, carbon monoxide (CO), and nitrogen dioxide (NO 2 ). While most existing obs4MIPs datasets consist of monthly-mean gridded data over the global domain, products with higher time resolution (e.g., daily) and/or regional products are now receiving more attention. Along with an increasing number of datasets, obs4MIPs has implemented a number of capability upgrades including (1) an updated obs4MIPs data specifications document that provides additional search facets and generally improves congruence with CMIP6 specifications for model datasets, (2) a set of six easily understood indicators that help guide users as to a dataset's maturity and suitability for application, and (3) an option to supply supplemental information about a dataset beyond what can be found in the standard metadata. With the maturation of the obs4MIPs framework, the dataset inclusion process, and the dataset formatting guidelines and resources, the scope of the observations being considered is expected to grow to include gridded in situ datasets as well as datasets with a regional focus, and the ultimate intent is to judiciously expand this scope to any observation dataset that has applicability for evaluation of the types of Earth system models used in CMIP.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Synoptic System for Capturing Ecosystem Control Points Across Terrestrial‐Aquatic Interfaces

Interconnected landscape features such as terrestrial‐aquatic interfaces play an outsized role in biogeochemical cycles as ecosystem control points, but it is notoriously challenging to characterize these. Here, we document a synoptic sensor network design that is (a) flexible to accommodate diverse ecosystem interfaces and gradients, (b) adaptable to monitoring and modeling needs of small and large projects alike, (c) standardized for intercomparability across sites and field experiments, and (d) adequately replicated to capture heterogeneity of each parameter monitored. This real‐time monitoring of surface water, groundwater, soil, and vegetation supports configuration and evaluation of models that span upland, wetland, open water strata, and transitions between them. We established the network at seven sites along the Chesapeake Bay and Lake Erie coastlines, including large‐scale flood manipulation experiments in both regions. A central design element is “one data logger program to rule them all”—a collection of sensor‐specific modules deployed on 40 loggers controlling ∼2,000 sensors, with the goal of streamlining maintenance, debugging, and reproducible data processing. The network generates ∼6 M observations per month, capturing system dynamics at the broad spatial and fine temporal scales needed to initialize and benchmark models; measurement frequency can be modified remotely to capture events. This network design has also revealed behaviors not represented in Earth system models, such as transient groundwater oxygen pulses. Completely documented and open source, this standardized, flexible, and efficient sensor network design can reduce barriers to understanding environmental changes and ecosystem responses across systems and scales.

Ward, Nicholas D. [Pacific Northwest National Labo↗
DOE OSTI2020

Assessing mechanical response of CO 2 storage into a depleted carbonate reef using a site-scale geomechanical model calibrated with field tests and InSAR monitoring data

Geomechanical risks of injection have raised concerns regarding secure CO 2 storage. In this work, a combined monitoring and modeling approach is used to assess the stress changes and surface uplift associated with CO 2 injection into a depleted carbonate reef of the Michigan basin. A site-scale geomechanical model is built by assigning mechanical properties of formations using well-log and experimental data. Gravity load is applied to the model to estimate the vertical component of stress as well as different lateral boundary displacement scenarios to estimate horizontal stresses. We used a poroelastic pressure-dependent model (instead of a linear elastic mechanical earth model) to calibrate initial stresses using hydraulic fracture test data measured at depleted reservoir status. Multi-phase fluid flow-geomechanical simulations are performed to estimate the poroelastic response during (1) primary depletion (2) field-scale CO 2 injection phase (3) a hypothetical forecast scenario in which well bottom hole pressure (BHP) reach 45000 KPa. The predicted surface uplift is less than 1 mm at the end of the field-scale CO 2 injection phase which is in good agreement with Interferometric Synthetic Aperture Radar (InSAR) uplift measurement. Although the InSAR data shows an insignificant uplift, hydromechanical modeling of injection shows that CO 2 injection still causes reservoir deformation emphasizing the role of carbonate overburden and reservoir formation mechanical properties and limited size of reef on diminishing the surface deformation. Modeling indicates poroelastic response of caprock matters to estimate uplift. The lower permeability of the top two layers provides additional barrier to large uplift. Also, history of subsidence due to production should be accounted to predict uplift due to a follow up injection correctly. This report shows the significance of combining a calibrated geomechanical model with field measured stresses and monitoring data to be used as a tool to ensure the safety of CO 2 storage.

42 ENGINEERING↗