Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mean”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A meta-analysis of global cropland soil carbon changes due to cover cropping

Including cover crops within agricultural rotations may increase soil organic carbon (SOC). However, contradictory findings generated by on-site experiments make it necessary to perform a comprehensive assessment of interactions between cover crops, environmental and management factors, and changes in SOC. In this study, we collected data from studies that compared agricultural production with and without cover crops, and then analyzed those data using meta-analysis and regression. Our results showed that including cover crops into rotations significantly increased SOC, with an overall mean change of 15.5% (95% confidence interval of 13.8% to 17.3%). Whereas medium-textured soils had highest SOC stocks with (mean = 39 Mg ha-1) and without (mean = 37 Mg ha-1) cover crops, fine-textured soils showed the greatest increase in SOC after the inclusion of cover crops (mean change of 39.5%). Coarse-textured (11.4%) and medium-textured soils (10.3%) had comparatively smaller changes in SOC. Cover crop mixtures resulted in greater increases in SOC compared to mono-species cover crops, and using legumes resulted in greater SOC increases than grasses and broadleaf species. Cover cropping was associated with significant SOC increases in shallow soils (= 30 cm), but not in subsurface soils (> 30 cm). The regression analysis revealed that SOC changes from cover cropping were correlated with improvements in soil quality, specifically decreased runoff and erosion and increased mineralizable carbon, mineralizable nitrogen, and soil nitrogen. Soil carbon change was also affected by annual temperature, number of years after start of cover crop usage, latitude, and initial SOC concentrations. Finally, the mean rate of carbon sequestration from cover cropping across all studies was 0.56 Mg ha-1 yr-1. If 15% of current global cropland were to adopt cover crops, this value would translate to 0.17 Pg of carbon sequestered per year or ~1-2% of current fossil fuels emissions. Altogether, these results indicated that inclusion of cover crops into agricultural rotations can enhance soil carbon concentrations, improve many soil quality parameters, and serve as a potential sink for atmosphere CO2.

Cover crop, soil organic carbon, soil quality, soi↗
DOE OSTI2020

Quantifying Progress Across Different CMIP Phases With the ESMValTool

Abstract More than 40 model groups worldwide are participating in the Coupled Model Intercomparison Project Phase 6 (CMIP6), providing a new and rich source of information to better understand past, present, and future climate change. Here, we use the Earth System Model Evaluation Tool (ESMValTool) to assess the performance of the CMIP6 ensemble compared to the previous generations CMIP3 and CMIP5. While CMIP5 models did not capture the observed pause in the increase in global mean surface temperature between 1998 and 2013, the historical CMIP6 simulations agree well with the observed recent temperature increase, but some models have difficulties in reproducing the observed global mean surface temperature record of the second half of the twentieth century. While systematic biases in annual mean surface temperature and precipitation remain in the CMIP6 multimodel mean, individual models and high‐resolution versions of the models show significant reductions in many long‐standing biases. Some improvements are also found in the vertical temperature, water vapor, and zonal wind speed distributions, and root‐mean‐square errors for selected fields are generally smaller with reduced intermodel spread and higher average skill in the correlation patterns relative to observations. An emerging property of the CMIP6 ensemble is a higher effective climate sensitivity with an increased range between 2.3 and 5.6 K. A possible reason for this increase in some models is improvements in cloud representation resulting in stronger shortwave cloud feedbacks than in their predecessor versions.

Bock, L.↗
DOE OSTI2024

On the subject of large-scale pool fires and turbulent boundary layer interactions

The role to which a realistic inflow turbulent boundary layer (TBL) influences transient and mean large-scale pool fire quantities of interest (QoIs) is numerically investigated. High-fidelity, low-Mach large-eddy simulations that activate low-dissipation, unstructured numerics are conducted using an unsteady flamelet combustion modeling approach with mutiphysics coupling to soot and participating media radiation transport. Three inlet profile configurations are exercised for a large-scale, high-aspect rectangular pool that is oriented perpendicular to the flow direction: a time-varying, TBL inflow profile obtained from a periodic precursor simulation, the time-mean of the transient TBL, and a steady power-law inflow profile that replicates the mean TBL crosswind velocity of 10.0 m/s at a vertical height of 10 m. Results include both qualitative transient flame evolution and quantitative flame shape with ground-level temperature and convective/radiative heat flux profiles. While transient fire events, which are driven by burst-sweep TBL coupling, such as blow-off and reattachment are vastly different in the TBL case (contributing to increased root mean square QoI fluctuation prediction and disparate flame lengths), mean surface QoI magnitudes are similar. Quadrant analysis demonstrates that the TBL configuration modifies burst-sweep phenomena at windward pool locations, while leeward recovery is found. In conclusion, positive fluctuations of convective heat flux correlate with fast moving fluid away from the pool surface due to intermittent combustion events.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The role of turbulent fluctuations in aerosol activation and cloud formation

