Search NASA⌕ Search

SEARCH · Search NASA

Results for “biophysical simulations”

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 91 records · Page 5

cabbi-bio/BECT-ABM

The emergence of the cellulosic bioeconomy requires not only adequate technological, economic, and policy advances, but also effective communication and coordination among the multiple stakeholders in the bioeconomy community. Aiming to facilitate stakeholder communication and collaboration for cellulosic bioeconomy growth, an agent-based model is developed as a community communication tool in this study. The tool simulates the behaviors and interactions of key stakeholders, as well as the responses and interactions between the stakeholders and biophysical, socioeconomic, and environmental conditions. The agent behavior rules are developed and validated based on behavior and economic theories and information attained from a stakeholder focus group study and a land use survey.

Yang, Pan↗

Quantum Hardware-Enabled Molecular Dynamics via Transfer Learning

The ability to perform ab initio molecular dynamics simulations using potential energy surfaces provided by quantum computers would open the door to virtually exact dynamics for a variety of chemical and biochemical systems, with impacts on catalysis and biophysics. Nonetheless, performing molecular dynamics on surfaces produced by quantum hardware has been hampered by the noisy energies typically produced by quantum computers and challenges associated with computing gradients and scaling to large systems interest. A recent set of advances in machine learning, known as transfer learning, provides a new path forward for molecular dynamics simulations on quantum hardware. Transfer learning offers a workaround, where one first trains models on larger, less accurate classical datasets and then refines them on smaller, more accurate quantum datasets. We explore this approach by training machine learning models to predict a molecule's potential energy based on its geometric structure using Behler-Parrinello neural networks. When successfully trained, the model enables energy gradient predictions necessary for dynamic simulations. To reduce the quantum resources needed, the model is initially trained with data derived from classical density functional theory and subsequently refined with a smaller dataset obtained from a variational quantum eigensolver optimization of the unitary coupled cluster ansatz. We show that this approach significantly reduces the size of the needed quantum training dataset while capturing the high accuracies needed within quantum chemistry simulations. The success of this two-step training method opens more opportunities to apply machine learning models on quantum data, a significant stride towards efficient quantum-classical hybrid computational models.

quantum computing↗

Hyperspectral Radiative Transfer Modeling to Explore the Combined Retrieval of Biophysical Parameters and Canopy Fluorescence from FLEX - Sentinel-3 Tandem Mission Multi-Sensor Data

The FLuorescence EXplorer (FLEX) satellite mission, selected as ESA's 8th Earth Explorer, has been designed forthe measurement of sun-induced fluorescence (F) spectra emitted by plants. This will be accomplished through amulti-sensor approach by placing it in a common orbit in tandem with the Sentinel-3 (S3) mission, which willhave two optical sensors on board, OLCI (Ocean and Land Colour Instrument) and SLSTR (Sea and Land SurfaceTemperature Radiometer) to complement FLEX. These S3 instruments will be used in combination with theimaging spectrometers on board FLEX to provide data useful for atmospheric correction of FLEX data. However,a fully synergetic approach, i.e. by exploiting the spectral and directional information from all tandem missioninstruments together, is an attractive alternative which is explored in this paper. By employing all combined topof-atmosphere (TOA) spectral radiance data, one can (i) characterize the relevant optical properties of the atmosphere,(ii) retrieve biophysical canopy properties including the associated reflectance anisotropy, and (iii)retrieve a more accurate and consistent canopy F.Regarding retrieval methods, Fraunhofer Line Depth (FLD) and Spectral Fitting (SF) are well-known techniquesapplied to hyperspectral data. Both methods depend on a high spectral resolution and assume aLambertian (isotropic) canopy reflectance. However, most vegetation canopies are non-Lambertian. This impliesthat, in particular when ignoring the anisotropic surface reflection, substantial retrieval errors can occur due tothe interaction between atmospheric absorption bands and surface reflectance anisotropy. In this paper, a novelmethod based on spectral radiative transfer (RT) modeling is proposed, in which coupled RT models are used tosimulate TOA radiance spectra. These are then matched with ‘measured' spectra in order to retrieve surfacefluorescence, along with a suite of biophysical parameters, by model inversion through optimization. By applyingcoupled RT models of the soil-leaf-canopy and the surface-atmosphere systems, TOA radiance spectra canbe simulated for all optical sensors of this tandem mission. In this way, complex effects due to surface reflectanceanisotropy and the spectral sampling by the various instruments, which are difficult to compensate for in the endproducts, are properly taken into account by their incorporation in the forward modeling. Next, by model inversionof TOA radiance data via optimization, the most accurate F retrievals can be achieved in a consistentmanner, along with important canopy level biophysical parameters that may help interpret the F spectrum, suchas chlorophyll content and leaf area index (LAI). The potential of this approach has been explored in a numericalexperiment, and the results are presented in this paper. We find that, with the assumed well-characterized andplausible FLEX/S3 instrument performances, the simultaneous retrieval of biophysical canopy parameters and Fspectra would be possible with a remarkable accuracy, provided the correct atmospheric characterization isavailable.

