Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “separated representations”

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

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

At least 415 records · Page 23

Assessing the Atmospheric Response to Subgrid Surface Heterogeneity in the Single-Column Community Earth System Model, Version 2 (CESM2)

Land-atmosphere interactions are central to the evolution of the atmospheric boundary layer and the subsequent formation of clouds and precipitation. Existing global climate models represent these connections with bulk approximations on coarse spatial scales, but observations suggest that small-scale variations in surface characteristics and co-located turbulent and momentum fluxes can significantly impact the atmosphere. Recent model development efforts have attempted to capture this phenomenon by coupling existing representations of subgrid-scale (SGS) heterogeneity between land and atmosphere models. Such approaches are in their infancy and it is not yet clear if they can produce a realistic atmospheric response to surface heterogeneity. Here, we implement a parameterization to capture the effects of SGS heterogeneity in the Community Earth System Model (CESM2), and compare single-column simulations against high-resolution Weather Research and Forecasting (WRF) large-eddy simulations (LESs), which we use as a proxy for observations. The CESM2 experiments increase the temperature and humidity variances in the lowest atmospheric levels, but the response is weaker than in WRF-LES. In part, this is attributed to an underestimate of surface heterogeneity in the land model due to a lack of SGS meteorology, a separation between deep and shallow convection schemes in the atmosphere, and a lack of explicitly represented mesoscale secondary circulations. These results highlight the complex processes involved in capturing the effects of SGS heterogeneity and suggest the need for parameterizations that communicate their influence not only at the surface but also vertically.

54 ENVIRONMENTAL SCIENCES↗

TDD LoRa and Delta Encoding in Low-Power Networks of Environmental Sensor Arrays for Temperature and Deformation Monitoring: Supporting Data

This dataset was used to assess the performance of a compression algorithm for Soil Deformation and Temperature Monitoring data. The development of this technology is presented in the article named time division duplex (TDD) "TDD LoRa and Delta Encoding in Low-Power Networks of Environmental Sensor Arrays for Temperature and Deformation Monitoring" and published in Springer's Journal on Signal Processing Systems. [Time Division Duplex (TDD); Long Range (LoRa)]Field experiment data includes: 1) soil temperature measurements, 2) snow temperature data, 3) soil temperature and deformation data. In these files, the first column contains the time (UTC), followed by the sensor's battery voltage and temperature and acceleration values (X, Y, Z) in subsequent columns. The measurements were acquired every 15 minutes. Lists included for probe, datalogger, and packaging components. The two temperature data sets were acquired at a site at the East River Watershed, Colorado, while the temperature/deformation data were collected at the Teller Road (mile 27) site near Nome, Alaska. The package includes 8 comma-separated files (.csv) and one .pdf file.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Classification of journal surfaces using surface topography parameters and software methods to compensate for stylus geometry

Measurements made with a stylus surface tracer which provides a digitized representation of a surface profile are discussed. Parameters are defined to characterize the height (e.g., RMS roughness, skewness, and kurtosis) and length (e.g., autocorrelation) of the surface topography. These are applied to the characterization of crank shaft journals which were manufactured by different grinding and lopping procedures known to give significant differences in crank shaft bearing life. It was found that three parameters (RMS roughness, skewness, and kurtosis) are necessary to adequately distinguish the character of these surfaces. Every surface specimen has a set of values for these three parameters. They can be regarded as a set coordinate in a space constituted by three characteristics axes. The various journal surfaces can be classified along with the determination of a proper wavelength cutoff (0.25 mm) by using a method of separated subspace. The finite radius of the stylus used for profile tracing gives an inherent measurement error as it passes over the fine structure of the surface. A mathematical model is derived to compensate for this error.

Li, C. J.↗

Nonadditive Interactions Unlock Small-Particle Mobility in Binary Colloidal Monolayers

