Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “rooting depth”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Factors and processes governing the C-14 content of carbonate in desert soils

A model is presented describing the factors and processes which determine the measured C-14 ages of soil calcium carbonate. Pedogenic carbonate forms in isotopic equilium with soil CO2. Carbon dioxide in soils is a mixture of CO2 derived from two biological sources: respiration by living plant roots and respiration of microorganisms decomposing soil humus. The relative proportion of these two CO2 sources can greatly affect the initial C-14 content of pedogenic carbonate: the greater the contribution of humus-derived CO2, the greater the initial C-14 age of the carbonate mineral. For any given mixture of CO2 sources, the steady-state (14)CO2 distribution vs. soil depth can be described by a production/diffusion model. As a soil ages, the C-14 age of soil humus increases, as does the steady-state C-14 age of soil CO2 and the initial C-14 age of any pedogenic carbonate which forms. The mean C-14 age of a complete pedogenic carbonate coating or nodule will underestimate the true age of the soil carbonate. This discrepancy increases the older a soil becomes. Partial removal of outer (and younger) carbonate coatings greatly improves the relationship between measured C-14 age and true age. Although the production/diffusion model qualitatively explains the C-14 age of pedogenic carbonate vs. soil depth in many soils, other factors, such as climate change, may contribute to the observed trends, particularily in soils older than the Holocene.

Amundson, Ronald↗
NASA NTRS2019

Assimilation of Satellite-Based Snow Cover and Freeze/Thaw Observations Over the High Mountain Asia

Toward qualifying hydrologic changes in the High Mountain Asia (HMA) region, this study explores the use of a hyper-resolution (1 km) land data assimilation (DA)framework developed within the NASA Land Information System using the Noah Multi-parameterization Land Surface Model (Noah-MP) forced by the meteorological boundary conditions from Modern-Era Retrospective analysis for Research and Applications, Version 2 data. Two different sets of DA experiments are conducted:(1) the assimilation of a satellite-derived snow cover map (MOD10A1) and (2) the assimilation of the NASA MEaSUREs landscape freeze/thaw product from 2007 to 2008. The performance of the snow cover assimilation is evaluated via comparisons with available remote sensing-based snow water equivalent product and ground-based snow depth measurements. For example, in the comparison against ground-based snow depth measurements, the majority of the stations (13 of 14) show slightly improved goodness-of-fit statistics as a result of the snow DA, but only four are statistically significant. In addition, comparisons to the satellite-based land surface temperature products (MOD11A1 and MYD11A1) show that freeze/thaw DA yields improvements (at certain grid cells) of up to 0.58 K in the root-mean-square error (RMSE) and 0.77K in the absolute bias (relative to model-only simulations). In the comparison against three ground-based soil temperature measurements along the Himalayas, the bias and the RMSE in the 0-10 cm soil temperature are reduced (on average) by 10 and 7%,respectively. The improvements in the top layer of soil estimates also propagate through the deeper soil layers, where the bias and the RMSE in the 10-40 cm soil temperature are reduced (on average) by 9 and 6%, respectively. However, no statistically significant skill differences are observed for the freeze/thaw DA system in the comparisons against ground-based surface temperature measurements at mid-to-low altitude. Therefore, the two proposed DA schemes show the potential of improving the predictability of snow mass, surface temperature, and soil temperature states across HMA, but more ground-based measurements are still required, especially at high-altitudes, in order to document a more statistically significant improvement as a result of the two DA schemes.

High Mountain Asia↗
DOE OSTI2023

Explosive Yield Estimation Using Regional Seismic Moment Tensors

Here, we use the Pasyanos and Chiang (2022) data set to calculate the seismic moment M 0 for each explosion and use the measured explosive yield W to validate the W~M 0 relationship in Denny and Johnson (1991; hereafter, DJ91). The M 0 is corrected by transforming to a potency tensor and applying more appropriate near-source geophysical parameter values in the moment estimate. The mean residual between observed and predicted yield is near zero; however, the standard deviation of the residuals results in an F-value (a 95% confidence factor) of about 5. We re-estimate the coefficients in the DJ91 model and find similar values and only a slight improvement in the F-value. Next, we embark on a similar model selection process as DJ91, allowing for non-cube-root yield scaling and other plausible near-source elastic moduli. As was found by DJ91, the yield dependence is not significantly different from unity, and a cube root assumption is valid. Therefore, we yield scale the seismic moment and test the significance of all plausible explanatory variables. Isotropic moment performs better in the response variable than total moment. The preference for isotropic moment could be due to its relationship to volume change, which would be more directly affected by explosive yield. Surprisingly, we find that the overburden pressure, which is a function of depth, is not a significant parameter in the model. We hypothesize that this is due to the competing depth effects on source asymmetry and the incorporation of depth in the Green’s functions used to calculate the seismic moment tensors. Importantly, this emphasizes that only seismic moment tensor-derived moments should be used in these models. After removing insignificant model parameters, we are left with a simple model to predict explosive yield $\widehat{W}$ in kt from isotropic moment M I in N·m, $\widehat{W}$=κ –1.4132 10 0.035626GP M I , in which κ and GP are the near-source bulk modulus and gas porosity in Pa and %, respectively. The F-value for this model is approximately 3.