Verhoef, Wouter↗

A 1D–0D–3D coupled model for simulating blood flow and transport processes in breast tissue

In this work, we present mixed dimensional models for simulating blood flow and transport processes in breast tissue and the vascular tree supplying it. These processes are considered, to start from the aortic inlet to the capillaries and tissue of the breast. Large variations in biophysical properties and flow conditions exist in this system necessitating the use of different flow models for different geometries and flow regimes. In total, we consider four different model types. First, a system of 1D nonlinear hyperbolic partial differential equations (PDEs) is considered to simulate blood flow in larger arteries with highly elastic vessel walls. Second, we assign 1D linearized hyperbolic PDEs to model the smaller arteries with stiffer vessel walls. The third model type consists of ODE systems (0D models). It is used to model the arterioles and peripheral circulation. Finally, homogenized 3D porous media models are considered to simulate flow and transport in capillaries and tissue within the breast volume. Sink terms are used to account for the influence of the venous and lymphatic systems. Combining the four model types, we obtain two different 1D–0D–3D coupled models for simulating blood flow and transport processes: The first model results in a fully coupled 1D–0D–3D model covering the complete path from the aorta to the breast combining a generic arterial network with a patient specific breast network and geometry. The second model is a reduced one based on the separation of the generic and patient specific parts. The information from a calibrated fully coupled model is used as inflow condition for the patient specific sub-model allowing a significant computational cost reduction. Several numerical experiments are conducted to calibrate the generic model parameters and to demonstrate realistic flow simulations compared to existing data on blood flow in the human breast and vascular system. Moreover, we use two different breast vasculature and tissue data sets to illustrate the robustness of our reduced sub-model approach.

multiscale blood flow models↗

A metabolic modeling platform for the computation of microbial ecosystems in time and space (COMETS)

Genome-scale stoichiometric modeling of metabolism has become a standard systems biology tool for modeling cellular physiology and growth. Extensions of this approach are emerging as a valuable avenue for predicting, understanding and designing microbial communities. Computation of microbial ecosystems in time and space (COMETS) extends dynamic flux balance analysis to generate simulations of multiple microbial species in molecularly complex and spatially structured environments. Here we describe how to best use and apply the most recent version of COMETS, which incorporates a more accurate biophysical model of microbial biomass expansion upon growth, evolutionary dynamics and extracellular enzyme activity modules. In addition to a command-line option, COMETS includes user-friendly Python and MATLAB interfaces compatible with the well-established COBRA models and methods, as well as comprehensive documentation and tutorials. Overall, this protocol provides a detailed guideline for installing, testing and applying COMETS to different scenarios, generating simulations that take from a few minutes to several days to run, with broad applicability to microbial communities across biomes and scales.

59 BASIC BIOLOGICAL SCIENCES↗

Identifying sequence perturbations to an intrinsically disordered protein that determine its phase-separation behavior

Phase separation of intrinsically disordered proteins (IDPs) commonly underlies the formation of membraneless organelles, which compartmentalize molecules intracellularly in the absence of a lipid membrane. Identifying the protein sequence features responsible for IDP phase separation is critical for understanding physiological roles and pathological consequences of biomolecular condensation, as well as for harnessing phase separation for applications in bioinspired materials design. To expand our knowledge of sequence determinants of IDP phase separation, we characterized variants of the intrinsically disordered RGG domain from LAF-1, a model protein involved in phase separation and a key component of P granules. Based on a predictive coarse-grained IDP model, we identified a region of the RGG domain that has high contact probability and is highly conserved between species; deletion of this region significantly disrupts phase separation in vitro and in vivo. We determined the effects of charge patterning on phase behavior through sequence shuffling. We designed sequences with significantly increased phase separation propensity by shuffling the wild-type sequence, which contains well-mixed charged residues, to increase charge segregation. This result indicates the natural sequence is under negative selection to moderate this mode of interaction. We measured the contributions of tyrosine and arginine residues to phase separation experimentally through mutagenesis studies and computationally through direct interrogation of different modes of interaction using all-atom simulations. Finally, we show that despite these sequence perturbations, the RGG-derived condensates remain liquid-like. Together, these studies advance our fundamental understanding of key biophysical principles and sequence features important to phase separation.

