SEARCH · Engineering Papers

Results for “Scale down”

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.

SourceResearch connectionFrom yearThrough year

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

At least 829 records · Page 46

All sources NASA NTRS DOE OSTI

DOE OSTI2025

On the emergent scale of bedrock groundwater contribution to headwater mountain streams

We investigated the contribution of bedrock groundwater to streamflow as a function of catchment scale in a headwater stream. Synoptic surveys were conducted during hydrologically important periods of the year using multiple environmental tracers in stream water, soil water, and bedrock groundwater, along a first-order montane stream, in west-central Montana. Sampled analytes included 222 Rn, used to constrain total subsurface flux, and major and minor elements, used in end-member mixing analysis (EMMA) to identify the contributions of soil and bedrock groundwater to the stream. Partitioning between soil-derived and bedrock-derived groundwater was then analyzed as a function of the incremental and accumulated sub-catchment sizes. Radon results indicated that subsurface water contributions accounted for the majority of streamflow at all surveyed times. EMMA results revealed that the bedrock groundwater contribution to streamflow varied between 26% during peak snowmelt and 44% during late summer. Streamflow generation was dominated by soil groundwater contribution along the entire reach, but the bedrock groundwater contribution increased consistently with accumulated sub-catchment size. However, groundwater contributions were not well-correlated with incremental sub-catchment size. The scale at which increased bedrock groundwater discharge can be correlated with sub-catchment size appears to be >1 km 2 for our study. Our results are consistent with a conceptual model where streamflow is predominantly generated by a 3D subsurface nested flow system. Local subsurface heterogeneities control the stream source at local scales but begin to average out at scales >2 km 2 . Our study indicates that, while soil groundwater is the dominant source, bedrock groundwater remains an important and predictable contributor to streamflow throughout the year, even in a snow-dominated, mountainous headwater catchment.

environmental tracers↗

DOE OSTI2022

Dynamical System Scaling of a Thermocline Thermal Storage System in the Thermal Energy Distribution System (TEDS) Facility

The purpose of this study was to develop a process to convert input signals from one facility into another by reflecting geometric and environmental settings. The Dynamic Energy Transport and Integration Laboratory (DETAIL) is one facility in development that aims to emulate the daily interactions among power production industry systems and be capable of receiving real-time data from those systems as inputs. To convert signals and ensure that the temporal sequences and magnitudes reflect the laboratory settings, the ability to scale and project data is essential. To demonstrate this ability, Dynamical System Scaling (a methodology that enables systems to scale and project or extrapolate datasets to desired environments while conserving the observed transient behavior based on first principles) was applied to DETAIL’s thermocline thermal storage system in the Thermal Energy Distribution System. The thermocline system was successfully scaled and a test case was conducted to generate a doubly accelerated energy charge and discharge in reference to past experimental data from the facility. The accelerated data were determined as able to conserve the amount of energy stored and the associated test boundary conditions were charge line maximum temperature, charge line velocity, and thermocline maximum temperature at 354 °C, 0.458 m/s, and 418 °C, respectively. The research results represented a case that required signals to be accelerated without altering the stored energy.

Yoshiura, Ramon (ORCID:0000000319781313)↗

DOE OSTI2023

Vegetation Management Cost and Maintenance Implications of Different Ground Covers at Utility-Scale Solar Sites

Utility-scale solar photovoltaics (PV) is the largest and fastest-growing sector of the solar energy market, and plays an important role in ensuring that state and local jurisdictions can meet renewable energy targets. Potential adverse environmental impacts of utility-scale solar PV are well-documented, and the effects of diverse mitigation and dual land use strategies under the banner of ’low-impact solar’ are justly receiving more attention; this article seeks to contribute to improving understanding of this topic. Capital costs for different PV configurations are well-documented; however, operation and maintenance (O&M) costs for vegetation management at low-impact utility-scale solar PV sites are not as well-understood, particularly as they compare to costs for sites that use more conventional ground cover practices, such as turfgrass or gravel. After a literature review of different vegetation strategies and O&M cost considerations, we collected data from utility-scale solar PV O&M stakeholders, including site owners/operators, O&M service providers, vegetation maintenance companies, and solar graziers, on costs and activities associated with vegetation management at low-impact, agrivoltaic, and conventional PV sites. In this paper, we perform data analysis to detail the per-activity and total O&M costs for vegetation management at PV sites with different ground covers and management practices, providing the most comprehensive and detailed assessment of PV vegetation O&M costs to date. For the 54 sites included in our analysis, we found that while the per-acre and per-kilowattdc (kWdc) costs for individual activities, such as mowing, trimming, and herbicide application at native or pollinator friendly ground covers, were lower than at turfgrass sites, the total combined vegetation O&M costs were slightly higher; this is presumably because more individual activities are required for the first 3–5 years of vegetation establishment. Qualitative results include recommendations from data providers for site and system design, and ongoing vegetation management operations.