Aerosol indirect effects are one of the leading contributors to cloud radiative properties relevant to climate. Aerosol particles become cloud droplets when the ambient relative humidity (saturation ratio) exceeds a critical value, which depends on the particle size and chemical composition. In the traditional formulation of this problem, only average, uniform saturation ratios are considered. Using experiments and theory, we examine the effects of fluctuations, produced by turbulence. Our measurements, from a multiphase, turbulent cloud chamber, show a clear transition from a regime in which the mean saturation ratio dominates to one in which the fluctuations determine cloud properties. The laboratory measurements demonstrate cloud formation in mean-subsaturated conditions (i.e., relative humidity <100%) in the fluctuation-dominant activation regime. The theoretical framework developed to interpret these measurements predicts a transition from a mean- to a fluctuation-dominated regime, based on the relative values of the mean and standard deviation of the environmental saturation ratio and the critical saturation ratio at which aerosol particles activate or become droplets. The theory is similar to the concept of stochastic condensation and can be used in the context of the atmosphere to explore the conditions under which droplet activation is driven by fluctuations as opposed to mean supersaturation. Overall, it provides a basis for future development of cloud droplet activation parameterizations that go beyond the internally homogeneous parcel calculations that have been used in the past.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Evaluation of historical precipitation interannual variability in CMIP6 over the United States

Interannual precipitation variability profoundly influences society via its effects on agriculture, water resources, infrastructure, and disaster risks. In this study, we use daily in situ precipitation observations from the global historical climatology network-daily (GHCN-D) to assess the ability of 21 Coupled Model Intercomparison Project Phase 6 (CMIP6) models, including the 50-member fifth-generation Canadian Earth System Model single model initial-condition large ensemble (CanESM5_SMILE), to realistically simulate historical interannual precipitation variability trends within 17 regions of the contiguous United States (CONUS). We assess how accurately the CMIP6 simulations align with observational data across annual, summer, and winter periods, focusing on four key hydrometeorological metrics, including interannual precipitation variability, relative interannual precipitation variability (coefficient of variation), annual mean precipitation, and annual wet day frequency. Our findings reveal that CMIP6 ensemble members generally reproduce the spatial patterns of observed trends in annual mean precipitation. In most regions, models agree well with the signs of observed changes in annual mean precipitation, though discrepancies in trend magnitude are evident. Further, observed trends in winter mean precipitation broadly exhibit a spatial pattern similar to that of the observed annual mean. However, analysis of the CanESM5_SMILE shows that trends in precipitation variability may primarily be the result of model-simulated internal variability, suggesting caution in interpreting multi-model single-realization ensemble results. Challenges in accurately simulating interannual precipitation variability underscore the need for ongoing model refinement and validation to enhance climate projections, especially in regions vulnerable to extreme precipitation events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Diabetes and hypertension are related to amyloid-beta burden in the population-based Rotterdam Study