We examine the organization and dynamics of binary colloidal monolayers composed of micron-scale silica particles interspersed with smaller-diameter silica particles that serve as minority component impurities. These binary monolayers are prepared at the surface of ionic liquid droplets over a range of size ratios (σ = 0.16–0.66) and are studied with low-dose minimally perturbative scanning electron microscopy (SEM). The high resolution of SEM imaging provides direct tracking of all particle coordinates over time, enabling a complete description of the microscopic state. In these bidisperse size mixtures, particle interactions are nonadditive because interfacial pinning to the droplet surface causes the equators of differently sized particles to lie in separate planes. By varying the size ratio, we control the extent of nonadditivity in order to achieve phase behavior inaccessible to additive 2D systems. Across the range of size ratios, we tune the system from a mobile small-particle phase (σ < 0.24) to an interstitial solid (0.24 < σ < 0.33) and furthermore to a disordered glass (σ > 0.33). These distinct phase regimes are classified through measurements of hexagonal ordering of the large-particle host lattice and the lattice’s capacity for small-particle transport. Altogether, we explain these structural and dynamic trends by considering the combined influence of interparticle interactions and the colloidal packing geometry. Our measurements are reproduced in molecular dynamics simulations of 2D nonadditive disks, suggesting an efficient method for describing confined systems with reduced dimensionality representations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An Intercomparison of 2-D Models Within a Common Framework

A model intercomparison among the Atmospheric and Environmental Research (AER) 2-D model, the Goddard Space Flight Center (GSFC) 2-D model, and the Lawrence Livermore National Laboratory 2-D model allows us to separate differences due to model transport from those due to the model's chemical formulation. This is accomplished by constructing two hybrid models incorporating the transport parameters of the GSFC and LLNL models within the AER model framework. By comparing the results from the native models (AER and e.g. GSFC) with those from the hybrid model (e.g. AER chemistry with GSFC transport), differences due to chemistry and transport can be identified. For the analysis, we examined an inert tracer whose emission pattern is based on emission from a High Speed Civil Transport (HSCT) fleet; distributions of trace species in the 2015 atmosphere; and the response of stratospheric ozone to an HSCT fleet. Differences in NO(y) in the upper stratosphere are found between models with identical transport, implying different model representations of atmospheric chemical processes. The response of O3 concentration to HSCT aircraft emissions differs in the models from both transport-dominated differences in the HSCT-induced perturbations of H2O and NO(y) as well as from differences in the model represent at ions of O3 chemical processes. The model formulations of cold polar processes are found to be the most significant factor in creating large differences in the calculated ozone perturbations

Weisenstein, Debra K.↗

An Improved Energy Deposition Model in MPACT and Explicit Heat Generation Coupling with CTF

The default energy deposition model in the CASL neutronics code MPACT assumes all fission energy is deposited locally in fuel rods. Furthermore, equilibrium delayed energy release is assumed for both steady-state and transient calculations. These approximations limit the accurate representation of the heat generation distribution in space and its variations over time, which are essential for power distribution and thermal-hydraulic coupling calculations. In this paper, an improved energy deposition model is presented in both the spatial and time domains. Spatially, the energy deposition through fission, neutron capture, and slowing-down reactions are explicitly modeled to account for the heat generation from all regions of a reactor core, and a gamma smearing scheme is developed that utilizes the gamma sources from neutron fission and capture. In the time domain, the delayed energy release is modeled by solving an additional equation of delayed heat emitters, similar to the equation of delayed neutron precursors.To allow the explicit heat generation coupling, the interfaces between MPACT and CTF were updated to transfer separate heat sources for different material regions (fuel, clad, moderator and guide tube). The results show that the distributions of the energy deposition between MPACT and MCNP agree very well for various 2-D assembly and quarter-core problems without TH feedback. The MPACT/CTF coupled calculation for the hot full power quarter-core case exhibited a reduced peak pin power by 2.3% and a reduced peak fuel centerline temperature by 17 K when using the explicit energy deposition and heat transfer. The new model also shows a maximum 100 pcm keff effect on assembly depletion problems and an increased overall energy release by 7% in a PWR reactivity-initiated accident (RIA) problem.

Liu, Yuxuan↗

Genome-wide association study identifies four pan-ancestry loci for suicidal ideation in the Million Veteran Program