14 SOLAR ENERGY↗

DOE OSTI2022

Cell-free Scaled Production and Adjuvant Addition to a Recombinant Major Outer Membrane Protein from Chlamydia muridarum for Vaccine Development

Subunit vaccines offer advantages over more traditional inactivated or attenuated whole-cell-derived vaccines in safety, stability, and standard manufacturing. To achieve an effective protein-based subunit vaccine, the protein antigen often needs to adopt a native-like conformation. This is particularly important for pathogensurface antigens that are membrane-bound proteins. Cell-free methods have been successfully used to produce correctly folded functional membrane protein through the co-translation of nanolipoprotein particles (NLPs), commonly known as nanodiscs. This strategy can be used to produce subunit vaccines consisting of membrane proteins in a lipid-bound environment. However, cell-free protein production is often limited to small scale (<1 mL). The amount of protein produced in small-scale production runs is usually sufficient for biochemical and biophysical studies. However, the cell-free process needs to be scaled up, optimized, and carefully tested to obtain enough protein for vaccine studies in animal models. Other processes involved in vaccine production, such as purification, adjuvant addition, and lyophilization, need to be optimized in parallel. This paper reports the development of a scaled-up protocol to express, purify, and formulate a membrane-bound protein subunit vaccine

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Parsec-scale Properties of Steep- and Flat-spectrum Extragalactic Radio Sources from a VLBA Survey of a Complete North Polar Cap Sample

We observed with the Very Long Baseline Array (VLBA) at 2.3 and 8.6 GHz a complete flux-density-limited sample of 482 radio sources with decl. > +75° brighter than 200 mJy at 1.4 GHz drawn from the NVSS catalog. A total of 34% of the sources show parsec-scale emission above the flux density detection limit of 30 mJy; their accurate positions and parsec-scale structure parameters are determined. Among all the sources detected at least at the shortest VLBA baselines, the majority, or 72%, have a steep single-dish spectrum. The fraction of the sources with a detectable parsec-scale structure is above 95% among the flat-spectrum objects and close to 25% among the steep-spectrum objects. We identified 82 compact steep-spectrum source candidates, which make up 17% of the sample; most of them are reported for the first time. The compactness and the brightness temperature of the sources in our sample show a positive correlation with single-dish and VLBA spectral indices. All the sources with a significant 8 GHz variability were detected by the VLBA snapshot observations, which independently confirmed their compactness. We demonstrated that 54% of the sources detected by the VLBA at 2.3 GHz in our sample have a steep VLBA spectrum. The compact radio emission of these sources is likely dominated by optically thin jets or mini-lobes, not by an opaque jet core. These results show that future VLBI surveys aimed at searching for new sources with parsec-scale structure should include not only flat-spectrum sources but also steep-spectrum ones in order to reach an acceptable level of completeness.

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

Four-year Cosmology Large Angular Scale Surveyor (CLASS) Observations: On-sky Receiver Performance at 40, 90, 150, and 220 GHz Frequency Bands

The Cosmology Large Angular Scale Surveyor (CLASS) observes the polarized cosmic microwave background (CMB) over the angular scales of 1° ≲ θ ≤ 90° with the aim of characterizing primordial gravitational waves and cosmic reionization. We report on the on-sky performance of the CLASS Q-band (40 GHz), W-band (90 GHz), and dichroic G-band (150/220 GHz) receivers that have been operational at the CLASS site in the Atacama desert since 2016 June, 2018 May, and 2019 September, respectively. We show that the noise-equivalent power measured by the detectors matches the expected noise model based on on-sky optical loading and lab-measured detector parameters. Using Moon, Venus, and Jupiter observations, we obtain power to antenna temperature calibrations and optical efficiencies for the telescopes. From the CMB survey data, we compute instantaneous array noise-equivalent-temperature sensitivities of 22, 19, 23, and 71 $\mu {{\rm{K}}}_{\mathrm{cmb}}\sqrt{{\rm{s}}}$ for the 40, 90, 150, and 220 GHz frequency bands, respectively. These noise temperatures refer to white noise amplitudes, which contribute to sky maps at all angular scales. Future papers will assess additional noise sources impacting larger angular scales.