Abstract Higher vascular disease burden increases the likelihood of developing dementia, including Alzheimer’s disease. Better understanding the association between vascular risk factors and Alzheimer’s disease pathology at the predementia stage is critical for developing effective strategies to delay cognitive decline. In this work, we estimated the impact of six vascular risk factors on the presence and severity of in vivo measured brain amyloid-beta (Aβ) plaques in participants from the population-based Rotterdam Study. Vascular risk factors (hypertension, hypercholesterolaemia, diabetes, obesity, physical inactivity and smoking) were assessed 13 (2004–2008) and 7 years (2009–2014) prior to 18F-florbetaben PET (2018–2021) in 635 dementia-free participants. Vascular risk factors were associated with binary amyloid PET status or continuous PET readouts (standard uptake value ratios, SUVrs) using logistic and linear regression models, respectively, adjusted for age, sex, education, APOE4 risk allele count and time between vascular risk and PET assessment. Participants’ mean age at time of amyloid PET was 69 years (range: 60–90), 325 (51.2%) were women and 190 (29.9%) carried at least one APOE4 risk allele. The adjusted prevalence estimates of an amyloid-positive PET status markedly increased with age [12.8% (95% CI 11.6; 14) in 60–69 years versus 35% (36; 40.8) in 80–89 years age groups] and APOE4 allele count [9.7% (8.8; 10.6) in non-carriers versus 38.4% (36; 40.8) to 60.4% (54; 66.8) in carriers of one or two risk allele(s)]. Diabetes 7 years prior to PET assessment was associated with a higher risk of a positive amyloid status [odds ratio (95% CI) = 3.68 (1.76; 7.61), P < 0.001] and higher standard uptake value ratios, indicating more severe Aβ pathology [standardized beta = 0.40 (0.17; 0.64), P = 0.001]. Hypertension was associated with higher SUVr values in APOE4 carriers (mean SUVr difference of 0.09), but not in non-carriers (mean SUVr difference 0.02; P = 0.005). In contrast, hypercholesterolaemia was related to lower SUVr values in APOE4 carriers (mean SUVr difference −0.06), but not in non-carriers (mean SUVr difference 0.02). Obesity, physical inactivity and smoking were not related to amyloid PET measures. The current findings suggest a contribution of diabetes, hypertension and hypercholesterolaemia to the pathophysiology of Alzheimer’s disease in a general population of older non-demented adults. As these conditions respond well to lifestyle modification and drug treatment, further research should focus on the preventative effect of early risk management on the development of Alzheimer’s disease neuropathology.

van Arendonk, Joyce (ORCID:0000000172941633)↗
DOE OSTI2021

Mass enhancement in 3 d and s - p perovskites from symmetry breaking

In some d-electron oxides, the measured effective mass m* exptl has long been known to be significantly larger than the model effective mass m* model deduced from mean-field band theory, i.e., m* exptl =βm* model , where β>1 is the “mass-enhancement” or “mass-renormalization” factor. Previous applications of density functional theory (DFT), based on a symmetry-restricted structure with the smallest number of possible magnetic, orbital, and structural degrees of freedom, missed such mass enhancement. This fact has been taken as evidence of strong electronic correlation, often described via the symmetry-restricted dynamic mean-field approach of the many-body theory, being the exclusive enabling physics. This paper uses instead a static density functional approach that does not restrict positional or spin degrees of freedom (symmetry-broken structures). This approach analyzes the contributions of different symmetry-broken modalities to mass enhancement for a few nominally highly correlated d-electron perovskites as well as the nominally uncorrelated, closed-shell s-p bonding perovskites. It shows that the energy-lowering symmetry-broken spin effects (e.g., nonzero local moment in the paramagnetic phase) and structural effects (e.g., atomic displacement) as described in mean-field DFT already manifest mass enhancement for both electrons and holes in a range of d-electron perovskites SrVO 3 , SrTiO 3 , BaTiO 3 , and LaMnO 3 , as well as p-electron perovskites CsPbI 3 and SrBiO 3 , including both metals (SrVO 3 ) and insulators (the rest). This is revealed only when enlarged unit cells of the same parent global symmetry, which are large enough to allow for symmetry-breaking distortions and concomitant variations in spin order, are explored for their ability to lower the total energy. Positional symmetry breaking that leads to mass enhancement includes octahedral rotation in halide perovskites such as CsPbI 3 , Jahn-Teller-like Q${+}\atop{2}$ distortion in LaMnO 3 , and bond disproportionation in SrBiO 3 , while magnetic symmetry breaking resulting in mass enhancement includes the formation of a distribution of local moments in SrVO3 that averages to zero in the paramagnetic phase. Not all symmetry breaking leads to significant mass enhancement, e.g., the rather small octahedral rotations in the nearly perfectly cubic SrTiO 3 cause negligible mass enhancement, as do the paraelectric displacements in BaTiO 3 . In principle, by ergodicity, the two descriptions, i.e., the symmetry-restricted dynamic approach with a single time-fluctuating magnetic moment and the symmetry-broken mean-field approach with a static distribution of spatially fluctuated local moments, are not mutually exclusive but are a choice of representation and consequently, a choice of computational efficiency. In approximate implementations, the symmetry-broken mean-field approach appears to remove much of what was strong correlation in dynamically correlated symmetry-restricted solutions, leaving smaller (“weak”) residual correlation with respect to the exact solution.