Suicidal ideation (SI) often precedes and predicts suicide attempt and death, is the most common suicidal phenotype and is over-represented in veterans. The genetic architecture of SI in the absence of suicide attempt (SA) is unknown, yet believed to have distinct and overlapping risk with other suicidal behaviors. We performed the first GWAS of SI without SA in the Million Veteran Program (MVP), identifying 99,814 SI cases from electronic health records without a history of SA or suicide death (SD) and 512,567 controls without SI, SA or SD. GWAS was performed separately in the four largest ancestry groups, controlling for sex, age and genetic substructure. Ancestry-specific results were combined via meta-analysis to identify pan-ancestry loci. Four genome-wide significant (GWS) loci were identified in the pan-ancestry meta-analysis with loci on chromosomes 6 and 9 associated with suicide attempt in an independent sample. Pan-ancestry gene-based analysis identified GWS associations with DRD2, DCC, FBXL19, BCL7C, CTF1, ANNK1, and EXD3. Gene-set analysis implicated synaptic and startle response pathways (q’s<0.05). European ancestry (EA) analysis identified GWS loci on chromosomes 6 and 9, as well as GWS gene associations in EXD3, DRD2, and DCC. No other ancestry-specific GWS results were identified, underscoring the need to increase representation of diverse individuals. The genetic correlation of SI and SA within MVP was high (r G = 0.87; p = 1.09e-50), as well as with post-traumatic stress disorder (PTSD; r G = 0.78; p = 1.98e-95) and major depressive disorder (MDD; r G = 0.78; p = 8.33e-83). Conditional analysis on PTSD and MDD attenuated most pan-ancestry and EA GWS signals for SI without SA to nominal significance, with the exception of EXD3 which remained GWS. Our novel findings support a polygenic and complex architecture for SI without SA which is largely shared with SA and overlaps with psychiatric conditions frequently comorbid with suicidal behaviors.

59 BASIC BIOLOGICAL SCIENCES↗

Refining Planetary Boundary Layer Height Retrievals From Micropulse‐Lidar at Multiple ARM Sites Around the World

Abstract Knowledge of the planetary boundary layer height (PBLH) is crucial for various applications in atmospheric and environmental sciences. Lidar measurements are frequently used to monitor the evolution of the PBLH, providing more frequent observations than traditional radiosonde‐based methods. However, lidar‐derived PBLH estimates have substantial uncertainties, contingent upon the retrieval algorithm used. In addressing this, we applied the Different Thermo‐Dynamic Stabilities (DTDS) algorithm to establish a PBLH data set at five separate Department of Energy's Atmospheric Radiation Measurement sites across the globe. Both the PBLH methodology and the products are subject to rigorous assessments in terms of their uncertainties and constraints, juxtaposing them with other products. The DTDS‐derived product consistently aligns with radiosonde PBLH estimates, with correlation coefficients exceeding 0.77 across all sites. This study delves into a detailed examination of the strengths and limitations of PBLH data sets with respect to both radiosonde‐derived and other lidar‐based estimates of the PBLH by exploring their respective errors and uncertainties. It is found that varying techniques and definitions can lead to diverse PBLH retrievals due to the inherent intricacy and variability of the boundary layer. Our DTDS‐derived PBLH data set outperforms existing products derived from ceilometer data, offering a more precise representation of the PBLH. This extensive data set paves the way for advanced studies and an improved understanding of boundary‐layer dynamics, with valuable applications in weather forecasting, climate modeling, and environmental studies.

54 ENVIRONMENTAL SCIENCES↗

African Easterly Jet: Structure and Maintenance

This article investigates the African Easterly Jet (AEJ), its structure and the forcings contributing to its maintenance, critically revisiting previous work which attributed the maintenance of the jet to soil moisture gradients over tropical Africa. A state-of-the-art global model in a high-end computer framework is used to produce a 3-member 73-year ensemble run forced by observed SST to represent the Control run. The AEJ as produced by the Control is compared with the representation of the AEJ in the European Center for Medium Range Forecast Reanalyses (ERA-40) and other observational data sets and found very realistic. Five Experiments are then performed, each represented by sets of 3-member 22 year long (1980-2001) ensemble runs. The goal of the Experiments is to investigate the role of meridional soil moisture gradients, different land surface properties and orography. Unlike previous studies, which have suppressed soil moisture gradients within a highly idealized framework (i.e., the so-called bucket model), terrestrial evaporation control is here achieved with a highly sophisticated landsurface treatment and with an extensively tested and complex methodology. The results show that the AEJ is suppressed by a combination of absence of meridional evaporation gradients over Africa and constant vegetation, even if the individual forcings taken separately do not lead to the AEJ disappearance, but only its modification. Moreover, the suppression of orography also leads to a different circulation in which there is no AEJ. This work suggests that it is not just soil moisture gradients, but a unique combination of geographical features present only in northern tropical Africa, which causes and maintains the jet.

Wu, Man-Li C.↗

Disentangling the Regional Climate Impacts of Competing Vegetation Responses to Elevated Atmospheric CO2

