Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial statistics”

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

Data and scripts associated with a manuscript on a meta-analysis synthesizing stream biogeochemical response to wildfires across space and time (v2)

This data package is associated with the publication “Catchment characteristics modulate the influence of wildfires on nitrate and dissolved organic carbon in lotic systems across space and time: A meta-analysis” submitted to Global Biogeochemical Cycles (Cavaiani et al. 2025). This study uses meta-analytical techniques to evaluate the effect of wildfire on in-stream responses in burned and unburned watersheds. The study aims to provide additional insight into the range of responses and net influences that wildfires have on hydro-biogeochemistry across broad spatial scales, burn extents, and the persistence of water-quality change. This study compiles data and metadata from 18 total publications that includes 1) surface water geochemistry data (dissolved organic carbon; nitrate), 2) climate classifications, 3) year of the wildfire, 4) the time lag between when the fire occurred and when the sampling occurred, and 5) study design of the publication. In total, this meta-analysis draws data that spans 8 climate guilds, 3 biomes, 62 watersheds, and 20 unique wildfires. See Sites_meta_data.csv for citations of the papers used in this meta-analysis. All R scripts and the associated data can also be found on GitHub at This data package was originally published in March 2024. It was updated in April 2025 (v2; new and modified files). See the change history section in the readme for more details. This data package contains five primary folders that include the following: (1) inputs; (2) output for analysis; (3) initial plots; (4) R scripts; and (5) GIS data. The data package also contains a data dictionary (dd) that provides column header definitions and a file-level metadata (flmd) file that describes every file. The “inputs” folder contains a list of all publications identified during the formal web search and an indication of whether each publication was included in the final analysis. Additionally, it includes site-level metadata, catchment characteristics, and GIS data for all publications included in the final analysis. The “Output_for_analysis” folder contains all data frames and figures generated from each R script used for additional data analysis. The “initial_plots” folder includes all exploratory figures that will be included in a supplemental and figures that will be submitted with the manuscript for publication. The “R_scripts” folder contains the scripts that perform all the data manipulations, statistical analyses, and plots. The “gis_data” folder includes shape files for each fire included in this meta-analysis. This data package contains the following file types: csv, pdf, jpeg, cpg, dbf, prj, shp, shp.ea.iso.xml, shp.iso.xml, shx.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The role of peer review and science for advancing public understanding

We live in a time of rapid change. The Covid-19 pandemic has united the world to understand and develop solutions to protect the health and wellbeing of society. At the same time, access to information on the internet and social media has exploded. People are exposed to information with many claims. Some of the information is relevant and accurate. A significant fraction may be inaccurate, biased, and potentially dangerous. With Covid-19, the information might point to cures that are toxic or risks that are not relevant. The scientific community has developed a process to develop and convey knowledge. The validation of results and their interpretation is critical to advancing knowledge. The core of the scientific process is based on developing a hypothesis-driven concept, designing and executing an experimental approach to prove or disprove the hypothesis, conducting a statistical analysis of the results to confirm their validity, and interpreting their meaning and impact. When their work is completed the scientist prepares to publish the work. Publication requires peer review. Hypotheses, experimental approaches, results, and interpretation of results are independently and anonymously reviewed by a scientific peer. The reviewers are selected based on their expertise in the field. The peer review process is essential for the integrity of the scientific process. Without a robust scientific process, society cannot discriminate between validated results and interpretations and biased, inaccurate, and potentially dangerous conjecture. With inaccurate science society and policy makers may make poor decisions that could hurt people and the world around them. While scientists should not censor information available to the public, scientists have a responsibility to inform the public. We must articulate what the scientific process is and why it is essential when considering the validity of information. Scientist are responsible for pointing out when public information has not been vetted through a scientific process and therefore cannot be considered conclusive. Scientists are responsible for highlighting information that has been generated and vetted with the scientific process. Over time, our goal is for the public to value the scientific process based on peer review. The public could become enlightened to the necessity for peer review to support decisions and action. Scientists need to contribute to advancing public knowledge. When invited to serve as a peer reviewer, scientists have a responsibility to contribute to the process. We all are extremely busy making peer review a distraction. While peer review does not directly benefit our own work and career, it benefits our chosen field of science. This is very important for the research area that our journal, Water-Energy Nexus, is focusing on. We cover the convergence of science and technology integrating across large spatial ranges (micro to global) in environmental pollution and the balance of remediation issues with energy use. Water-Energy Nexus wants to provide scientists with a platform for scientific debate and discussion on this rich and complicated topic. With the manuscripts submitted to WEN cover integrated research topics across a broad range of subjects from environment sciences through energy engineering, the role of peer review is essential. Peer review on a new submission to WEN should not be biased towards a specific science field. It should be done considering the complex nature of issues in the water-energy nexus. Scientists need to contribute to the public dialog. To accomplish these goals, we need to be honest, accurate, articulate, and use understandable language. The impact of WEN is very diverse affecting society in many ways. The scientific audience for WEN is diverse and readers may not fully understand the content of the papers outside their research expertise. Therefore, scientists, especially those who want to contribute to WEN, should work with communications specialists to convey message to benefit society. The audience includes the general public and researchers outside the author’s core specialty. While the peer-reviewed journal article is the cornerstone of our careers and critical for advancing scientific knowledge, a general interest story could benefit society even more. The author cannot lose sight that our work has two audiences, the science community, and more importantly, our family, friends, neighbors, and policy makers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Efficient and generalizable nested Fourier-DeepONet for three-dimensional geological carbon sequestration