36 MATERIALS SCIENCE↗
DOE OSTI2021

State densities of heavy nuclei in the static-path plus random-phase approximation

We report that nuclear state densities are important inputs to statistical models of compound-nucleus reactions. State densities are often calculated with self-consistent mean-field approximations that do not include important correlations and must be augmented with empirical collective enhancement factors. Here, we benchmark the static-path plus random-phase approximation (SPA + RPA) to the state density in a chain of samarium isotopes 148–155 Sm against exact results (up to statistical errors) obtained with the shell-model Monte Carlo (SMMC) method. The SPA + RPA method incorporates all static fluctuations beyond the mean field together with small-amplitude quantal fluctuations around each static fluctuation. Using a pairing plus quadrupole interaction, we show that the SPA + RPA state densities agree well with the exact SMMC densities for both the even- and odd-mass isotopes. For the even-mass isotopes, we also compare our results with mean-field state densities calculated with the finite-temperature Hartree-Fock-Bogoliubov (HFB) approximation. We find that the SPA + RPA repairs the deficiencies of the mean-field approximation associated with broken rotational symmetry in deformed nuclei and with the violation of particle-number conservation in the pairing condensate. In particular, in deformed nuclei the SPA + RPA reproduces the rotational enhancement of the state density relative to the mean-field state density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Photometric redshifts and intrinsic alignments: Degeneracies and biases in the 3 × 2 pt analysis

We present a systematic study of cosmological parameter bias in weak lensing and large-scale structure analyses for upcoming imaging surveys induced by the interplay of intrinsic alignments (IA) and photometric redshift (photo-z) model misspecification error. We first examine the degeneracies between the parameters of the tidal alignment–tidal torquing (TATT) model for IA and of a photo-z model including a mean shift ( Δ z ¯ ) and variance ( σ z ) for each tomographic bin of lenses and sources, under a variety of underlying true IA behaviors. We identify strong degeneracies between: (1) the redshift scaling of the tidal alignment amplitude and the mean shift and variances of source bins, (2) the redshift scaling of the tidal torquing amplitude and the variance of the lowest- z source bin, and (3) the IA source density weighting and the mean shift and variance of several source bins. We then use this information to guide our exploration of the level of cosmological parameter bias which can be induced given incorrect modeling of IA, photo-z, or both. We find that marginalizing over all the parameters of TATT is generally sufficient to preclude cosmological parameter bias in the scenarios we consider. However, this does not necessarily mean that IA and photo-z parameters are themselves unbiased, nor does it mean that the best-fit model is a good fit to the data. We also find scenarios where the inferred parameters produce χ d . o . f . 2 values indicative of a good fit but cosmological parameter bias is significant, particularly when the IA source density weighting parameter is not marginalized over. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Once-in-a-lifetime encounter models for neutrino media. II. Quasisteady states and miscidynamic flavor evolution

We extended the once-in-a-lifetime encounter (OILE) model to stochastic interactions among neutrinos. As in the original OILE model, the new model reproduces the mean-field behavior of a dense neutrino gas for time 𝑡 ≲ (𝜇⁢𝛾) −1 , where 𝜇 measures the strength of the mean-field neutrino self-interaction potential and is proportional to the neutrino density, and the dimensionless “impact parameter” 𝛾 is a measure of the change in the flavor quantum state of a neutrino during interaction with another neutrino when the wave packets of the two neutrinos overlap. As in the mean-field case, the OILE model with random neutrino velocities experiences kinetic flavor decoherence as the flavor quantum states of the neutrinos diverge from each other. Unlike the mean-field case, however, the OILE model has a “collision term” due to the quantum entanglement among neutrinos. For 𝛾 ≪ 1, this incoherent effect can drive the neutrinos into a quasisteady state that is similar to the collective precession mode in a homogeneous and isotropic neutrino gas in the mean-field approximation. Subsequently, the collision term drives the neutrino gas adiabatically through different quasisteady states and eventually to flavor equilibration. This process is an example of miscidynamic flavor evolution, with the mixing equilibria being the quasisteady precession states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Quantifying the effects of varietal types × management on the spatial variability of sorghum biomass across US environments