58 GEOSCIENCES↗
DOE OSTI2024

Low flow characteristics for regulated and unregulated streams in North Carolina and prediction using climate signals

In this work, low flow statistics of regulated and unregulated streams in the state of North Carolina were updated with streamflow data through 2019. About 22% of the streams considered show a significant downward trend, but considerable low-frequency variability confounding trends. The relationship between regional index time series of groundwater-depth low flow and the Atlantic Multidecadal Oscillation (AMO) was examined to determine whether low-frequency climate modes can account for the long-term pattern in low flows. Consequently, a significant correlation was found between AMO and groundwater-depth low flow, such that positive AMO is associated with lower groundwater-depth low flow and vice versa, particularly for the Piedmont region. Predictive equations for annual low flows at the ecoregion level shows that springtime average streamflow and AMO were selected as the primary predictors of low flow for coastal and Piedmont regions, whereas springtime average streamflow and the November–December–January average Oceanic Niño Index were used as the primary predictors for the mountain region. The relative root mean square error (RMSE) of the disaggregated predictions to US Geologic Survey gauge locations was <23% at 79% of the stations, between 24% and 43% at 10% of the stations, and greater than 44% at 1% of the stations. The remaining 10% of stations showed large RMSEs. This latter percentage is characterized by smaller drainage basins and intermittent flows, suggesting the prediction models are not applicable to drainage basins smaller than roughly 20 square km and intermittent streams.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

Advances in Assimilation of Satellite-Based Passive Microwave Observations for Soil-Moisture Estimation

Satellite-based microwave measurements have long shown potential to provide global information about soil moisture. The European Space Agency (ESA) Soil Moisture and Ocean Salinity (SMOS, [1]) mission as well as the future National Aeronautics and Space Administration (NASA) Soil Moisture Active and Passive (SMAP, [2]) mission measure passive microwave emission at L-band frequencies, at a relatively coarse (40 km) spatial resolution. In addition, SMAP will measure active microwave signals at a higher spatial resolution (3 km). These new L-band missions have a greater sensing depth (of -5cm) compared with past and present C- and X-band microwave sensors. ESA currently also disseminates retrievals of SMOS surface soil moisture that are derived from SMOS brightness temperature observations and ancillary data. In this research, we address two major challenges with the assimilation of recent/future satellite-based microwave measurements: (i) assimilation of soil moisture retrievals versus brightness temperatures for surface and root-zone soil moisture estimation and (ii) scale-mismatches between satellite observations, models and in situ validation data.

De Lannoy, Gabrielle J. M.↗
DOE OSTI2021

Diverse ecophysiological adaptations of subsurface Thaumarchaeota in floodplain sediments revealed through genome-resolved metagenomics

Abstract The terrestrial subsurface microbiome contains vastly underexplored phylogenetic diversity and metabolic novelty, with critical implications for global biogeochemical cycling. Among the key microbial inhabitants of subsurface soils and sediments are Thaumarchaeota, an archaeal phylum that encompasses ammonia-oxidizing archaea (AOA) as well as non-ammonia-oxidizing basal lineages. Thaumarchaeal ecology in terrestrial systems has been extensively characterized, particularly in the case of AOA. However, there is little knowledge on the diversity and ecophysiology of Thaumarchaeota in deeper soils, as most lineages, particularly basal groups, remain uncultivated and underexplored. Here we use genome-resolved metagenomics to examine the phylogenetic and metabolic diversity of Thaumarchaeota along a 234 cm depth profile of hydrologically variable riparian floodplain sediments in the Wind River Basin near Riverton, Wyoming. Phylogenomic analysis of the metagenome-assembled genomes (MAGs) indicates a shift in AOA population structure from the dominance of the terrestrial Nitrososphaerales lineage in the well-drained top ~100 cm of the profile to the typically marine Nitrosopumilales in deeper, moister, more energy-limited sediment layers. We also describe two deeply rooting non-AOA MAGs with numerous unexpected metabolic features, including the reductive acetyl-CoA (Wood-Ljungdahl) pathway, tetrathionate respiration, a form III RuBisCO, and the potential for extracellular electron transfer. These MAGs also harbor tungsten-containing aldehyde:ferredoxin oxidoreductase, group 4f [NiFe]-hydrogenases and a canonical heme catalase, typically not found in Thaumarchaeota. Our results suggest that hydrological variables, particularly proximity to the water table, impart a strong control on the ecophysiology of Thaumarchaeota in alluvial sediments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Experimental and Computational Studies of Stress Corrosion Cracking of Alloys 308/309 and 82/182 Weldments in Corrosive and Radiation Environment