Observational cosmology↗

DOE OSTI2022

Dissecting Nearby Galaxies with piXedfit. II. Spatially Resolved Scaling Relations among Stars, Dust, and Gas

We study spatially resolved scaling relations among stars, dust, and gas in ten nearby spiral galaxies. In a preceding paper, we have derived spatially resolved properties of the stellar population and dust by a panchromatic spectral energy distribution fitting using piXedfit. Now, we investigate resolved star formation (${{\rm{\Sigma }}}_{{{\rm{H}}}_{2}}$–ΣSFR–Σ * ) and dust scaling relations. While the relations with all subgalactic regions of the galaxies are reasonably tight (σ ≲ 0.3 dex), we find that most of the scaling relations exhibit galaxy-to-galaxy variations in normalization and shape. Only two relations of Σdust–Σgas and Σdust–${{\rm{\Sigma }}}_{{{\rm{H}}}_{2}}$ do not show noticeable galaxy-to-galaxy variations among our sample galaxies. We further investigate the correlations among the scaling relations. We find significant correlations among the normalization of the ${{\rm{\Sigma }}}_{{{\rm{H}}}_{2}}$–ΣSFR–Σ * relations, which suggest that the galaxies with higher levels of resolved H 2 fraction (${f}_{{{\rm{H}}}_{2}}$) tend to have higher levels of resolved star formation efficiency (SFE) and specific star formation rate (sSFR). We also observe that the galaxies with higher levels of resolved dust-to-stellar mass ratios tend to have higher levels of resolved sSFR, SFE, and ${f}_{{{\rm{H}}}_{2}}$. Moreover, we find that the galaxies with higher global sSFR and less compact morphology tend to have higher levels of the resolved sSFR, SFE, and ${f}_{{{\rm{H}}}_{2}}$, which can explain the variations in the normalization of the ${{\rm{\Sigma }}}_{{{\rm{H}}}_{2}}$–ΣSFR–Σ * relationships. Overall, we observe indications of the contributions of both global and local factors in governing the star formation process in galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Constraints on the Correlation of IceCube Neutrinos with a Tracer of Nearby Large-scale Structure

The IceCube Neutrino Observatory has observed extragalactic astrophysical neutrinos with an apparently isotropic distribution. Only a small fraction of the observed astrophysical neutrinos can be explained by known sources. Neutrino production is thought to occur in energetic environments that are ultimately powered by the gravitational collapse of dense regions of the large-scale mass distribution in the universe. Whatever their identity, neutrino sources likely trace this large-scale mass distribution. The clustering of neutrinos with a tracer of the large-scale structure may provide insight into the distribution of neutrino sources with respect to redshift and the identity of neutrino sources. We implement a two-point angular cross correlation of the Northern sky track events with an infrared galaxy catalog derived from the Wide-field Infrared Survey Explorer (WISE) and Two Micron All Sky Survey (2MASS) source catalogs, which trace the nearby large-scale structure. No statistically significant correlation is found between the neutrinos and this infrared galaxy catalog. We find that ≤54% of the diffuse muon neutrino flux can be attributed to sources correlated with the galaxy catalog with 90% confidence. Additionally, when assuming that the neutrino source comoving density evolves following a power law in redshift, dN s /dV ∝ (1 + z) k , we find that sources with negative evolution, in particular k < −1.75, are disfavored at the 90% confidence level.