36 MATERIALS SCIENCE↗

NASA Space Radiation Program Integrative Risk Model Toolkit

NASA Space Radiation Program Element scientists have been actively involved in development of an integrative risk models toolkit that includes models for acute radiation risk and organ dose projection (ARRBOD), NASA space radiation cancer risk projection (NSCR), hemocyte dose estimation (HemoDose), GCR event-based risk model code (GERMcode), and relativistic ion tracks (RITRACKS), NASA radiation track image (NASARTI), and the On-Line Tool for the Assessment of Radiation in Space (OLTARIS). This session will introduce the components of the risk toolkit with opportunity for hands on demonstrations. The brief descriptions of each tools are: ARRBOD for Organ dose projection and acute radiation risk calculation from exposure to solar particle event; NSCR for Projection of cancer risk from exposure to space radiation; HemoDose for retrospective dose estimation by using multi-type blood cell counts; GERMcode for basic physical and biophysical properties for an ion beam, and biophysical and radiobiological properties for a beam transport to the target in the NASA Space Radiation Laboratory beam line; RITRACKS for simulation of heavy ion and delta-ray track structure, radiation chemistry, DNA structure and DNA damage at the molecular scale; NASARTI for modeling of the effects of space radiation on human cells and tissue by incorporating a physical model of tracks, cell nucleus, and DNA damage foci with image segmentation for the automated count; and OLTARIS, an integrated tool set utilizing HZETRN (High Charge and Energy Transport) intended to help scientists and engineers study the effects of space radiation on shielding materials, electronics, and biological systems.

Kim, Myung-Hee Y.↗

Solvent-Mediated Control of Nanocellulose Dispersion: An Integrated Computational and Experimental Investigation

Fibrillated cellulose derived from forestry feedstocks represents a renewable and high-strength materials platform for circular bioeconomies. However, its practical implementation is hindered by the irreversible aggregation of nanocellulose architectures, including cellulose nanofibers (CNFs). Solvent-based dispersion offers a simple and practical route to prevent CNF aggregation. Here, in this work, we integrate classical and enhanced sampling molecular dynamics (MD) simulations with experimental suspension rheology and atomic force microscopy (AFM) to elucidate how solvent environments tune CNF–CNF interactions and dispersion stability. CNF–CNF contact free energies computed from MD simulations reveal reduced aggregation in acetone/water, γ-valerolactone (GVL)/water, and tetrahydrofuran (THF)/water and pure acetone compared with pure water, reflecting stronger CNF-solvent relative to inter-CNF interactions. Correspondingly, CNF-solvent suspensions in these solvent systems exhibit stronger inter-fibril network structures and enhanced recovery compared to water, indicating improved CNF-solvent affinity. Liquid cell AFM imaging in acetone–water mixtures and in pure acetone further confirm the presence of well-dispersed CNFs. By combining multiscale computation with targeted experiments, this study establishes a rational framework for solvent design to achieve stable nanocellulose dispersions for high-strength biobased materials and efficient bioenergy conversion.

cellulose↗

Soil Surface Runoff Scheme for Improving Land-Hydrology and Surface Fluxes in Simple SiB (SSiB)

Evapotranspiration on land is hard to measure and difficult to simulate. On the scale of a GCM grid, there is large subgrid-scale variability of orography, soil moisture, and vegetation. Our hope is to be able to tune the biophysical constants of vegetation and soil parameters to get the most realistic space-averaged diurnal cycle of evaporation and its climatology. Field experiments such as First ISLSCP Field Experiment (FIFE), Boreal Ecosystem-Atmosphere Study (BOREAS), and LBA help a great deal in improving our evapotranspiration schemes. However, these improvements have to be matched with, and coupled to, consistent improvement in land-hydrology; otherwise, the runoff problems will intrinsically reflect on the soil moisture and evapotranspiration errors. Indeed, a realistic runoff simulation also ensures a reasonable evapotranspiration simulation provided the precipitation forcing is reliable. We have been working on all of the above problems to improve the simulated hydrologic cycle. Through our participation in the evaluation and intercomparison of land-models under the behest of Global Soil Wetness Project (GSWP), we identified a few problems with Simple SiB (SSIB; Xue et al., 1991) hydrology in regions of significant snowmelt. Sud and Mocko (1999) show that inclusion of a separate snowpack model, with its own energy budget and fluxes with the atmosphere aloft and soil beneath, helps to ameliorate some of the deficiencies of delayed snowmelt and excessive spring season runoff. Thus, much more realistic timing of melt water generation was simulated with the new snowpack model in the subsequent GSWP re-evaluations using 2 years of ISLSCP Initiative I forcing data for 1987 and 1988. However, we noted an overcorrection of the low meltwater infiltration of SSiB. While the improvement in snowmelt timing was found everywhere, the snowmelt infiltration has became excessive in some regions, e.g., Lena river basin. This leads to much reduced runoff in many basins as compared to observations. We believe this is a consequence of neglect of the influence of subgrid-scale variations in orography that affects the production of surface runoff.