The goal of the project was to determine factors that influence SCC and IASCC in weldments found in LWR nuclear power plants. We focused on a SA508-304L SS weldment fabricated by EPRI using gas tungsten arc welding and used an aggressive BWR normal water chemistry (NWC) immersion environment. This weldment used 309L butter and 308L groove filler material. The microstructure of the 309L butter was non-uniform, exhibiting a 20–30μm thick martensitic layer closest to the SA508 interface, a 1–4 mm thick single γ austenite phase dilution zone, and a γ–δ duplex region extending to the 308L groove filler. The 308L groove filler had an entirely γ–δ duplex microstructure. Two approximately 1-inch thick 304L and SA508 plates approximately 12 by 6 square inches were joined using standard nuclear grade welding techniques. This included a post weld heat treatment of the 309L butter after application, a 0.32 cm fit-up root opening, and 57 bead lines of 308L groove filler applied in 18 layers. A 60 degree weld bevel angle was used and the 309L butter was 1.5 cm thick. Displacement cascade damage was induced using proton irradiation at the Michigan Ion Beam Laboratory. The incident proton energy was 2 MeV, the sample temperature was 360 ºC, and the calculated dpa value at 10 μm (60% of the Bragg peak depth of ~18 μm) was 5 dpa using the quick Kinchin-Pease model. Proton irradiation to this damage level required approximately 125 hours of beam time. Two types of samples were irradiated, tensile specimens and TEM bars. Samples were selected from all regions of the weldment, including the SA508-309L butter interface, the 309L-308L interface, and 308L-304L interface. The stainless steel alloys (304L, 308L, and 309L) within the heat affected zone are characterized by a duplex skeletal morphology of δ-ferrite and γ-austenite resulting from the recrystallization associated with weld fabrication. Approximately 7 to 8 mm of length along the specimen was irradiated. The gauge volume surfaces were mechanically polished and then electro-polished to remove mechanical damage from the mechanical polishing step prior to irradiation. Immersion tests were performed in a recirculating autoclave under BWR NWC conditions (2000 ppb wt. dissolved oxygen, neutral pH, 288 ºC, 10 MPa, and inlet water conductivity <100 nS/cm) to accelerate corrosion. Constant strain rate tests were performed either to failure or to approximately 5% strain. Strain rates of 10 -7 to 10 -6 mm/mm/s were used and typical immersion testing required four to six weeks to achieve failure or strains near 5%. Analysis primarily used advanced electron microscopy techniques of FIB lift out specimens.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Aboveground Rather Than Belowground Productivity Drives Variability in Miscanthus × giganteus Net Primary Productivity

Quantifying the carbon (C) uptake of Miscanthus × giganteus ( M × g ) in both aboveground and belowground structures (e.g., net primary productivity (NPP)) and differences among methodological approaches is crucial. Our objectives were to directly measure Mxg NPP and evaluate the effects of nitrogen application, location, and belowground biomass sampling methods. We hypothesize that increased nitrogen application increases the overall NPP of M × g and that quantifying rhizome biomass using excavations will produce the lowest variability between replicates. We collected biomass from mature M × g stands from three locations in Iowa with three nitrogen application rates and one site in Illinois. We destructively sampled at two time points, when rhizome mass is anticipated to be at a minimum (initial) and anticipated to be at its maximum (peak). Biomass was collected from 1 × 1 m quadrats in which one in-clump and one beside-clump cores were collected and then excavated to 30 cm depth to extract all rhizomes. We found that aboveground M × g NPP ranged from 15.4 Mg DM ha –1 year –1 to 36.4 Mg DM ha –1 year–1 and belowground M × g NPP ranged from 4.4 Mg DM ha –1 year –1 to 19.6 Mg DM ha –1 year –1 . M × g NPP varied across sites, fertilization, and calculation assumptions. Aboveground NPP (yield) was on average 68.7% of the total NPP. Root-to-shoot ratios at peak biomass decreased with nitrogen application rate, from an average of 1.9 for 0 N plots to 0.89 for 224 N fertilized plots. There was more variation in core data than from excavations; however, when in-clump and beside-clump cores were averaged together, core and excavation averages were not different. Overall, these results show that the range of mature M × g NPP is driven by aboveground productivity, influenced by nitrogen application and site. Our results provide useful data to constrain agro-ecosystem models and provide crucial insights for future perennial belowground sampling.