Geological carbon sequestration (GCS) involves injecting CO2 into subsurface geological formationsfor permanent storage. Numerical simulations could guide decisions in GCS projects by predictingCO 2 migration pathways and the pressure distribution in storage formation. However, these simula-tions are often computationally expensive due to highly coupled physics and large spatial-temporalsimulation domains. Surrogate modelling with data-driven machine learning has become a promis-ing alternative to accelerate physics-based simulations. Among these, the Fourier neural operator(FNO) has been applied to three-dimensional synthetic subsurface models. Despite its good accuracyin simulating CO 2 plume migration, it requires large computational resources in training and alsolacks generalizability. Here, to further improve performance, we have developed a nested Fourier-DeepONet by combining the expressiveness of the FNO with the modularity of a deep operatornetwork (DeepONet). This new framework is twice as efficient as a nested FNO for training and has atleast 80% lower GPU memory requirement due to its flexibility to treat temporal coordinates sepa-rately. These performance improvements are achieved without compromising prediction accuracy.In addition, the generalization and extrapolation ability of nested Fourier-DeepONet beyond thetraining range has been thoroughly evaluated. Nested Fourier-DeepONet outperformed the nestedFNO for extrapolation in time with more than 50% reduced error. It also exhibited good extrapolationaccuracy beyond the training range in terms of reservoir properties, number of wells, and injectionrate.