Sud, Y. C.↗

Modeling Functional Organic Chemistry in Arctic Rivers: An Idealized Siberian System

Rivers of the Arctic will become ever more important for the global climate, since they carry a majority of continental dissolved organic carbon flux into the rapidly changing polar ocean. Aqueous organics comprise a wide array of functional groups, several of which are likely to impact coastal and open water biophysical properties. Light attenuation, interfacial films, aerosol formation, gas release and momentum exchange can all be cited. We performed Lagrangian kinetic modeling for the evolution of riverine organic chemistry as the molecules in question make their way from the highlands to Arctic outlets. Classes as diverse as the proteins, sugars, lipids, re-condensates, humics, bio-tracers and small volatiles are all included. Our reduced framework constitutes an idealized northward flow driving a major hydrological discharge rate and primarily representing the Russian Lena. Mountainous, high solute and tundra sources are all simulated, and they meet up at several points between soil and delta process reactors. Turnover rates are parameterized beginning with extrapolated coastal values imposed along a limited tributary network, with connections between different terrestrial sub-ecologies. Temporal variation of our total dissolved matter most closely resembles the observations when we focus on the restricted removal and low initial carbon loads, suggesting relatively slow transformation along the water course. Thus, channel combinations and mixing must play a dominant role. Nevertheless, microbial and photochemical losses help determine the final concentrations for most species. Chemical evolution is distinct for the various functionalities, with special contributions from pre- and post-reactivity in soil and delta waters. Several functions are combined linearly to represent the collective chromophoric dissolved matter, characterized here by its absorption. Tributaries carry the signature of lignin phenols to segregate tundra versus taiga sources, and special attention is paid to the early then marine behaviors of low molecular weight volatiles. Heteropolycondensates comprise the largest percentage of reactive carbon in our simulations due to recombination/accumulation, and they tend to be preeminent at the mouth. Outlet concentrations of individual structures such as amino acids and absorbers lie above threshold values for biophysical influence, on the monolayer and light attenuation. The extent of coastal spreading is examined through targeted regional box modeling, relying on salinity and color for calibration. In some cases, plumes reach the scale of peripheral arctic seas, and amplification is expected during upcoming decades. Conclusions are mapped from the Lena to other boreal discharges, and future research questions are outlined regarding the bonding type versus mass release as permafrost degrades. Dynamic aqueous organic coupling is recommended for polar system models, from headwaters to coastal diluent.

54 ENVIRONMENTAL SCIENCES↗

CO 2 fertilization of crops offsets yield losses due to future surface ozone damage and climate change

Tropospheric ozone (O 3 ) is harmful to plant productivity and negatively impacts crop yields. O3 concentrations are projected to decrease globally in the optimistic Representative Concentration Pathway of 2.6 W m –2 (RCP2.6) but increase globally following the high-emission scenario under the RCP8.5, with substantial implications for global food security. The damaging effect of O 3 on future crop yield is affected by CO 2 fertilization and climate change, and their interactions for RCP scenarios have yet to be quantified. In this study, we used the Joint UK Land Environment Simulator modified to include crops (JULES-crop) to quantify the impacts, and relative importance of present-day and future O3, CO 2 concentration and meteorology on crop production at the regional scale until 2100 following RCP2.6 and RCP8.5 scenarios. We focus on eight major crop-producing regions that cover the production of wheat, soybean, maize, and rice. Our results show that CO 2 alone has the largest effect on regional yields, followed by climate and O 3 . However, the CO 2 fertilization effect is offset by the negative impact of tropospheric O 3 in regions with high O 3 concentrations, such as South Asia and China. Simulated crop yields in 2050 were compared with Food and Agriculture Organisation (FAO) statistics to investigate the differences between a socioeconomic and a biophysical process-based approach. Results showed that FAO estimates are closer to our JULES-crop RCP8.5 scenario. This study demonstrates that air pollution could be the biggest threat to future food production and highlights an urgent policy need to mitigate the threat of climate change and O 3 pollution on food security.