Hartman, Theodore [Univ. of Illinois at Urbana-Cha↗
DOE OSTI2023

FUN-BioCROP model with litter decomposition parameters derived from the LIDET dataset

This repository contains the code and data necessary to run the FUN-BioCROP (Fixation and Uptake of Nitrogen-Bioenergy Carbon, Rhizosphere, Organisms, and Protection) model with litter decomposition parameters derived from a modified Monte Carlo simulation that used the Long-term Intersite Decomposition Experiment Team dataset. Related publication:Juice, S.M., Ridgeway, J.R., Hartman, M.D., Parton, W.J., Berardi, D.M., Sulman, B.N., Allen, K.E., & Brzostek, E.R. Reparameterizing litter decomposition using a simplified Monte Carlo method improves litter decay simulated by a microbial model and alters bioenergy soil carbon estimates. Description of Files: FUNBioCROP_LIDET Study.Rmd R code with FUN-BioCROP model that can be run with 10 different sets of parameters for litter decomposition (Baseline, LIDET, or eight other best parameter sets identified in the modified Monte Carlo simulation. CORPSE Functions_Bioenergy_V2.R Code with CORPSE model functions, called by FUNBioCROP_LIDET Study.Rmd Model Input Data: bulk.csv, bulk_till.csv, rhizo.csv, rhizo_till.csv, litter.csv Initial C and N (kg C or N/m2) pool values for each soil compartment, final values from spin up. All five files have the same columns: (Column - Description - Units) uFastC - Unprotected fast decomposing carbon - kg carbon/m2 uSlowC - Unprotected slow decomposing carbon - kg carbon/m2 uNecroC - Unprotected necromass carbon - kg carbon/m2 pFastC - Protected fast decomposing carbon - kg carbon/m2 pSlowC - Protected slow decomposing carbon - kg carbon/m2 pNecroC - Protected necromass carbon - kg carbon/m2 livingMicrobeC - Carbon in living microbial biomass - kg carbon/m2 uFastN - Unprotected fast decomposing nitrogen - kg nitrogen/m2 uSlowN - Unprotected slow decomposing nitrogen - kg nitrogen/m2 uNecroN - Unprotected necromass nitrogen - kg nitrogen/m2 pFastN - Protected fast decomposing nitrogen - kg nitrogen/m2 pSlowN - Protected slow decomposing nitrogen - kg nitrogen/m2 pNecroN - Protected necromass nitrogen - kg nitrogen/m2 inorganicN - Inorganic nitrogen - kg nitrogen/m2 CO2 - Carbon in carbon dioxide - kg carbon/m2 livingMicrobeN - Nitrogen in living microbial biomass - kg nitrogen/m2 Model Input Data: FluxTower_AvgSoilT.csv: Average daily soil temperature (oC) at 10 cm depth at University of Illinois Urbana-Champaign (UIUC) Energy Farm flux tower from 7/2008-3/2016. (One year of averaged data) Model Input Data: FluxTower_AvgSoilVWC.csv: Average daily soil volumetric water content (VWC) at 10 cm depth at UIUC Energy Farm flux tower from 7/2008-3/2016. (One year of averaged data) Model Input Data: input_CCS_LIDET Study.csv: This file has daily data to run FUN-BioCROP (Column - Description - Units): yr - calendar year - year doy - day of year (1 to 365) (no leap year) - day anpp - aboveground NPP (DayCent) - kg C/m2/day bnpp - belowground NPP (DayCent) - kg C/m2/day aglivc - live aboveground biomass carbon (DayCent) - kg C/m2 bglivcj - live juvenile fine root biomass carbon (DayCent) - kg C/m2 bglivcm - live mature fine root biomass carbon (DayCent) - kg C/m2 aglivn - live aboveground biomass nitrogen (DayCent) - kg N/m2 bglivnj - live juvenile fine root biomass nitrogen (DayCent) - kg N/m2 bglivnm - live mature fine root biomass nitrogen (DayCent) - kg N/m2 nyr - simulation year - year cult - indicates a cultivation event (0 or 1) crop - indicates a new crop (0 or 1) fert - indicates a fertilizer event (0 or 1) frst - indicates the first day of the growing season (0 or 1) harv - indicates a harvest event (0 or 1) last - indicates the end of the growing season (0 or 1) croptype - crop type (0=none; 1=alfalfa; 2=corn; 3=grass clover pasture; 4=soybean; 5=wheat) cropsrl - crop specific root length - mm/g root cultrhizmix - fraction of rhizosphere mixed with bulk soil during cultivation (0.0-1.0) - fraction cultlitmix - fraction of litter mixed with bulk soil during cultivation (0.0-1.0) - fraction harvremov - fraction of above ground biomass removed during harvest (0.0-1.0) - fraction fertamt - fertilization amount - g N/m2 lifehist - plant life history (0 = annual, 1 = perennial) froot_turnover_c - amount of C in fine root turnover - kg C/m2 froot_turnover_n - amount of N in fine root turnover - kg N/m2 agrd_turnover_c - amount of C in aboveground biomass turnover - kg C/m2 agrd_turnover_n - amount of N in aboveground biomass turnover - kg N/m2 leaf_litter_fastfrac - Fast decomposing fraction of leaf litter (0.0-1.0) - fraction root_litter_fastfrac - Fast decomposing fraction of root litter (0.0-1.0) - fraction root_diameter - root diameter - mm root_length - root length - mm root/m2 rhizo_frac - fraction of total soil volume that is rhizosphere (0.0 - 1.0) - fraction date - date in format YYYY-MM-DD Instructions: Save the model code ("FUN-BioCROP_LIDET Study.Rmd") and accompanything files (data streams and CORPSE function code) in the same folder. In model code "Chunk 3: Load CORPSE Data Streams" set the working directory (setwd) to the folder with the files saved in step #1. In "Chunk 5: Define LIDET parameter sets" select the litter decomposition parameter set to be used in the run, and comment out all other sets. If changing any parameter values, edit them in "Chunk 6: Load parameters." Run all chunks up to and including "Chunk 10: Prepare Data for Export." In "Chunk 11: Export Output Data" edit data frames for export and filenames, as necessary. "Chunk 12: Graph Total Soil C" makes a figure of C remaining over the model run period. Description of each model chunk (in file FUN-BioCROP_LIDET Study.Rmd): Chunk 1: Remove all functions, clear memory. Removes all functions from R environment, clears the memory. Chunk 2: Load Packages. Loads packages necessary to run the code. Chunk 3: Load CORPSE Data Streams. Sets the working directory and loads the data files necessary to run CORPSE. Chunk 4: Load CORPSE Functions. Loads the R script with CORPSE functions from the working directory, "CORPSE Functions_Bioenergy_V2.R". Chunk 5: Define LIDET parameter sets. Has ten different parameter sets for litter decomposition tested in this study: Baseline parameters, LIDET parameters, and the other 8 best performing parameter sets identified in the modified Monte Carlo. To run the model, all but one parameter set must be commented out. Chunk 6: Load Parameters. Loads all fixed parameters to run the model. Data frame with definitions of parameters is in the CORPSE function script "CORPSE Functions_Bioenergy_V2.R" Chunk 7: Prepare Data Streams. Takes data streams loaded in Chunk 3 and puts them in the format necessary to run the model. The model is coded to run at least two sites at a time, so if only one site is being run it must be run in duplicate. Individual data tables of daily values are created in this chunk from the input data file. Chunk 8: Set Initial Conditions. Creates data tables of soil C and N pools for each soil compartment (rhizo_till, rhizo, bulk_till, bulk, litter) and loads initial values into the data tables. Creates lists for each soil compartment to hold model output. Chunk 9: Load FUN Data and Set Up Matrices. Uses DayCent data to calculate FUN input data: root and leaf N demand, total N demand, plant CN, leaf N available for retranslocation, and litter production. Creates matrices for FUN model outputs. Chunk 10: Run Model. Runs the model. Chunk 11: Prepare Data for Export. Combines data from each day saved as lists into data frames for each soil compartment. Adds values from all soil compartments together to calculate total soil values, creates separate data frames for each soil C and N pool (e.g., protected slow C) for the total soil value. Adds different C and N pools together to calculate total soil C and N for all layers. Creates data frame of ratio of protected to unprotected SOC. Organizes FUN data for export. Chunk 12: Export Results. Exports CSV files of model results to the working directory. Chunk 13: Graph Total Soil C. Makes figure of C remaining over time. Related Links: Original FUN-BioCROP model: https://github.com/BrzostekEcologyLab/FUN-BioCROP LIDET dataset: https://andlter.forestry.oregonstate.edu/data/abstract.aspx?dbcode=TD023

Juice, Stephanie↗
DOE OSTI2020

Degradation of REBCO coated conductors due to a combination of epoxy impregnation, thermal cycles, and quench: Characteristics and a method of alleviation

REBa2Cu3Ox (REBCO) coated conductors are an enabling technology for >1 GHz nuclear magnetic resonance spectrometers that need >23.5 T solenoids, compact high-field fusion reactors beyond the International Thermonuclear Experimental Reactor, and 20 T class accelerator magnets. An in-depth understanding of stress and strain management is critically needed for these applications. Here, we report experimental data and an analysis of damages of REBCO-coated conductors in several scenarios relevant for their applications in high-field magnets and unique to their multi-layer thin film structure. The examined scenarios include failure induced by quench, failure induced by epoxy impregnation, and failure induced by a combination of epoxy impregnation and quench. The root cause is film blistering and peeling induced by thermally interfacial stresses characteristic of multilayered heteroepitaxial thin film structures on thick substrates utilized. The quantitative stress analysis is introduced using simple and easy to use formulas and clearly indicates how materials and structural characteristics affect the magnitude and the distribution of stresses and deflections. Finally, a simple method is introduced for removing degradation of REBCO coated conductors due to epoxy impregnation. Our data and analysis help our understanding of the practical limit of REBCO-coated conductors and can be utilized as a guidance for the physical design of REBCO coated conductors and high-current cables assembled from them and for minimizing damages.

Yin, Shijian (ORCID:0000000224312274)↗
DOE OSTI2022

PRODeepSyn: predicting anticancer synergistic drug combinations by embedding cell lines with protein–protein interaction network

Abstract Although drug combinations in cancer treatment appear to be a promising therapeutic strategy with respect to monotherapy, it is arduous to discover new synergistic drug combinations due to the combinatorial explosion. Deep learning technology holds immense promise for better prediction of in vitro synergistic drug combinations for certain cell lines. In methods applying such technology, omics data are widely adopted to construct cell line features. However, biological network data are rarely considered yet, which is worthy of in-depth study. In this study, we propose a novel deep learning method, termed PRODeepSyn, for predicting anticancer synergistic drug combinations. By leveraging the Graph Convolutional Network, PRODeepSyn integrates the protein–protein interaction (PPI) network with omics data to construct low-dimensional dense embeddings for cell lines. PRODeepSyn then builds a deep neural network with the Batch Normalization mechanism to predict synergy scores using the cell line embeddings and drug features. PRODeepSyn achieves the lowest root mean square error of 15.08 and the highest Pearson correlation coefficient of 0.75, outperforming two deep learning methods and four machine learning methods. On the classification task, PRODeepSyn achieves an area under the receiver operator characteristics curve of 0.90, an area under the precision–recall curve of 0.63 and a Cohen’s Kappa of 0.53. In the ablation study, we find that using the multi-omics data and the integrated PPI network’s information both can improve the prediction results. Additionally, the case study demonstrates the consistency between PRODeepSyn and previous studies.

Wang, Xiaowen↗
DOE OSTI2021

Assimilation of a Coordinated Fleet of Uncrewed Aircraft System Observations in Complex Terrain: EnKF System Design and Preliminary Assessment

Abstract Uncrewed aircraft system (UAS) observations collected during the 2018 Lower Atmospheric Process Studies at Elevation—a Remotely Piloted Aircraft Team Experiment (LAPSE-RATE) field campaign were assimilated into a high-resolution configuration of the Weather Research and Forecasting Model using an ensemble Kalman filter. The benefit of UAS observations was assessed for a terrain-driven (drainage and upvalley) flow event that occurred within Colorado’s San Luis Valley (SLV) using independent observations. The analysis and prediction of the strength, depth, and horizontal extent of drainage flow from the Saguache Canyon and the subsequent transition to upvalley and up-canyon flow were improved relative to that obtained both without data assimilation (benchmark) and when only surface observations were assimilated. Assimilation of UAS observations greatly improved the analyses of vertical variations in temperature, relative humidity, and winds at multiple locations in the northern portion of the SLV, with reductions in both bias and the root-mean-square error of roughly 40% for each variable relative to the benchmark run. Despite these noted improvements, some biases remain that were tied to measurement error and/or the impact of the boundary layer parameterization on vertically spreading the observations, both of which require further exploration. The results presented here highlight how observations obtained with a fleet of profiling UAS improve limited-area, high-resolution analyses and short-term forecasts in complex terrain.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2015

Spatiotemporal Prediction of Fine Particulate Matter Using High-Resolution Satellite Images in the Southeastern US 2003-2011

Numerous studies have demonstrated that fine particulate matter (PM(sub 2.5), particles smaller than 2.5 micrometers in aerodynamic diameter) is associated with adverse health outcomes. The use of ground monitoring stations of PM(sub 2.5) to assess personal exposure, however, induces measurement error. Land-use regression provides spatially resolved predictions but land-use terms do not vary temporally. Meanwhile, the advent of satellite-retrieved aerosol optical depth (AOD) products have made possible to predict the spatial and temporal patterns of PM(sub 2.5) exposures. In this paper, we used AOD data with other PM(sub 2.5) variables, such as meteorological variables, land-use regression, and spatial smoothing to predict daily concentrations of PM(sub 2.5) at a 1 sq km resolution of the Southeastern United States including the seven states of Georgia, North Carolina, South Carolina, Alabama, Tennessee, Mississippi, and Florida for the years from 2003 to 2011. We divided the study area into three regions and applied separate mixed-effect models to calibrate AOD using ground PM(sub 2.5) measurements and other spatiotemporal predictors. Using 10-fold cross-validation, we obtained out of sample R2 values of 0.77, 0.81, and 0.70 with the square root of the mean squared prediction errors of 2.89, 2.51, and 2.82 cu micrograms for regions 1, 2, and 3, respectively. The slopes of the relationships between predicted PM2.5 and held out measurements were approximately 1 indicating no bias between the observed and modeled PM(sub 2.5) concentrations. Predictions can be used in epidemiological studies investigating the effects of both acute and chronic exposures to PM(sub 2.5). Our model results will also extend the existing studies on PM(sub 2.5) which have mostly focused on urban areas because of the paucity of monitors in rural areas.

PM2.5↗
DOE OSTI2020

On the oceanic origin for the enhanced seasonal cycle of SST in the midlatitudes under global warming

Climate models project an enhancement in SST seasonal cycle over the midlatitude oceans under global warming. However, the maintaining mechanism for this enhancement is still under debate. This study investigates the SST seasonal cycle enhancement under quadrupling CO2 using a set of partially coupled experiments, in which the contribution from direct CO2 forcing (without the effect of wind changes) and wind effects can be isolated. Results indicate that both the direct CO2 and total wind effects contribute to the enhancement in the SST seasonal cycle, with the former (latter) effect being more important in the Northern Hemisphere (Southern Hemisphere). Further decomposition of the wind effect into the wind stress feedback and wind-evaporation-SST (WES) feedback through the change of wind speed reveals the importance of the wind stress-driven ocean response in the change of SST seasonal cycle, a result in contrast to a previous study that ascribed the midlatitude SST seasonal cycle change to the thermodynamic WES feedback. The direct CO2 effect regulates SST seasonal cycle mainly through the mediating effects of the ocean mixed layer depth (MLD): i) the CO2 induced warming leads to MLD shoaling, thus reducing its effective heat capacity; ii) the climatological MLD is greater in winter than summer, given the same amount of energy flux from the surface, SST warming in winter should be smaller. Notably, the surface wind seasonal cycle change due solely to the direct CO2 effect is found to bear a great resemblance to the full wind response, implicating that the direct CO2 effect is the root source for the enhancement of the midlatitude SST seasonal cycle. This is further supported by an ocean-alone experiment that replicates the SST seasonal cycle enhancement under a spatially and temporally homogeneous surface thermal forcing into the ocean.