Abbasi, R. [Loyola University, Chicago, IL (United↗

DOE OSTI2022

Bridging the Bondi and Event Horizon Scales: 3D GRMHD Simulations Reveal X-shaped Radio Galaxy Morphology

X-shaped radio galaxies (XRGs) produce misaligned X-shaped jet pairs and make up ≲10% of radio galaxies. XRGs are thought to emerge in galaxies featuring a binary supermassive black hole (SMBH), SMBH merger, or large-scale ambient medium asymmetry. We demonstrate that XRG morphology can naturally form without such special, preexisting conditions. Our 3D general-relativistic magnetohydrodynamic (GRMHD) simulation for the first time follows magnetized rotating gas from outside the SMBH sphere of influence of radius RB to the SMBH of gravitational radius Rg at the largest scale separation, RB/Rg = 103, to date. Initially, our axisymmetric system of constant-density hot gas contains a weak vertical magnetic field and rotates in the equatorial plane of a rapidly spinning SMBH. We seed the gas with small-scale 2% level pressure perturbations. Infalling gas forms an accretion disk, and the SMBH launches relativistically magnetized collimated jets reaching well outside RB. Under the pressure of the infalling gas, the jets intermittently turn on and off, erratically wobble, and inflate pairs of cavities in different directions, resembling an X-shaped jet morphology. Synthetic X-ray images reveal multiple pairs of jet-powered shocks and cavities. Large-scale magnetic flux accumulates on the SMBH, becomes dynamically important, and leads to a magnetically arrested disk state. The SMBH accretes at 2% of the Bondi rate ($\dot{M}\simeq 2.4\times {10}^{-3}{M}_{\odot }\,{\mathrm{yr}}^{-1}$ for M87*) and launches twin jets at η = 150% efficiency. These jets are powerful enough (P jets ≃ 2 × 10 44 erg s -1 ) to escape along the SMBH spin axis and end the short-lived intermittent jet state, whose transient nature can account for the rarity of XRGs.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

As a Matter of Dynamical Range – Scale Dependent Energy Dynamics in MHD Turbulence

Abstract Magnetized turbulence is ubiquitous in many astrophysical and terrestrial plasmas but no universal theory exists. Even the detailed energy dynamics in magnetohydrodynamic (MHD) turbulence are still not well understood. We present a suite of subsonic, super-Alfvénic, high plasma beta MHD turbulence simulations that only vary in their dynamical range, i.e., in their separation between the large-scale forcing and dissipation scales, and their dissipation mechanism (implicit large eddy simulation, ILES, and direct numerical simulation (DNS)). Using an energy transfer analysis framework we calculate the effective numerical viscosities and resistivities, and demonstrate that all ILES calculations of MHD turbulence are resolved and correspond to an equivalent visco-resistive MHD turbulence calculation. Increasing the number of grid points used in an ILES corresponds to lowering the dissipation coefficients, i.e., larger (kinetic and magnetic) Reynolds numbers for a constant forcing scale. Independently, we use this same framework to demonstrate that—contrary to hydrodynamic turbulence—the cross-scale energy fluxes are not constant in MHD turbulence. This applies both to different mediators (such as cascade processes or magnetic tension) for a given dynamical range as well as to a dependence on the dynamical range itself, which determines the physical properties of the flow. We do not observe any indication of convergence even at the highest resolution (largest Reynolds numbers) simulation at 2048 3 cells, calling into question whether an asymptotic regime in MHD turbulence exists, and, if so, what it looks like.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Bridging Scales in Black Hole Accretion and Feedback: Magnetized Bondi Accretion in 3D GRMHD

Fueling and feedback couple supermassive black holes (SMBHs) to their host galaxies across many orders of magnitude in spatial and temporal scales, making this problem notoriously challenging to simulate. We use a multi-zone computational method based on the general relativistic magnetohydrodynamic (GRMHD) code KHARMA that allows us to span 7 orders of magnitude in spatial scale, to simulate accretion onto a non-spinning SMBH from an external medium with a Bondi radius of R B ≈ 2 × 10 5 GM•/c 2 , where M• is the SMBH mass. For the classic idealized Bondi problem, spherical gas accretion without magnetic fields, our simulation results agree very well with the general relativistic analytic solution. Meanwhile, when the accreting gas is magnetized, the SMBH magnetosphere becomes saturated with a strong magnetic field. The density profile varies as ~r -1 rather than r -3/2 and the accretion rate $\dot{M}$ is consequently suppressed by over 2 orders of magnitude below the Bondi rate $\dot{M}$ B . Here, we find continuous energy feedback from the accretion flow to the external medium at a level of $\thicksim$ 10 -2 $\dot{M}$c 2 $\thicksim$ 5 X 10 -5 $\dot{M}$ B c 2 . Energy transport across these widely disparate scales occurs via turbulent convection triggered by magnetic field reconnection near the SMBH. Thus, strong magnetic fields that accumulate on horizon scales transform the flow dynamics far from the SMBH and naturally explain observed extremely low accretion rates compared to the Bondi rate, as well as at least part of the energy feedback.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Quantum Scales of Galaxies from Self-interacting Ultralight Dark Matter

We derive the characteristic scales for physical quantities of dwarf galaxies, such as mass, size, acceleration, and angular momentum, within the self-interacting ultralight dark matter (ULDM) model. Due to the small mass of ULDM, even minor self-interactions can drastically alter these scales in the Thomas-Fermi limit. We suggest that these characteristic scales are connected to mysteries of observed galaxies. Oscillation of ULDM field can explain the current cosmological density of dark matter. Many cosmological constraints suggest that the energy scale $\tilde{m}$ for self-interacting ULDM is typically of the order $10~eV$, whereas the mass $m$ for the non-interacting case is around $10^{-21}~eV$. Self-interacting ULDM provides the better explanation for cosmological observations than the non-interacting case.

Astrophysics of Galaxies (astro-ph.GA)↗

DOE OSTI2025

Holography vs. Scale Separation

In this work, we point out a contradiction between holography and scale-separated AdS (i.e. parametrically large mass gap) in string theory, making the standard assumption that the holographic CFT describes the IR degrees of freedom on a brane that decouple from gravity. We show that the CFT can only decouple from gravity if the scalar potential in the dual AdS satisfies a certain criterion. Namely, there must exist a scalar field trajectory that follows the gradient of the scalar potential to the asymptotic region of the scalar field space in which limit $\partial_ϕ\ln(V)\partial_ϕ\ln(Λ_s)\leq2/(d-2)$, where $Λ_s$ is the quantum gravity cut-off. This condition, which generically implies lack of scale-separation, is satisfied in the standard examples of AdS/CFT. However, proposed attempts at achieving scale separation, such as DGKT, employ scalar potentials that violate this condition. We therefore conclude that the CFT duals of DGKT vacua cannot exist in string theory. Barring fine-tuning, our conclusions apply to other Ricci-flat flux compactifications including the KKLT scenario which relies on scale separation to obtain a metastable de Sitter uplift.

FOS: Physical sciences↗

DOE OSTI2022

Using aircraft measurements to characterize subgrid-scale variability of aerosol properties near the Atmospheric Radiation Measurement Southern Great Plains site

Abstract. Complex distributions of aerosol properties evolve in space and time as a function of emissions, new particle formation, coagulation, condensational growth, chemical transformation, phase changes, turbulent mixing and transport, removal processes, and ambient meteorological conditions. The ability of chemical transport models to represent the multi-scale processes affecting the life cycle of aerosols depends on their spatial resolution since aerosol properties are assumed to be constant within a grid cell. Subgrid-scale-dependent processes that affect aerosol populations could have a significant impact on the formation of particles, their growth to cloud condensation nuclei (CCN) sizes, aerosol–cloud interactions, dry deposition and rainout and hence their burdens, lifetimes, and radiative forcing. To address this issue, we characterize subgrid-scale variability in terms of measured aerosol number, size, composition, hygroscopicity, and CCN concentrations made by repeated aircraft flight paths over the Atmospheric Radiation Measurement (ARM) program's Southern Great Plains (SGP) site during the Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystem (HI-SCALE) campaign. Subgrid variability is quantified in terms of both normalized frequency distributions and percentage difference percentiles using grid spacings of 3, 9, 27, and 81 km that represent those typically used by cloud-system-resolving models as well as the current and next-generation climate models. Even though the SGP site is a rural location, surprisingly large horizontal gradients in aerosol properties were frequently observed. For example, 90 % of the 3, 9, and 27 km cell mean organic matter concentrations differed from the 81 km cell around the SGP site by as much as ∼ 46 %, large spatial variability in aerosol number concentrations and size distributions were found during new particle formation events, and consequently 90 % of the 3, 9, and 27 km cell mean CCN number concentrations differed from the 81 km cell mean by as much as ∼ 38 %. The spatial variability varied seasonally for some aerosol properties, with some having larger spatial variability during the spring and others having larger variability during the late summer. While measurements at a single surface site cannot reflect the surrounding variability of aerosol properties at a given time, aircraft measurements that are averaged within an 81 km cell were found to be similar to many, but not all, aerosol properties measured at the ground SGP site. This analysis suggests that it is reasonable to directly compare most ground SGP site aerosol measurements with coarse global climate model predictions. In addition, the variability quantified by the aircraft can be used as an uncertainty range when comparing the surface point measurements with model predictions that use coarse grid spacings.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Cloud adjustments from large-scale smoke–circulation interactions strongly modulate the southeastern Atlantic stratocumulus-to-cumulus transition

Smoke from southern Africa blankets the southeastern Atlantic Ocean from June to October, producing strong and competing aerosol radiative effects. Smoke effects on the transition between overcast stratocumulus and scattered cumulus clouds are investigated along a Lagrangian (air-mass-following) trajectory in regional climate and large eddy simulation models. Results are compared with observations from three recent field campaigns that took place in August 2017: ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES), CLouds and Aerosol Radiative Impacts and Forcing: Year 2017 (CLARIFY), and Layered Atlantic Smoke Interactions with Clouds (LASIC). The case study is set up around the joint ORACLES–CLARIFY flight that took place near Ascension Island on 18 August 2017. Smoke sampled upstream on an ORACLES flight on 15 August 2017 likely entrained into the marine boundary layer later sampled during the joint flight. The case is first simulated with the WRF-CAM5 regional climate model in three distinct setups: (1) FireOn, in which smoke emissions and any resulting smoke–cloud–radiation interactions are included; (2) FireOff, in which no smoke emissions are included; (3) RadOff, in which smoke emissions and their microphysical effects are included but aerosol does not interact directly with radiation. Over the course of the Lagrangian trajectory, differences in free tropospheric thermodynamic properties between FireOn and FireOff are nearly identical to those between FireOn and RadOff, showing that aerosol–radiation interactions are primarily responsible for the free tropospheric effects. These effects are non-intuitive: in addition to the expected heating within the core of the smoke plume, there is also a “banding” effect of cooler temperature (~1–2 K) and greatly enhanced moisture (>2 g kg -1 ) at the plume top. This banding effect is caused by a vertical displacement of the former continental boundary layer in the free troposphere in the FireOn simulation resulting from anomalous diabatic heating due to smoke absorption of sunlight that manifests primarily as a few hundred meters per day reduction in large-scale subsidence over the ocean. A large eddy simulation (LES) is then forced with free tropospheric fields taken from the outputs for the WRF-CAM5 FireOn and FireOff runs. Cases are run by selectively perturbing one variable (e.g., aerosol number concentration, temperature, moisture, vertical velocity) at a time to better understand the contributions from different indirect (microphysical), “large-scale” semi-direct (above-cloud thermodynamic and subsidence changes), and “local” semi-direct (below-cloud smoke absorption) effects. Despite a more than 5-fold increase in cloud droplet number concentration when including smoke aerosol concentrations, minimal differences in cloud fraction evolution are simulated by the LES when comparing the base case with a perturbed aerosol case with identical thermodynamic and dynamic forcings. A factor of 2 decrease in background free tropospheric aerosol concentrations from the FireOff simulation shifts the cloud evolution from a classical entrainment-driven “deepening–warming” transition to trade cumulus to a precipitation-driven “drizzle-depletion” transition to open cells, however. The thermodynamic and dynamic changes caused by the WRF-simulated large-scale adjustments to smoke diabatic heating strongly influence cloud evolution in terms of both the rate of deepening (especially for changes in the inversion temperature jump and in subsidence) and in cloud fraction on the final day of the simulation (especially for the moisture “banding” effect). Such large-scale semi-direct effects would not have been possible to simulate using a small-domain LES model alone.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Global 1 km land surface parameters for kilometer-scale Earth system modeling

Earth system models (ESMs) are progressively advancing towards the kilometer scale (“k-scale”). However, the surface parameters for land surface models (LSMs) within ESMs running at the k-scale are typically derived from coarse-resolution and outdated datasets. This study aims to develop a new set of global land surface parameters with a resolution of 1 km for multiple years from 2001 to 2020, utilizing the latest and most accurate available datasets. Specifically, the datasets consist of parameters related to land use and land cover, vegetation, soil, and topography. Differences between the newly developed 1 km land surface parameters and conventional parameters emphasize their potential for higher accuracy due to the incorporation of the most advanced and latest data sources. To demonstrate the capability of these new parameters, we conducted 1 km resolution simulations using the E3SM Land Model version 2 (ELM2) over the contiguous United States. Our results demonstrate that land surface parameters contribute to significant spatial heterogeneity in ELM2 simulations of soil moisture, latent heat, emitted longwave radiation, and absorbed shortwave radiation. On average, about 31 % to 54 % of spatial information is lost by upscaling the 1 km ELM2 simulations to a 12 km resolution. Using eXplainable Machine Learning (XML) methods, the influential factors driving the spatial variability and spatial information loss of ELM2 simulations were identified, highlighting the substantial impact of the spatial variability and information loss of various land surface parameters, as well as the mean climate conditions. The comparison against four benchmark datasets indicates that ELM generally performs well in simulating soil moisture and surface energy fluxes. The new land surface parameters are tailored to meet the emerging needs of k-scale LSM and ESM modeling with significant implications for advancing our understanding of water, carbon, and energy cycles under global change. The 1 km land surface parameters are publicly available at https://doi.org/10.5281/zenodo.10815170 (Li et al., 2024).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

WRF-Comfort: simulating microscale variability in outdoor heat stress at the city scale with a mesoscale model

Abstract. Urban overheating and its ongoing exacerbation due to global warming and urban development lead to increased exposure to urban heat and increased thermal discomfort and heat stress. To quantify thermal stress, specific indices have been proposed that depend on air temperature, mean radiant temperature (MRT), wind speed, and relative humidity. While temperature and humidity vary on scales of hundreds of meters, MRT and wind speed are strongly affected by individual buildings and trees and vary on the meter scale. Therefore, most numerical thermal comfort studies apply microscale models to limited spatial domains (commonly representing urban neighborhoods with building blocks) with resolutions on the order of 1 m and a few hours of simulation. This prevents the analysis of the impact of city-scale adaptation and/or mitigation strategies on thermal stress and comfort. To solve this problem, we develop a methodology to estimate thermal stress indicators and their subgrid variability in mesoscale models – here applied to the multilayer urban canopy parameterization BEP-BEM within the Weather Research and Forecasting (WRF) model. The new scheme (consisting of three main steps) can readily assess intra-neighborhood-scale heat stress distributions across whole cities and for timescales of minutes to years. The first key component of the approach is the estimation of MRT in several locations within streets for different street orientations. Second, mean wind speed and its subgrid variability are downscaled as a function of the local urban morphology based on relations derived from a set of microscale LES and RANS simulations across a wide range of realistic and idealized urban morphologies. Lastly, we compute the distributions of two thermal stress indices for each grid square, combining all the subgrid values of MRT, wind speed, air temperature, and absolute humidity. From these distributions, we quantify the high and low tails of the heat stress distribution in each grid square across the city, representing the thermal diversity experienced in street canyons. In this contribution, we present the core methodology as well as simulation results for Madrid (Spain), which illustrate strong differences between heat stress indices and common heat metrics like air or surface temperature both across the city and over the diurnal cycle.

Geology↗

DOE OSTI2025

Modeling commercial-scale CO 2 storage in the gas hydrate stability zone with PFLOTRAN v6.0

Abstract. Safe and secure carbon dioxide (CO2) storage is likely to be critical for mitigating some of the most dangerous effects of climate change. In the last decade, there has been a significant increase in activity associated with reservoir characterization and site selection for large-scale CO2 storage projects across the globe. These prospective storage sites tend to be selected for their optimal structural, petrophysical, and geochemical trapping potential. However, it has also been suggested that storing CO2 in reservoirs within the CO2 hydrate stability zone (GHSZ), characterized by high pressures and low temperatures (e.g., Arctic or marine environments), could provide a natural thermodynamic barrier to gas leakage. Evaluating the prospect of commercial-scale, long-term storage of CO2 in the GHSZ requires reservoir-scale modeling capabilities designed to account for the unique physics and thermodynamics associated with these systems. We have developed the HYDRATE flow mode and the accompanying fully implicit parallel well model in the massively parallel subsurface flow and reactive transport simulator PFLOTRAN to model CO2 injection into the marine GHSZ. We have applied these capabilities to a set of CO2 injection scenarios designed to reveal the challenges and opportunities for commercial-scale CO2 storage in the GHSZ.

carbon storage↗