54 ENVIRONMENTAL SCIENCES↗

Large-scale experimental technology with remote sensing in land surface hydrology and meteorology

Two field experiments to study atmospheric and land surface processes and their interactions are summarized. The Hydrologic-Atmospheric Pilot Experiment, which tested techniques for measuring evaporation, soil moisture storage, and runoff at scales of about 100 km, was conducted over a 100 X 100 km area in France from mid-1985 to early 1987. The first International Satellite Land Surface Climatology Program field experiment was conducted in 1987 to develop and use relationships between current satellite measurements and hydrologic, climatic, and biophysical variables at the earth's surface and to validate these relationships with ground truth. This experiment also validated surface parameterization methods for simulation models that describe surface processes from the scale of vegetation leaves up to scales appropriate to satellite remote sensing.

Brutsaert, Wilfried↗

Data for A Hybrid Biophysical-Machine Learning Framework for Diurnal Surface Energy Flux Estimation Using Proximal Sensing

Thermal infrared-based remote sensing of surface energy fluxes has traditionally relied on high spatial resolution satellite data with revisit frequencies on the order of weeks. In this study, we evaluate a biophysics-based analytical surface energy balance model for predicting latent energy (LE) and sensible heat (H) fluxes using proximal sensing observations. The Surface Temperature Initiated Closure (STIC1.2) model has been extensively validated across a wide range of spatial and temporal scales using various satellite-derived thermal infrared data sets. Here we extend this validation by applying STIC at sub-hourly temporal resolution over multiple growing seasons for four distinct agricultural systems. We further develop and evaluate novel STIC variants that incorporate machine learning (ML) techniques to eliminate the need for surface energy balance observations, specifically net radiation and soil heat flux, thereby enhancing model applicability in data-sparse settings. The integration of a ML component to estimate surface available energy is shown to have strong predictive performance for both LE (R2 = 0.81–0.94) and H (R2 = 0.46–0.72) across all agricultural systems examined here, demonstrating the potential of hybrid biophysical-machine learning approaches for surface energy balance modeling with minimal data requirements. This study concludes with a novel application of explainable machine learning (exML) to diagnose sources of model error. This exML framework attributes residual prediction errors to both model input variables and environmental drivers not explicitly included in the simulation experiments. This approach provides a new pathway for improving model design and integrating previously overlooked yet influential variables into future model iterations.

AI/ML↗

BOREAS Level-2 MAS Surface Reflectance and Temperature Images in BSQ Format

The BOReal Ecosystem-Atmosphere Study (BOREAS) Staff Science Aircraft Data Acquisition Program focused on providing the research teams with the remotely sensed aircraft data products they needed to compare and spatially extend point results. The MODIS Airborne Simulator (MAS) images, along with other remotely sensed data, were collected to provide spatially extensive information over the primary study areas. This information includes biophysical parameter maps such as surface reflectance and temperature. Collection of the MAS images occurred over the study areas during the 1994 field campaigns. The level-2 MAS data cover the dates of 21-Jul-1994, 24-Jul-1994, 04-Aug-1994, and 08-Aug-1994. The data are not geographically/geometrically corrected; however, files of relative X and Y coordinates for each image pixel were derived by using the C130 navigation data in a MAS scan model. The data are provided in binary image format files.

Hall, Forrest G.↗

Multi-modal characterization and simulation of human epileptic circuitry

Temporal lobe epilepsy is the fourth most common neurological disorder, with about 40% of patients not responding to pharmacological treatment. Increased cellular loss is linked to disease severity and pathological phenotypes such as heightened seizure propensity. While the hippocampus is the target of therapeutic interventions, the impact of the disease at the cellular level remains unclear. Here, we show that hippocampal granule cells change with disease progression as measured in living, resected hippocampal tissue excised from patients with epilepsy. We show that granule cells increase excitability and shorten response latency while also enlarging in cellular volume and spine density. Single-nucleus RNA sequencing combined with simulations ascribes the changes to three conductances: BK, Cav2.2, and Kir2.1. In a network model, we show that these changes related to disease progression bring the circuit into a more excitable state, while reversing them produces a less excitable, “early-disease-like” state.

60 APPLIED LIFE SCIENCES↗