Lee, Jonathan E. [Department of Chemical and Envir↗
DOE OSTI2021

High-Quality Revision of the Israeli Seismic Bulletin

Seismic bulletins, with trustworthy phase picks, origin times, and source locations are key for regional seismic studies, such as travel-time (TT) tomography, attenuation tomography, and anisotropy studies. To lay the groundwork for such studies in Israel, we revised the seismic bulletin of Israel and the surrounding area and obtained a trustworthy TT data set. From the earthquake and explosion bulletins of the Geophysical Institute of Israel, we compiled a starting data set of about 123,000 earthquakes and explosions that occurred during the past 40 yr. After screening out the poorly recorded events, we were left with a data set of ~38,000 well-recorded events. We then revised the remaining data set in two consecutive steps. In the first, we reviewed and updated station metadata, including changes in station metadata parameters over time. In the second step, we jointly relocated a list of selected seismic events, using the Bayesian hierarchical location software package (BayesLoc) of Myers et al. (2007) that performs joint relocation of multiple events. We observed striking dissimilarities between the spatial distributions of the newly relocated catalog and the initial locations. Although the depth distribution of the starting catalog is trimodal with peaks at 0, 5, and 10 km, the distribution in this study is unimodal, with a broad peak between 7.5 and 12.5 km. By differencing the observed arrival times and the origin times obtained through relocation with BayesLoc, we obtained a revised TT database that consists of 261,336 Pg, 132,876 Pn, 114,816 Sg, and 60,394 Sn arrivals, from a set of 30,458 jointly relocated seismic sources. In this work, we compared prerevision and postrevision TTs as a function of epicentral distance and concluded that the revised data set contains far fewer outliers and inconsistencies than the original data set. The revised TT data set may be used for seismic studies, such as TT tomography, attenuation tomography, and anisotropy studies.

58 GEOSCIENCES↗
DOE OSTI2026

Analyzing historical snow trends in interior Alaska

Study region The Chena River watershed in Interior Alaska, USA Study focus This study examines 40 years (water years 1982–2021) of snowpack characteristics to consider its hydrological implications in the 5350 km² Chena River basin. Using observations and a fine-scale physics model, we analyzed trends of snow water equivalent (SWE), snow onset and disappearance, and snow cover duration (SCD). New hydrological insights for the region Results indicate a decline in SWE across the modeled domain, averaging a decrease of 3 mm per decade, with larger decreases (up to 10 mm per decade) at lower elevations. While domain-averaged SWE trends were not statistically significant, observed SCD showed statistically significant decreases: −5.2, −5.0, and −4.4 days per decade at Teuchet Creek, Fairbanks F.O., and Little Chena Ridge, respectively. Notably, observations at SNOTEL stations and modeling revealed no statistically significant change in domain-averaged Rain-on-Snow (ROS) events over the 40-year period, contrasting some regional future estimates of increased ROS frequency. Peak streamflow did not consistently correlate with peak SWE levels, suggesting that other environmental factors such as ROS events and rapid temperature increases (e.g., a 10°C spike observed in 1992) are key drivers of hydrological outcomes. These findings improve understanding of complex subarctic hydrological processes impacting permafrost and highlight the need for adaptive water resource management to mitigate multi-factor risks like flooding and wildfire, requiring proactive planning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Data-centric framework for crystal structure identification in atomistic simulations using machine learning

Atomic-level modeling performed at large scales enables the investigation of mesoscale materials properties with atom-by-atom resolution. The spatial complexity of such cross-scale simulations renders them unsuitable for simple human visual inspection. Instead, specialized structure characterization techniques are required to aid interpretation. These have historically been challenging to construct, requiring significant intuition and effort. Here we propose an alternative framework for a fundamental structural characterization task: classifying atoms according to the crystal structure to which they belong. Our approach is data-centric and favors the employment of Machine Learning over heuristic rules of classification. A group of data-science tools and simple local descriptors of atomic structure are employed together with an efficient synthetic training set. We also introduce the first standard and publicly available benchmark data set for evaluation of algorithms for crystal-structure classification. Further, it is demonstrated that our data-centric framework outperforms all of the most popular heuristic methods—especially at high temperatures when lattices are the most distorted—while introducing a systematic route for generalization to new crystal structures. Moreover, through the use of outlier detection algorithms our approach is capable of discerning between amorphous atomic motifs (i.e., noncrystalline phases) and unknown crystal structures, making it uniquely suited for exploratory materials synthesis simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Initial Mobility Analysis for ORNL VA-EDH Synthetic Populations

Travel burdens are a major barrier to healthcare access among US Veteran patient populations, particularly those residing in rural areas. Spatial accessibility to points of care for US Veteran populations is commonly assessed in two ways. The first approach uses open data from the US Census to represent collective travel burdens, for example the distance between population-weighted census tract centroids and VHA points of care. The second approach uses restricted-access VHA patient data to measure travel costs (e.g., distance, time) for accessing points of care with respect to geolocated patient addresses and real or approximated transportation networks. While the advantage of the open data approach lies in its reproducibility, it has notable limitations in its tendency to infer individual travel behavior from aggregate population characteristics, a problem known as ecological fallacy. Conversely, while the patient data approach is able to account for individual travel behavior, its ability to account for localized access disparities (e.g., a neighborhood with exceptionally high transportation costs) and patient demographics is limited as protecting individual patient data requires their storage in closed systems with limited capacity for adequately modeling real-world travel patterns or for supplementing patient attributes. Additionally, the patient data approach cannot account for veterans who are not enrolled in the VHA system but who may be eligible for care. These challenges limit the ability to perform “what if” analyses on the effects of place-specific interventions on veteran populations with high access barriers to healthcare. To address these challenges, we explore the application of realistic synthetic populations to examine travel burdens and spatial accessibility issues among veteran patient populations. Synthetic populations provide a virtual, individually-resolved and cross-sectional representation of the veteran patient population that enables investigation of spatial access to points of care in ways in which aggregate data and patient data do not. First, synthetic populations allow one to directly assess how individuals access points of care, from synthesized residential locations to outpatient facilities on real-world transportation networks. Modeling access to points of care at the individual scale addresses the ecological fallacy problem associated with using aggregated census data to represent veteran populations and patterns of movement. Second, synthetic populations provide a means of completely representing an area’s veteran population using only publicly available, anonymized census microdata from the American Community Survey (ACS) to ensure the privacy of real-world individuals. Generating synthetic populations from the ACS also expands descriptive characteristics beyond what patient data typically offers to include socio-demographic, economic, housing, and mobility attributes. More detailed profiles of both VHA patient populations and veterans not enrolled in the VA system will provide a comprehensive picture of groups that may benefit from interventions or outreach. As an initial exercise for using synthetic populations to measure veteran travel burdens to VA care, we apply Oak Ridge National Laboratory’s (ORNL) UrbanPop capability to generate a series of synthetic VHA patient populations for 9 Veterans Integrated Services Networks (VISN) market areas in 9 Census Divisions across the continental United States, which are listed in Table 1. We use UrbanPop to produce synthetic populations for the VISN markets selected for each US Census Division, then assign VA outpatient clinic destinations to synthetic VHA patients based on travel about each VISN market’s road network. To demonstrate using the synthetic populations to evaluate healthcare travel burdens, we compare the time-based impedance between simulated home locations and VA outpatient clinics in each VISN market. We then perform validation exercises on the synthetic populations with respect to neighborhood (block group) demographic composition as well as patient mobility, comparing aggregate origin-destination statistics for the synthetic population to outpatient visits available in restricted patient data from the VA’s Corporate Data Warehouse (CDW) database.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Pooling Data Improves Multimodel IDF Estimates over Median-Based IDF Estimates: Analysis over the Susquehanna and Florida

Traditional multimodel methods for estimating future changes in precipitation intensity, duration, and frequency (IDF) curves rely on mean or median of models’ IDF estimates. Such multimodel estimates are impaired by large estimation uncertainty, shadowing their efficacy in planning efforts. Here, assuming that each climate model is one representation of the underlying data generating process, i.e., the Earth system, we propose a novel extension of current methods through pooling model data: (i) evaluate performance of climate models in simulating the spatial and temporal variability of the observed annual maximum precipitation (AMP), (ii) bias-correct and pool historical and future AMP data of reasonably performing models, and (iii) compute IDF estimates in a nonstationary framework from pooled historical and future model data. Pooling enhances fitting of the extreme value distribution to the data and assumes that data from reasonably performing models represent samples from the “true” underlying data generating distribution. Through Monte Carlo simulations with synthetic data, we show that return periods derived from pooled data have smaller biases and lesser uncertainty than those derived from ensembles of individual model data. We apply this method to NA-CORDEX models to estimate changes in 24-h precipitation intensity–frequency (PIF) estimates over the Susquehanna watershed and Florida peninsula. Our approach identifies significant future changes at more stations compared to median-based PIF estimates. The analysis suggests that almost all stations over the Susquehanna and at least two-thirds of the stations over the Florida peninsula will observe significant increases in 24-h precipitation for 2–100-yr return periods.

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

A theory of localized excitations in supercooled liquids

Here, a new connection between the structure and dynamics in glass-forming liquids is presented. We show how the origin of spatially localized excitations, as defined by the dynamical facilitation (DF) theory, can be understood from a structure-based framework. This framework is constructed by associating excitation events in the DF theory to hopping events between energy minima in the potential energy landscape (PEL). By reducing the PEL to an equal energy well picture and applying a harmonic approximation, we develop a field theory to describe elastic fluctuations about inherent states, which are energy minimizing configurations of the PEL. We model an excitation as a shear transformation zone (STZ) inducing a localized pure shear deformation onto an inherent state. We connect STZs to T1 transition events that break the elastic bonds holding the local structure of an inherent state. A formula for the excitation energy barrier, denoted as J σ , is obtained as a function of inherent-state elastic moduli and the radial distribution function. The energy barrier from the current theory is compared to the one predicted by the DF theory where good agreement is found in various two-dimensional continuous poly-disperse atomistic models of glass formers. These results strengthen the role of structure and elasticity in driving glassy dynamics through the creation and relaxation of localized excitations.

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

A review on recent machine learning applications for imaging mass spectrometry studies

Imaging mass spectrometry (IMS) is a powerful analytical technique widely used in biology, chemistry, and materials science fields that continue to expand. IMS provides a qualitative compositional analysis and spatial mapping with high chemical specificity. The spatial mapping information can be 2D or 3D depending on the analysis technique employed. Due to the combination of complex mass spectra coupled with spatial information, large high-dimensional datasets (hyperspectral) are often produced. Therefore, the use of automated computational methods for an exploratory analysis is highly beneficial. The fast-paced development of artificial intelligence (AI) and machine learning (ML) tools has received significant attention in recent years. These tools, in principle, can enable the unification of data collection and analysis into a single pipeline to make sampling and analysis decisions on the go. There are various ML approaches that have been applied to IMS data over the last decade. Here, in this review, we discuss recent examples of the common unsupervised (principal component analysis, non-negative matrix factorization, k-means clustering, uniform manifold approximation and projection), supervised (random forest, logistic regression, XGboost, support vector machine), and other methods applied to various IMS datasets in the past five years. The information from this review will be useful for specialists from both IMS and ML fields since it summarizes current and representative studies of computational ML-based exploratory methods for IMS.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Clustering of CODEX clusters

The clustering of galaxy clusters links the spatial nonuniformity of dark matter halos to the growth of the primordial spectrum of perturbations. The amplitude of the clustering signal is widely used to estimate the halo mass of astrophysical objects. The advent of cluster mass calibrations enables using clustering in cosmological studies. Aims. We analyze the autocorrelation function of a large contiguous sample of galaxy clusters, the Constrain Dark Energy with X-ray (CODEX) sample, in which we take particular care of cluster definition. These clusters were X-ray selected using the ROentgen SATellite All-Sky Survey and then identified as galaxy clusters using the code redMaPPer run on the photometry of the Sloan Digital Sky Survey. We develop methods for precisely accounting for the sample selection effects on the clustering and demonstrate their robustness using numerical simulations. Methods. Using the clean CODEX sample, which was obtained by applying a redshift-dependent richness selection, we computed the two-point autocorrelation function of galaxy clusters in the 0.1 < z < 0.3 and 0.3 < z < 0.5 redshift bins. We compared the bias in the measured correlation function with values obtained in numerical simulations using a similar cluster mass range. Results. By fitting a power law, we measured a correlation length r 0 = 18.7 ± 1.1 and slope γ = 1.98 ± 0.14 for the correlation function in the full redshift range. By fixing the other cosmological parameters to their nine-year Wilkinson Microwave Anisotropy Probe values, we reproduced the observed shape of the correlation function under the following cosmological conditions: Ω m0 = $0.22_{-0.03}^{+0.04}$ and S 8 = σ 8 (Ω m0 /0.3)0.5 = $0.85_{-0.08}^{+0.10}$ with estimated additional systematic errors of σ Ω m0 = 0.02 and σS 8 = 0.20. We illustrate the complementarity of clustering constraints by combining them with CODEX cosmological constraints based on the X-ray luminosity function, deriving Ωm = 0.25 ± 0.01 and σ 8 = $0.81_{-0.02}^{+0.01}$ with an estimated additional systematic error of σ Ω m0 = 0.07 and σ σ8 = 0.04. The mass calibration and statistical quality of the mass tracers are the dominant source of uncertainty.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Reconstructing large-scale temperature profiles around z ∼ 6 quasars

ABSTRACT High-redshift quasars ionize He ii into He iii around them, heating the intergalactic medium in the process and creating large regions with elevated temperature. In this work, we demonstrate a method based on a convolutional neural network (CNN) to recover the spatial profile for T0, the temperature at the mean cosmic density, in quasar proximity zones. We train the neural network with synthetic spectra drawn from a Cosmic Reionization on Computers simulation. We discover that the simple CNN is able to recover the temperature profile with an accuracy of ≈1400 K in an idealized case of negligible observational uncertainties. We test the robustness of the CNN and discover that it is robust against the uncertainties in quasar host halo mass, quasar continuum, and ionizing flux. We also find that the CNN has good generality with regard to the hardness of quasar spectra. This shows that with noiseless spectra, one could use a simple CNN to distinguish gas inside or outside the He iii region created by the quasar. Because the size of the He iii region is closely related to the total quasar lifetime, this method has great potential in constraining the quasar lifetime on ∼Myr time-scales. However, noise poses a big problem for accuracy and could downgrade the accuracy to ≈2340 K even for very high signal-to-noise (≳50) spectra. Future studies are needed to reduce the error associated with noise to constrain the lifetimes of reionization epoch quasars with currently available data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Effect of dust on rainfall over the Red Sea coast based on WRF-Chem model simulations

Water is the single most important element of life. Rainfall plays an important role in the spatial and temporal distribution of this precious natural resource, and it has a direct impact on agricultural production, daily life activities, and human health. One of the important elements that govern rainfall formation and distribution is atmospheric aerosol, which also affects the Earth's radiation balance and climate. Therefore, understanding how dust compositions and distributions affect the regional rainfall pattern is crucial, particularly in regions with high atmospheric dust loads such as the Middle East. Although aerosol and rainfall research has garnered increasing attention as both an independent and interdisciplinary topic in the last few decades, the details of various direct and indirect pathways by which dust affects rainfall are not yet fully understood. Here, we explored the effects of dust on rainfall formation and distribution as well as the physical mechanisms that govern these phenomena, using high-resolution WRF-Chem simulations (~1.5 km × 1.5 km) configured with an advanced double-moment cloud microphysics scheme coupled with a sectional eight-bin aerosol scheme. Our model-simulated results were realistic, as evaluated from multiple perspectives including vertical profiles of aerosol concentrations, aerosol size distributions, vertical profiles of air temperature, diurnal wind cycles, and spatio-temporal rainfall patterns. Rainfall over the Red Sea coast is mainly caused by warm rain processes, which are typically confined within a height of ~6 km over the Sarawat mountains and exhibit a strong diurnal cycle that peaks in the evening at approximately 18:00 local time under the influence of sea breezes. Numerical experiments indicated that dust could both suppress or enhance rainfall. The effect of dust on rainfall was calculated as total, indirect, and direct effects, based on 10-year August-average daily-accumulated rainfall over the study domain covering the eastern Red Sea coast. For extreme rainfall events (domain-average daily-accumulated rainfall of ≥ 1.33 mm), the net effect of dust on rainfall was positive or enhancement (6.05 %), with the indirect effect (4.54 %) and direct effect (1.51 %) both causing rainfall increase. At a 5 % significance level, the total and indirect effects were statistically significant whereas the direct effect was not. For normal rainfall events (domain-average daily-accumulated rainfall < 1.33 mm), the indirect effect enhanced rainfall (4.76 %) whereas the direct effect suppressed rainfall (-5.78 %), resulting in a negative net suppressing effect (-1.02 %), all of which were statistically significant. We investigated the possible physical mechanisms of the effects and found that the rainfall suppression by dust direct effects was mainly caused by the scattering of solar radiation by dust. The surface cooling induced by dust weakens the sea breeze circulation, which decreases the associated landward moisture transport, ultimately suppressing rainfall. For extreme rainfall events, dust causes net rainfall enhancement through indirect effects as the high dust concentration facilitates raindrops to grow when the water vapor is sufficiently available. Our results have broader scientific and environmental implications. Specifically, although dust is considered a problem from an air quality perspective, our results highlight the important role of dust on sea breeze circulation and associated rainfall over the Red Sea coastal regions. Our results also have implications for cloud seeding and water resource management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Dissect two-halo galactic conformity effect for central galaxies: the dependence of star formation activities on the large-scale environment

We investigate the two-halo galactic conformity effect for central galaxies, which is the spatial correlation of the star formation activities for central galaxies to several Mpcs, by studying the dependence of the star formation activities of central galaxies on their large-scale structure in our local Universe using the SDSS data. Here we adopt a novel environment metric using only central galaxies quantified by the distance to the nth nearest central galaxy. This metric measures the environment within an aperture from ∼1 to ≳ 10 Mpc, with a median value of ∼4 Mpc. We found that two kinds of conformity effects in our local Universe. The first one is that low-mass central galaxies are more quenched in high-density regions, and we found that this effect mainly comes from low-mass centrals that are close to a more massive halo. A similar trend is also found in the IllustrisTNG simulation, which can be entirely explained by backsplash galaxies. The second conformity effect is that massive central galaxies in low-density regions are more star-forming. This population of galaxies also possesses a higher fraction of spiral morphology and lower central stellar velocity dispersion, suggesting that their low quiescent fraction is due to less-frequent major merger events experienced in the low-density regions and, as a consequence, less-massive bulges and central black holes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Field-level reconstruction from foreground-contaminated 21-cm maps

Current and upcoming 21-cm experiments will soon be able to map 21-cm spatial fluctuations in three dimensions for a wide range of redshifts. However, bright foreground contamination and the nature of radio interferometry create significant challenges, making it difficult to access rich cosmological information from the Fourier modes that lie within the “foreground wedge”. Here, in this work, we introduce two approaches aiming to reconstruct the full 21-cm density field, including the missing modes in the wedge: (a) a field-level inference under an effective field theory (EFT) framework; (b) a diffusion-based deep generative model trained on simulations. Under the EFT framework, we implement a fully differentiable forward model that maps the initial conditions of matter fluctuations to the observed, foreground-filtered 21-cm maps. This enables a gradient-based sampler to simultaneously sample the initial conditions and bias parameters, allowing a physically motivated mode reconstruction. Alternatively, we apply a variational diffusion model to perform 21-cm density reconstruction at the map level. Our model is trained on semi-numerical simulations over a wide range of astrophysical parameters. Our results from both approaches should provide improved cosmological constraints from the field level and also enable cross-correlation between experiments that have little or no overlapping modes.

cosmological perturbation theory↗
DOE OSTI2025

Incorporating the molecular-scale into a hydrodynamic description of confined aqueous systems

Hydrodynamics provides a continuum-level description of fluid motion, but its applicability at the nanoscale becomes uncertain due to the emerging importance of molecular-level effects such as spatial heterogeneity. Hydrodynamic boundary conditions that incorporate molecular details allow us to partition the system into a near-wall region and a bulk fluid region. We identify a hydrodynamic wall located inside the fluid that determines where slip begins. By extending the hydrodynamic wall with the slip length, the position of the extrapolated wall is established. This offers a unified description of both slip and stagnant flow behaviors, with wall hydrophobicity characterized by the relative location of the extrapolated wall with respect to the physical wall. Employing this concept in analyses of equilibrium molecular dynamics (MD) and non-equilibrium MD simulations of Couette and Poiseuille flows, our results demonstrate consistency between equilibrium and non-equilibrium approaches across different flow types and confinement levels. This demonstrates the robust nature of linear response theory. We then explore the effects of fluid-wall and bulk fluid interactions on the hydrodynamic properties. These findings enhance the effectiveness of molecular-based simulations for investigating complex confined systems in nanofluidics, biology, and colloidal science, offering a complementary molecular-scale perspective to traditional continuum approaches.

flow simulations↗
DOE OSTI2020

A Causal Inference Model Based on Random Forests to Identify the Effect of Soil Moisture on Precipitation

Soil moisture influences precipitation mainly through its impact on land–atmosphere interactions. Understanding and correctly modeling soil moisture–precipitation (SM–P) coupling is crucial for improving weather forecasting and subseasonal to seasonal climate predictions, especially when predicting the persistence and magnitude of drought. However, the sign and spatial structure of SM–P feedback are still being debated in the climate research community, mainly due to the difficulty in establishing causal relationships and the high degree of nonlinearity in land–atmosphere processes. To this end, we developed a causal inference model based on the Granger causality analysis and a nonlinear machine learning model. This model includes three steps: nonlinear anomaly decomposition, nonlinear Granger causality analysis, and evaluation of the quality of SM–P feedback, which eliminates the nonlinear response of interannual and seasonal variability and the memory effects of climatic factors and isolates the causal relationship of local SM–P feedback. We applied this model by using National Climate Assessment–Land Data Assimilation System (NCA-LDAS) datasets over the United States. Here, the results highlight the importance of nonlinear atmosphere responses in land–atmosphere interactions. In addition, the strong feedback over the southwestern United States and the Great Plains both highlight the impacts of topographic factors rather than only the sensitivity of evapotranspiration to soil moisture. Furthermore, the SM–P index defined by our framework is used to benchmark Earth system models (ESMs), which provides a new metric for efficiently identifying potential model biases in modeling local land–atmosphere interactions and may help the development of ESMs in improving simulations of water cycle variability and extremes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The MaNGA firefly Value-Added-Catalogue: resolved stellar populations of 10,010 nearby galaxies

We present the Mapping Nearby Galaxies at Apache Point Observatory survey (MaNGA) FIREFLY Value Added Catalogue (VAC) – a catalogue of ~3.7 million spatially resolved stellar population properties across 10 010 nearby galaxies from the final data release of the MaNGA survey. The full spectral fitting code FIREFLY is employed to derive parameters such as stellar ages, metallicities, stellar and remnant masses, star formation histories, star formation rates, and dust attenuation. In addition to Voronoi-binned measurements, our VAC also provides global properties, such as central values and radial gradients. Two variants of the VAC are available: presenting the results from fits using the M11-MILES and the novel MaStar stellar population models. MaStar allows to constrain the fit over the whole MaNGA wavelength range, extends the age-metallicity parameter space, and uses empirical spectra from the same instrument as MaNGA. The fits employing MaStar models find, on average, slightly younger ages, higher mass-weighted metallicities, and smaller colour excesses. These differences are reduced when matching the wavelength range and converging template grids. We further report that FIREFLY stellar masses are systematically lower by ~0.3 dex than masses from the MaNGA PCA and Pipe3D VACs, but match masses from the NSA best with only ~0.1-dex difference. Finally, we show that FIREFLY stellar ages correlate with spectral index age indicators H δ A and D n (4000), though with a clear additional metallicity dependence.

79 ASTRONOMY AND ASTROPHYSICS↗