Abstract Regional‐scale estimations of sorghum biomass production allow identification of optimum genotype×environment×management (G×E×M) combinations for bioenergy generation. The objective of this study was to determine the degree of contributions of G, E, and M toward variability in sorghum biomass in the United States. Using the Agricultural Production Systems sIMulator in a grid computing platform, biomass was simulated for irrigated and rainfed conditions for 30 years across the United States for four sorghum varietal types (grain—GS, sudangrass—SS, photosensitive—PS, and photo‐insensitive—PI). Simulated biomass was assessed by environments clustered using the sum of intercepted solar radiation ( ir ), mean of temperature stress factor ( tp ) and water stress factor ( sw ). Simulated biomass ranged from 5.8 t ha −1 (GS‐rainfed) to 27.5 t ha −1 (PI‐irrigated). Under high‐temperature environments (mean annual temperature = 25°C), rainfed biomass between 40 and 80 days after planting (DAP) was strongly correlated with sw ( r = 0.64–0.86) and irrigated biomass with ir ( r = 0.68–0.81). Under low‐temperature environments (mean annual temperature = 18°C) after 40 DAP, tp and ir had greater effects than sw ( r = 0.55–0.82). Biomass variance was mainly explained by varietal type (50%–76%) in all environments×irrigation combinations, except in the high‐ and mid‐temperature environments under rainfed conditions where rainfall had the major effect (25%–45%). However, when mean temperature during the growing season decreased from 25°C (high environments) to 18°C (low environments), the contribution of mean temperature to biomass variance increased from 7% to 34% (rainfed) and from 4% to 36% (irrigated). Varietal type had the larger interactions with other factors independently of the environment and irrigation. We demonstrated a need to quantify (i) the main G×E×M drivers of biomass variability based on environmental stress factors and (ii) the variance contribution of these drivers on sorghum biomass. Our regional‐scale estimations are key inputs for future robust biomass projections of energy sorghum genotypes integrating G×E×M under climate change scenarios.

Ojeda, Jonathan J.↗
DOE OSTI2023

Temperature uncertainty modelling with proxy structural data as geostatistical constraints for well siting: an example applied to Granite Springs Valley, NV, USA

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

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Toward Understanding the Differences between Mesoscale and Large-Eddy Simulations of Tropical Cyclones

In this work, we investigate the ability of mesoscale and large-eddy simulation (LES) model configurations to predict the mean wind speed profile within the boundary layer of tropical cyclones (TCs). To this end, we perform idealized simulations of five hypothetical intense storms ranging from categories 1 to 5 on the Saffir–Simpson scale and extract time-averaged quantities near the eyewall region. We compare the model-generated data against mean wind speed profiles compiled from dropsondes launched from reconnaissance aircraft operating in the North Atlantic basin. Our analysis shows that mesoscale- and LES-generated mean wind fields display important differences in the boundary layer, including the magnitude of shear as well as the height where their low-level wind speed maxima are located. In addition, a comparison between the two model configurations with the dropsonde data shows that both modeling approaches are unable to capture the typical structure of mean winds in the lower part of the TC boundary layer (10–500 m), calling into question the use of simulations of near-axisymmetric storms for investigating the wind structure of past events. To better understand these differences, we conduct a momentum-budget analysis and show that modeled turbulent fluxes are underestimated in the mesoscale boundary layer parameterization compared to the LES model. Based on the analysis of the horizontal turbulent fluxes and their potential impact on mean flow quantities, a TC-specific boundary layer parameterization may be needed.

17 WIND ENERGY↗
DOE OSTI2025

Tropical High Cloud Feedback Relationships to Climate Sensitivity