Liu, Fukai↗
NASA NTRS2016

Relationship of EUV Irradiance Coronal Dimming Slope and Depth to Coronal Mass Ejection Speed and Mass

Extreme ultraviolet (EUV) coronal dimmings are often observed in response to solar eruptive events. These phenomena can be generated via several different physical processes. For space weather, the most important of these is the temporary void left behind by a coronal mass ejection (CME). Massive, fast CMEs tend to leave behind a darker void that also usually corresponds to minimum irradiance for the cooler coronal emissions. If the dimming is associated with a solar are, as is often the case, the are component of the irradiance light curve in the cooler coronal emission can be isolated and removed using simultaneous measurements of warmer coronal lines. We apply this technique to 37dimming events identified during two separate two-week periods in 2011, plus an event on 2010 August 7 analyzed in a previous paper, to parameterize dimming in terms of depth and slope. We provide statistics on which combination of wavelengths worked best for the flare-removal method, describe the fitting methods applied to the dimming light curves, and compare the dimming parameters with corresponding CME parameters of mass and speed. The best linear relationships found are nu(sub CME) [km/s] approx. equals 2.36 x 10 6 [km/%] x s(sub dim) [%/s] m(sub CME) [g] approx. equals 2.59 x 10(exp.15 [g/%] x the square root of d(sub dim) [%].These relationships could be used for space weather operations of estimating CME mass and speed using near-real-time irradiance dimming measurements.

Mason, James Paul↗
DOE OSTI2025

Mapping the gas density with the kinematic Sunyaev-Zel’dovich and patchy screening effects: A self-consistent comparison

The secondary anisotropies of the cosmic microwave background (CMB) provide a wealth of astrophysical and cosmological information. Pairing measurements of the CMB temperature map obtained by DR5 of the Atacama Cosmology Telescope (ACT) with the imaging survey conducted by the Dark Energy Spectroscopic Instrument for the purposes of target selection, DECaLS DR9, we investigate two effects that are sensitive to the gas density 𝜏: kinematic Sunyaev-Zel’dovich (kSZ) and patchy screening or anisotropic screening (resulting from the Thomson scattering of CMB photons away from the line-of-sight by free electrons). In particular, we measure the stacked profiles of the gas density around luminous red galaxies (LRGs) at a mean redshift of 𝑧 ≈ 0.7. We detect the kSZ signal at 7.2⁢𝜎, and we find a signal at ∼ 4.1⁢𝜎 for the patchy screening estimator, which is in excess relative to the kSZ signal. We attribute this excess to contamination from CMB lensing. Here, we demonstrate the effect of lensing using 𝑁-body simulations, and we show that the screening signal is dominated by it. Accounting for lensing, our measurement places a 95% upper bound on the optical depth of the Extended DESI LRG sample of 𝜏 < 2.5 10 −4 for a mean value of the sample of 𝜏 ≈ 1.6 10 −4 . Furthermore, via hydro simulations, we show that the underlying optical depth signal measured by both effects (after removing the CMB lensing contribution) is in perfect agreement when adopting either a compensated aperture photometry (CAP) filter or a high-pass filter. Consistent with previous measurements, we see evidence for excess baryonic feedback around DESI LRGs in the patchy screening measurement. In the future, when both effects can be measured with high signal-to-noise, one can measure the amplitude ratio between them, which is proportional to the root-mean-square velocity of the host halo sample, and even place constraints on velocity-sensitive models such as modified gravity and phantom dark energy.

Astrophysical & cosmological simulations↗
DOE OSTI2020

Using respiration quotients to track changing sources of soil respiration seasonally and with experimental warming

Developing a more mechanistic understanding of soil respiration is hampered by the difficulty in determining the contribution of different organic substrates to respiration and in disentangling autotrophic-versus-heterotrophic and aerobic-versus-anaerobic processes. Here, we use a relatively novel tool for better understanding soil respiration: the apparent respiration quotient (ARQ). The ARQ is the amount of CO 2 produced in the soil divided by the amount of O 2 consumed, and it changes according to which organic substrates are being consumed and whether oxygen is being used as an electron acceptor. We investigated how the ARQ of soil gas varied seasonally, by soil depth, and by in situ experimental warming ( +4 ° C) in a coniferous-forest whole-soil-profile warming experiment over 2 years. We then compared the patterns in ARQ to those of soil δ 13 CO 2 . Our measurements showed strong seasonal variations in ARQ, from ≈0.9 during the late spring and summer to ≈0.7 during the winter. This pattern likely reflected a shift from respiration being fueled by oxidized substrates like sugars and organic acids derived from root and root respiration during the growing season to more reduced substrates such as lipids and proteins derived from microbial necromass during the winter. This interpretation was supported by δ 13 CO 2 values, which were lower, like lipids, in the winter and higher, like sugars, in the summer. Furthermore, experimental warming significantly changed how both ARQ and δ 13 CO 2 responded to soil temperature. Wintertime ARQ and δ 13 CO 2 values were higher in heated than in control plots, probably due to the warming-driven increase in microbial activity that may have utilized oxidized carbon substrates, while growing-season values were lower in heated plots. Experimental warming and phenology change the sources of soil respiration throughout the soil profile. The sensitivity of ARQ to these changes demonstrates its potential as a tool for disentangling the biological sources contributing to soil respiration.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2021

Extensive Characterization of Aerosol Optical Properties and Chemical Component Concentrations: Application of the GRASP/Component Approach to Long-term AERONET Measurements

A recently developed GRASP/Component approach (GRASP: Generalized Retrieval of Atmosphere and Surface Properties) was applied to AERONET (Aeronet Robotic Network) sun photometer measurements in this study. Unlike traditional aerosol component retrieval, this approach allows the inference of some information about aerosol composition directly from measured radiance, rather than indirectly through the inversion of optical parameters, and has been integrated into the GRASP algorithm. The newly developed GRASP/Component approach was applied to 13 AERONET sites for different aerosol types under the assumption of aerosol internal mixing rules to analyze the characteristics of aerosol components and their distribution patterns. The results indicate that the retrievals can characterize well the spatial and temporal variability of the component concentration for different aerosol types. A reasonable agreement between GRASP BC retrievals and MERRA-2 BC products is found for all different aerosol types. In addition, the relationships between aerosol component content and aerosol optical parameters such as aerosol optical depth (AOD), fine-mode fraction (FMF), absorption Ångström exponent (AAE), scattering Ångström exponent (SAE), and single scattering albedo (SSA) are also analyzed for indirect verifying the reliability of the component retrieval. It was demonstrated the GRASP/Component optical retrievals are in good agreement with AERONET standard products [e.g., correlation coefficient (R) of 0.93–1.0 for AOD, fine-mode AOD (AODF), coarse-mode AOD (AODC) and Ångström exponent (AE); R = ~ 0.8 for absorption AOD (AAOD) and SSA; RMSE (root mean square error) < 0.03 for AOD, AODF, AODC, AAOD and SSA]. Thus, it is demonstrated the GRASP/Component approach can provide aerosol optical products with comparable accuracy as the AERONET standard products from the ground-based sun photometer measurements as well as some additional important inside on aerosol composition.

Xindan Zhang↗