Biophysical vegetation responses to elevated atmospheric carbon dioxide (CO2) affect regional hydroclimate through two competing mechanisms. Higher CO2 increases leaf area (LAI), thereby increasing transpiration and water losses. Simultaneously, elevated CO2 reduces stomatal conductance and transpiration, thereby increasing rootzone soil moisture. Which mechanism dominates in the future is highly uncertain, partly because these two processes are difficult to explicitly separate within dynamic vegetation models. We address this challenge by using the GISS ModelE global climate model to conduct a novel set of idealized 2×CO2 sensitivity experiments to: evaluate the total vegetation biophysical contribution to regional climate change under high CO2; and quantify the separate contributions of enhanced LAI and reduced stomatal conductance to regional hydroclimate responses. We find that increased LAI exacerbates soil moisture deficits across the sub‐tropics and more water‐limited regions, but also attenuates warming by ∼0.5–1°C in the US Southwest, Central Asia, Southeast Asia, and northern South America. Reduced stomatal conductance effects contribute ∼1°C of summertime warming. For some regions, enhanced LAI and reduced stomatal conductance produce nonlinear and either competing or mutually amplifying hydroclimate responses. In northeastern Australia, these effects combine to exacerbate radiation‐forced warming and contribute to year‐round water limitation. Conversely, at higher latitudes these combined effects result in less warming than would otherwise be predicted due to nonlinear responses. These results highlight substantial regional variation in CO2‐driven vegetation responses and the importance of improving model representations of these processes to better quantify regional hydroclimate impacts.

Regional climate impacts↗

Automated characterization of spatial and dynamical heterogeneity in supercooled liquids via implementation of machine learning

Abstract A computational approach by an implementation of the principle component analysis (PCA) with K -means and Gaussian mixture (GM) clustering methods from machine learning algorithms to identify structural and dynamical heterogeneities of supercooled liquids is developed. In this method, a collection of the average weighted coordination numbers ( W C N s ‾ ) of particles calculated from particles’ positions are used as an order parameter to build a low-dimensional representation of feature (structural) space for K -means clustering to sort the particles in the system into few meso-states using PCA. Nano-domains or aggregated clusters are also formed in configurational (real) space from a direct mapping using associated meso-states’ particle identities with some misclassified interfacial particles. These classification uncertainties can be improved by a co-learning strategy which utilizes the probabilistic GM clustering and the information transfer between the structural space and configurational space iteratively until convergence. A final classification of meso-states in structural space and domains in configurational space are stable over long times and measured to have dynamical heterogeneities. Armed with such a classification protocol, various studies over the thermodynamic and dynamical properties of these domains indicate that the observed heterogeneity is the result of liquid–liquid phase separation after quenching to a supercooled state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Three-Dimensional Analysis of Deep Space Network Antenna Coverage

There is a need to understand NASA s Deep Space Network (DSN) coverage gaps and any limitations to provide redundant communication coverage for future deep space missions, especially for manned missions to Moon and Mars. The DSN antennas are required to provide continuous communication coverage for deep space flights, interplanetary missions, and deep space scientific observations. The DSN consists of ground antennas located at three sites: Goldstone in USA, Canberra in Australia, and Madrid in Spain. These locations are not separated by the exactly 120 degrees and some DSN antennas are located in the bowl-shaped mountainous terrain to shield against radiofrequency interference resulting in a coverage gap in the southern hemisphere for the current DSN architecture. To analyze the extent of this gap and other coverage limitations, simulations of the DSN architecture were performed. In addition to the physical properties of the DSN assets, the simulation incorporated communication forward link calculations and azimuth/elevation masks that constrain the effects of terrain for each DSN antenna. Analysis of the simulation data was performed to create coverage profiles with the receiver settings at a deep space altitudes ranging from 2 million to 10 million km and a spherical grid resolution of 0.25 degrees with respect to longitude and latitude. With the results of these simulations, two- and three-dimensional representations of the area without communication coverage and area with coverage were developed, showing the size and shape of the communication coverage gap projected in space. Also, the significance of this communication coverage gap is analyzed from the simulation data.

Kegege, Obadiah↗

Freshwater Input and Vertical Mixing in the Canada Basin’s Seasonal Halocline: 1975 versus 2006–12