Clouds constitute a large portion of uncertainty in predictions of equilibrium climate sensitivity (ECS). While low cloud feedbacks have been the focus of intermodel studies due to their high variability among global climate models, tropical high cloud feedbacks also exhibit considerable uncertainty. Here, we apply the cloud radiative kernel technique of Zelinka et al. to 22 models across the CMIP5 and CMIP6 ensembles to survey tropical high cloud feedbacks and analyze their relationship to ECS. We find that the net high cloud feedback and its altitude and optical depth feedback components are significantly positively correlated with ECS in the tropical mean. On the other hand, the tropical mean high cloud amount feedback is not correlated with ECS. These relationships are most pronounced outside of areas of strong climatological ascent, suggesting the importance of thin cirrus feedbacks. Finally, we explore connections between high cloud feedbacks, climate sensitivity, and mean state high cloud properties. In general, high ECS models are cloudier in the upper troposphere but have a thinner high cloud population. Furthermore, we find that having more thin cirrus in the mean state relates to more positive high cloud altitude and optical depth feedbacks, and it either amplifies or dampens the high cloud amount feedback depending on the large-scale dynamical regime (amplifying in descent and dampening in ascent). In summary, our analysis highlights the importance of tropical high cloud feedbacks for driving intermodel spread in ECS and suggests that mean state high cloud characteristics might provide a unique opportunity for observationally constraining high cloud feedbacks.

atmosphere↗
DOE OSTI2024

Surface-to-Interior Transport Timescales and Ventilation Patterns in a Time-Dependent Circulation Driven by Sustained Climate Warming

Abstract The effect of climate warming in response to rising atmospheric CO 2 on the ventilation of the ocean remains uncertain. Here we make theoretical advances in elucidating the relationship between ideal age and transit time distribution (TTD) in a time-dependent flow. Subsequently, we develop an offline tracer-transport model to characterize the ventilation patterns and time scales in the time-evolving circulation for the 1850–2300 period as simulated with the Community Earth System Model version 1 (CESMv1) under a business-as-usual warming scenario. We found that by 2300 2.1% less water originates from the high-latitude deep water formation regions (both hemispheres) compared to 1850. In compensation, there is an increase in the water originating from the subantarctic. We also found that slowing meridional overturning circulation causes a gradual increase in mean age during the 1850–2300 period, with a globally averaged mean-age increase of ∼110 years in 2300. Where and when the water will be re-exposed to the atmosphere depends on the post-2300 circulation. For example, if we assume that the circulation persists in its year-2300 state (scenario 1), the mean interior-to-surface transit time in year 1850 is ∼1140 years. In contrast, if we assume that the circulation abruptly recovers to its year-1850 state (scenario 2), the mean interior-to-surface transit time in 1850 is only ∼740 years. By 2300, these differences become even larger; in scenario 1, the mean interior-to-surface transit time increases by ∼200 years, whereas scenario 2 decreases by ∼80 years. The dependence of interior-to-surface transit time on the future ocean circulation produces an additional unavoidable uncertainty in the long-term durability of marine carbon dioxide removal strategies. Significance Statement The ocean’s circulation, when altered by climate warming, can affect its capacity to absorb heat and CO 2 , which are crucial for the global climate. In our study, we investigated how global warming, caused by rising CO 2 levels, might impact the ocean circulation—the way water moves from deep ocean to the surface and vice versa. We discovered that by 2300, if we continue on our current warming trajectory, the origins of water within the ocean will shift, with less coming from deep, cold zones near the poles and more from subantarctic regions. As a result, deep water will take longer time before it resurfaces than shallow water. How quickly this water travels from deep regions to the surface could change, depending on the state of future ocean circulation. If the circulation remains as predicted in 2300, this journey will take longer. Conversely, if it reverts to the pattern in 1850, the process will be quicker. This variability introduces added uncertainty to strategies aimed at mitigating climate change by storing CO 2 in the ocean. Our work highlights the intricate ways in which climate change can influence our oceans, potentially affecting our plans to mitigate global warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Soil Temperature and Moisture within the Kougarok Fire Complex, Kougarok Road Mile Marker 86, Seward Peninsula, Alaska, 2019-2023

Daily averages of soil temperature and moisture measured once every hour at different heights located at Intensive Monitoring Stations within the Kougarok Fire Complex, Kougarok Road Mile Marker 86 site. Data were retrieved annually from 2019-2023. Package contains 21 *.CSV data files plus a file level metadata *.CSV, data dictionary *.CSV, data file inventory *.CSV, and sensor location site map *.JPG. Data files have header rows, NaN fields indicate invalid or missing data, and negative vertical offsets are above ground.The Kougarok tundra fire complex (KFC) is located north of Nome and the Kigluaik Mountains, near Quartz Creek and the Kougarok River. The site is accessed by foot from the end of the Nome-Taylor Highway (mile marker 86; also called the Kougarok or Beam Road). The KFC burned in six major fires in the decades since 1950 (Alaska Interagency Coordination Center, unpublished data). Lightning ignited five of these fires (1971, 1997, 2015, and 2019) and one was human caused (2002). The mosaic of overlapping fire scars allows for the study of repeat fires in the tundra which, until recently, was not a common phenomenon outside the boreal forest in Alaska. Our reference unburned tundra fire site is south of the KFC located at mile marker 80 of the Nome-Taylor Highway.The two most recent fires are the Mingvk Lake (2015; 21,698 acres burned from 7/27/2015 to 9/28/2015) and Garfield Creek (2019; 422 acres burned from 7/31/2019 to 8/20/19). The Mingvk Lake fire scar includes areas that burned 1-4x (1971, 1997, 2002), while the entirety of the Garfield Creek fire scar has burned 2x previously (1971, 2002).Previous research at the KFC focused on permafrost (Liljedahl et al. 2007; Narita et al. 2015; Iwahana et al. 2016; Tsuyuzaki, Iwahana, and Saito 2017) and vegetation (Narita et al. 2015; Hollingsworth et al. 2021) response to fire. The central Seward Peninsula is characterized by continuous permafrost with a thickness of 15 to 30 m and a mean active layer thickness of 56 cm (Hinzman et al. 2003). Sloping hills with mixed shrub–tussock tundra and tussock tundra vegetation in the uplands are characteristic of the region. Three micrometeorological towers near the Kougarok field site recorded a mean annual temperature of −2.4°C, mean January temperature of −23.1°C, mean July temperature of +11°C, and mean summer rainfall (June–August) of 94 mm from 2000 to 2006 (Liljedahl et al. 2007).The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

3D CFD Model Validation Using Benchmark Data of 1/16th Scaled VHTR Upper Plenum and Development of Wall Heat-Transfer Correlation For Laminar Flow

With support from the U.S. Department of Energy-Office of Nuclear Energy’s (DOE-NE’s) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, an effort has been pursued to support high-temperature gas-cooled reactor (HTGR) technology development and its modeling and simulation needs. There is a particular need for advanced modeling and simulation tools to predict thermal-fluid behavior in the nuclear reactor primary system, especially in the core and the lower and upper plena, during safety-related transients. In this report, two main such activities are presented relevant to the HTGRs: (1) three-dimensional (3D) computational fluid dynamics (CFD) validation using benchmark data from the upper plenum of Texas A&M University’s 1/16th scaled very-high-temperature gas-cooled reactor (VHTR), and (2) development of wall heat-transfer correlation for laminar flow in a wall-heated pipe. The CFD tool validation exercises can be helpful to choose the models and CFD tools to simulate and design specific components of the HTRGs such as upper plenum where jet mixing is a complex phenomenon. In a loss of forced circulation event, the laminar flow can be observed during the development of natural circulation flow. This work includes the development and validation of heat transfer correlations for laminar flow using the Nek5000 CFD code due to limited available experimental data for laminar flow conditions to guide low-order models (1D). In this report, the flow characteristics of a single isothermal jet discharging into the upper plenum was investigated using the Nek5000 Large-Eddy Simulation (LES) CFD tool. Several numerical simulations were performed for various jet-discharged Reynolds numbers ranging from 3,413 to 12,819. A grid-independent study was performed. The numerical results of mean velocity, root-mean-square fluctuating velocity, and Reynolds stress were compared against the benchmark data. Good agreement was obtained between simulated and measured data for axial mean velocities, except near the upper plenum hemisphere. The maximum predicted errors for axial mean velocities at various normalized coolant channel diameter heights of 1, 5, and 10 are 1.56%, 1.88%, and 3.82%, respectively. In addition, the predicted root-mean-square fluctuating velocity and Reynolds stress are qualitatively in agreement with the experimental data. The Nek5000 code was used to develop wall-heat transfer correlation for laminar flow in a cylindrical tube. Several simulations were performed for various Reynolds flow and wall-heat fluxes. A new heat transfer correlation was developed using data from Nek5000 simulation results and regression functions in Matlab. The developed heat transfer correlation is valid for various Reynolds flows from 200 to 2000. The predicted R² value for model fit was 0.875, which ensures that 87.5% of the model data lies on the Nek5000 data. Moreover, a machine learning (ML) tool was used to train and test the Nek5000 data. A good fit of the ML-based model was observed with the test data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Rare Earth Elements (REE) and Critical Minerals (CM) in Middle Pennsylvanian-Age Coals and Associated Sediments in the Central Appalachian Basin, Eastern U.S.A.