Abstract The Arctic seasonal halocline impacts the exchange of heat, energy, and nutrients between the surface and the deeper ocean, and it is changing in response to Arctic sea ice melt over the past several decades. Here, we assess seasonal halocline formation in 1975 and 2006–12 by comparing daily, May–September, salinity profiles collected in the Canada Basin under sea ice. We evaluate differences between the two time periods using a one-dimensional (1D) bulk model to quantify differences in freshwater input and vertical mixing. The 1D metrics indicate that two separate factors contribute similarly to stronger stratification in 2006–12 relative to 1975: 1) larger surface freshwater input and 2) less vertical mixing of that freshwater. The larger freshwater input is mainly important in August–September, consistent with a longer melt season in recent years. The reduced vertical mixing is mainly important from June until mid-August, when similar levels of freshwater input in 1975 and 2006–12 are mixed over a different depth range, resulting in different stratification. These results imply that decadal changes to ice–ocean dynamics, in addition to freshwater input, significantly contribute to the stronger seasonal stratification in 2006–12 relative to 1975. These findings highlight the need for near-surface process studies to elucidate the impact of lateral processes and ice–ocean momentum exchange on vertical mixing. Moreover, the results may provide insight for improving the representation of decadal changes to Arctic upper-ocean stratification in climate models that do not capture decadal changes to vertical mixing.

54 ENVIRONMENTAL SCIENCES↗

Transcriptomics-based Machine Learning (ML) Analysis Predicts Space-Exposed Murine Livers

Limited sample sizes, high data dimensionality, and sensitivity to technical and biological variability of next generation sequencing (NGS), typically limits machine learning (ML) approaches in spaceflight studies that include radiation effects. However, pooling smaller studies while addressing intra- and inter-study variabilities allows for ML predictive modeling. Here, integration methods were applied to whole transcriptome shotgun sequencing (RNA-seq) data from six mouse liver GeneLab datasets (GLDS) (n ranging from 6 to 39 samples) from with a total of 81 spaceflight and ground-control samples to determine top features (i.e. genes) relevant to spaceflight including the effect of radiation exposure. RNASeq counts were normalized for each study, then merged and scaled across all datasets. Data dimensionality was reduced using a minimum redundancy maximum relevance (MRMR) methodology. Redundancy and relevance were computed using the Pearson correlation and F-statistic, respectively. The top 100 MRMR features were used to predict spaceflight vs. ground-control samples using Random Forest (RF), Support Vector Machine (SVM), and Linear Discriminant Analysis (LDA) classifiers with 5-fold cross validation (CV). Principal component analysis (PCA) on the complete feature set versus the MRMR features shows separation between spaceflight samples and ground controls (Figure 1A). The ML-based gene sets were compared against differential gene expression results obtained with DESeq2 from individual GLDS. Using all features or randomly sampled subsets at matching set sizes with MRMR, a maximum classifier accuracy of 69% on the test set over 5 folds. For all classifiers, CV training using at least the top 30 MRMR genes show minimum 89% accuracy and 0.95 AUC value on the test set over 5 folds (Figure 1B). Baseline set analysis on differentially expressed genes (DEGs) identified using padj ≤ 0.05 show 295 DEGs that overlap at least two studies and 13 DEGs that overlap three studies (Figure 1C). Set analysis between the top 100 MRMR features and the DEGs showed 47 genes that overlap at least one study and 24 genes that overlap two studies. Over-representation analysis showed overlapping biological processes related to fatty acid and lipid metabolism which may indicate these processes in the response to spaceflight stressors. MRMR feature selection for the selected ML methods improve performance relative to a classifier built on all features or randomly sampled subsets. Permutation feature importance within the decorrelated MRMR features showed concordance in feature ranking between ML methods. A challenge of applying ML methods across heterogeneous NGS data is accounting for signal:noise. Here, signal validation across studies was shown by intersecting sets between top MRMR genes and DEGs from DESeq2 analysis. Non-intersecting sets introduce opportunity to explore genes relevant to differentiating space flight exposed groups and implementing ML methods across existing NGS datasets may overcome sample size limitations.

Machine Learning↗

Transcriptomics-based Machine Learning Analysis Predicts Space-Exposed Murine Livers

Limited sample sizes, high data dimensionality, and sensitivity to technical and biological variability of next generation sequencing (NGS), typically limits machine learning (ML) approaches in spaceflight studies that include radiation effects. However, pooling smaller studies while addressing intra- and inter-study variabilities allows for ML predictive modeling. Here, integration methods were applied to whole transcriptome shotgun sequencing (RNA-seq) data from six mouse liver GeneLab datasets (GLDS) (n ranging from 6 to 39 samples) from with a total of 81 spaceflight and ground-control samples to determine top features (i.e. genes) relevant to spaceflight including the effect of radiation exposure. RNASeq counts were normalized for each study, then merged and scaled across all datasets. Data dimensionality was reduced using a minimum redundancy maximum relevance (MRMR) methodology. Redundancy and relevance were computed using the Pearson correlation and F-statistic, respectively. The top 100 MRMR features were used to predict spaceflight vs. ground-control samples using Random Forest (RF), Support Vector Machine (SVM), and Linear Discriminant Analysis (LDA) classifiers with 5-fold cross validation (CV). Principal component analysis (PCA) on the complete feature set versus the MRMR features shows separation between spaceflight samples and ground controls (Figure 1A). The ML-based gene sets were compared against differential gene expression results obtained with DESeq2 from individual GLDS. Using all features or randomly sampled subsets at matching set sizes with MRMR, a maximum classifier accuracy of 69% was shown on the test set over 5 folds. For all classifiers, CV training using at least the top 30 MRMR genes show minimum 89% accuracy and 0.95 AUC value on the test set over 5 folds (Figure 1B). Baseline set analysis on differentially expressed genes (DEGs) identified using padj ≤ 0.05 show 295 DEGs that overlap at least two studies and 13 DEGs that overlap three studies (Figure 1C). Set analysis between the top 100 MRMR features and the DEGs showed 47 genes that overlap at least one study and 24 genes that overlap two studies. Over-representation analysis showed overlapping biological processes related to fatty acid and lipid metabolism which may indicate these processes in the response to spaceflight stressors. MRMR feature selection for the selected ML methods improve performance relative to a classifier built on all features or randomly sampled subsets. Permutation feature importance within the decorrelated MRMR features showed concordance in feature ranking between ML methods. A challenge of applying ML methods across heterogeneous NGS data is accounting for signal:noise. Here, signal validation across studies was shown by intersecting sets between top MRMR genes and DEGs from DESeq2 analysis. Non-intersecting sets introduce opportunity to explore genes relevant to differentiating space flight exposed groups and implementing ML methods across existing NGS datasets may overcome sample size limitations.

Machine Learning↗

Transcriptomics-based Machine Learning Analysis Predicts Space-Exposed Murine Livers

Limited sample sizes, high data dimensionality, and sensitivity to technical and biological variability of next generation sequencing (NGS), typically limits machine learning (ML) approaches in spaceflight studies that include radiation effects. However, pooling smaller studies while addressing intra- and inter-study variabilities allows for ML predictive modeling. Here, integration methods were applied to whole transcriptome shotgun sequencing (RNA-seq) data from six mouse liver GeneLab datasets (GLDS) (n ranging from 6 to 39 samples) from with a total of 81 spaceflight and ground-control samples to determine top features (i.e. genes) relevant to spaceflight including the effect of radiation exposure. RNASeq counts were normalized for each study, then merged and scaled across all datasets. Data dimensionality was reduced using a minimum redundancy maximum relevance (MRMR) methodology. Redundancy and relevance were computed using the Pearson correlation and F-statistic, respectively. The top 100 MRMR features were used to predict spaceflight vs. ground-control samples using Random Forest (RF), Support Vector Machine (SVM), and Linear Discriminant Analysis (LDA) classifiers with 5-fold cross validation (CV). Principal component analysis (PCA) on the complete feature set versus the MRMR features shows separation between spaceflight samples and ground controls (Figure 1A). The ML-based gene sets were compared against differential gene expression results obtained with DESeq2 from individual GLDS. Using all features or randomly sampled subsets at matching set sizes with MRMR, a maximum classifier accuracy of 69% was shown on the test set over 5 folds. For all classifiers, CV training using at least the top 30 MRMR genes show minimum 89% accuracy and 0.95 AUC value on the test set over 5 folds (Figure 1B). Baseline set analysis on differentially expressed genes (DEGs) identified using padj ≤ 0.05 show 295 DEGs that overlap at least two studies and 13 DEGs that overlap three studies (Figure 1C). Set analysis between the top 100 MRMR features and the DEGs showed 47 genes that overlap at least one study and 24 genes that overlap two studies. Over-representation analysis showed overlapping biological processes related to fatty acid and lipid metabolism which may indicate these processes in the response to spaceflight stressors. MRMR feature selection for the selected ML methods improve performance relative to a classifier built on all features or randomly sampled subsets. Permutation feature importance within the decorrelated MRMR features showed concordance in feature ranking between ML methods. A challenge of applying ML methods across heterogeneous NGS data is accounting for signal:noise. Here, signal validation across studies was shown by intersecting sets between top MRMR genes and DEGs from DESeq2 analysis. Non-intersecting sets introduce opportunity to explore genes relevant to differentiating space flight exposed groups and implementing ML methods across existing NGS datasets may overcome sample size limitations.