The Central Appalachian Basin (CAB) of Kentucky, Tennessee, West Virginia and Virginia has a long history of coal mining and oil and gas extraction that has empowered the regional and national economies, the development of infrastructure, and a highly trained energy resources work force. As our societal demands for advanced technologies have rapidly increased in recent years, coal-related materials are viewed as an important new unconventional domestic source of critical minerals (CM) that are required for telecommunications, aerospace and transportation industries, electronics, the transition to low-carbon emissions energy production, and many consumer products. Coal-related materials encompass coal, associated sediments, coal mining waste materials, produced waters, and ash residues from coal-fired power plants. An important objective of the Evolve Central Appalachia Project (Evolve CAPP), sponsored by the U.S. Department of Energy (DOE) National Energy Technology Laboratory (NETL) is to assess the quantity and distribution of CM resources in the CAB region. The rare earth elements (REE) are considered highest priority, although other important CM such as niobium, gallium, and zirconium are known to occur in the Middle Pennsylvanian-age coals and sediments. Working with coal industry partners who provided access to drill cores, coal-related sediments and waste materials, over 600 samples have been collected for laboratory analysis, and over 730 materials have been scanned using portable x-ray fluorescence (pXRF) equipment. The application of pXRF provides the means for real-time semiquantitative analysis of CM content at very close spacings (typically 2-3 inch intervals) along drill core and in-situ channel samples that span the roof, coal seam, and floor rock. The comparison of pXRF geochemical data with laboratory results, geologic data, and downhole spectral gamma logs can provide high resolution input to lithologic and depositional models for future CM resource evaluations. The preliminary findings show that pXRF is capable of accurately measuring low concentrations of many of the CM with a high level of confidence (Ba, Cr, Ga, K, Nb, Rb, Sr, Th, Y), whereas for others (La, Ce, Co, Mn, Nd, Ni, Sc, Ti, V) the detection limits are very high or spectral interferences increase the uncertainty. Notably, the mean abundances of Y (34 ppm), La (86 ppm), Ga (40 ppm), and V (146 ppm) in coal underclays in the CAB region are up to 7X enriched compared with the overlying coal. These values also exceed the reported concentrations in published reference materials for upper continental crustal rocks (Rudnick and Gao, 2003), North American Shale Composite (Gromet et al., 1984), and North American coal (Finkelman, 1993). The pXRF data are in part verified by laboratory results that indicate the mean Y abundance (36 ppm) is highly correlated (R2 = 0.807) with total REE (ΣREE). The correlation is even higher (R2 = 0.957) with heavy REE (ΣHREE). Applying these correlations to the pXRF data for the coal underclays, the mean estimated values for ΣREE+Y and ΣHREE+Y are 270 ppm and 58 ppm, respectively. Although these average values are not considered high, the range of Y measured by pXRF in the coal underclays extended as high as 114 ppm, which would suggest ΣREE+Y equal to 847 ppm. The mean abundances of Zr (192 ppm) and Th (21 ppm) in coal underclays are also enriched compared with the overlying coal and these results likely reflect the presence of resistant detrital heavy minerals such as monazite, xenotime, and zircon in the underclay matrix. Several of the profiled coal seams and associated wall rocks contained thin volcanic ash layers up to 4-5 inches in thickness. The extent to which these ash fall layers provided a source for CM under the paleoenvironmental conditions that resulted in coal deposits in the CAB remains to be fully studied. Continuing investigations in the Evolve CAPP study area will include laboratory determinations of mineralogic and clay compositions, and evaluations of CM geochemical mobility in the coal and coal underclays.

Lassetter, Billy↗