Machine Learning↗

Limitations of Separate Cloud and Rain Categories in Parameterizing Collision-Coalescence for Bulk Microphysics Schemes

Warm rain collision-coalescence has been persistently difficult to parameterize in bulk microphysics schemes. We use a flexible bulk microphysics scheme with bin scheme process parameterizations, called AMP, to investigate reasons for the difficulty. AMP is configured in a variety of ways to mimic bulk schemes and is compared to simulations with the bin scheme upon which AMP is built. We find that an important limitation in traditional bulk schemes is the use of separate cloud and rain categories. When the drop size distribution is instead represented by a continuous distribution, the simulation of cloud-to-rain conversion is substantially improved. We also find large sensitivity to the threshold size to distinguish cloud and rain in traditional schemes; substantial improvement is found by decreasing the threshold from 40 to 25 μm. Neither the use of an assumed functional form for the size distribution nor the choice of predicted distribution moments has a large impact on the ability of AMP to simulate rain production. When predicting four total moments of the liquid drop size distribution, either with a traditional two-category, two-moment scheme with a reduced size threshold, or a four-moment single-category scheme, errors in the evolution of mass and the cloud size distribution are similar, but the single-category scheme has a substantially better representation of the rain size distribution. Optimal moment combinations for the single-category approach are investigated and appear to be linked more to the information content they provide for constraining the size distributions than to their correlation with collision-coalescence rates.

54 ENVIRONMENTAL SCIENCES↗

Impact of Asian Aerosols on Precipitation Over California: An Observational and Model Based Approach

Dust and pollution emissions from Asia are often transported across the Pacific Ocean to over the western United States. Therefore, it is essential to fully understand the impact of these aerosols on clouds and precipitation forming over the eastern Pacific and western United States, especially during atmospheric river events that account for up to half of California's annual precipitation and can lead to widespread flooding. In order for numerical modeling simulations to accurately represent the present and future regional climate of the western United States, we must account for the aerosol-cloud-precipitation interactions associated with Asian dust and pollution aerosols. Therefore, we have constructed a detailed study utilizing multi-sensor satellite observations, NOAA-led field campaign measurements, and targeted numerical modeling studies where Asian aerosols interacted with cloud and precipitation processes over the western United States. In particular, we utilize aerosol optical depth retrievals from the NASA Moderate Resolution Imaging Spectroradiometer (MODIS), NOAA Geostationary Operational Environmental Satellite (GOES-11), and Japan Meteorological Agency (JMA) Multi-functional Transport Satellite (MTSAT) to effectively detect and monitor the trans-Pacific transport of Asian dust and pollution. The aerosol optical depth (AOD) retrievals are used in assimilating the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem) in order to provide the model with an accurate representation of the aerosol spatial distribution across the Pacific. We conduct WRF-Chem model simulations of several cold-season atmospheric river events that interacted with Asian aerosols and brought significant precipitation over California during February-March 2011 when the NOAA CalWater field campaign was ongoing. The CalWater field campaign consisted of aircraft and surface measurements of aerosol and precipitation processes that help extensively validate our WRF-Chem model simulations. After validating the capability of the WRF-Chem in realistically simulating the aerosol-cloud precipitation interactions, we conduct sensitivity studies where the AOD is doubled to diagnose whether an increasing concentration of Asian aerosols over the western United States will lead to further impacts on the cloud and precipitation processes over California. We also perform sensitivity studies where the aerosols will be partitioned into dust-only and pollution-only in order to separate the impacts of the differing Asian aerosol species. The results of our WRF-Chem model simulations aim to show that the trans-Pacific transport of Asian aerosols influence the precipitation associated with atmospheric river events that can ultimately impact the regional climate of the western United States. 1 University

Naeger, Aaron R.↗