Search NASA⌕ Search

SEARCH · Search NASA

Results for “native state model”

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 19 records

Utah FORGE Phase 3 Native State Model: 2022 Update

This is the Phase 3 native state model update. The Phase 3 numerical model represents a significant subsurface volume below the FORGE site footprint. The model domain of 4.0 km x 4.0 km x 4.2 km is located approximately between depths of 4000 to 4200 meters below land surface. This data archive consists of 10 files, 4 of which are simulation input files and the remaining 6 are simulation output files. There is an included readme.txt file that contains details on each of the data files. The input files include meshes, FALCON code inputs, tabulated data of water properties, temperature values, and model boundaries. The output files include simulation outfiles and point data of modeled material properties.

15 GEOTHERMAL ENERGY↗

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2023 Circulation Test Data

Over the past several years, many new data sets have become available regarding the characterization of the Utah FORGE reservoir. These include, but are not limited to, the stimulation of Well 16A, the drilling and completion of Well 16B, and interwell circulation confirmatory testing. As part of the characterization efforts, conceptual models of the reservoir are re-examined as new data become available. As part of the planning for FORGE activities, numerical models are often used to predict the reservoir response to the planned testing. Stochastic methods are often employed to bound uncertainty and allow for evaluation of comprehensive ranges of key reservoir parameters. For the most recent interwell circulation confirmatory testing (July 2023), a priori numerical model predictions did bound the observed behavior (Xinj et al., 2023), but key deviations from expected behavior prompted the FORGE team to reevaluate our conceptual model of the reservoir. In early October 2023, key members of the development, testing, and monitoring teams met for 2 days to review newly collected data and discuss ‘interesting’ or ‘key’ observations. From these discussions, 15 Key Observations were documented, with several significant ones being that the discrete fracture network developed from the 16A stimulation data may not be appropriate and that the early time pressure data obtained during the summer 2023 reservoir testing were best described using radial solutions. In July 2023, two campaigns of interwell confirmatory testing were conducted, the first set of tests occurred on July 4-5 and the second set on July 18-19. The second set of circulation tests conducted at the Utah FORGE site between the injection well 16A(78)-32 and production well 16B(78)-32 on July 18 and 19, 2023 are used to calibrate material properties in a thermal-hydraulic-mechanical (THM) simulation of the discrete fracture network connecting the wells. The spatially and temporally varying reservoir properties are calibrated to match the time dependent pressure and production profiles from the circulation tests. In future work, this calibrated model will be coupled to the native state THM model of the FORGE reservoir to predict surface deformation and strains resulting from pumping schedules.

58 GEOSCIENCES↗

Numerically Testing Conceptual Models of the Utah FORGE Reservoir Using July 2024 Circulation Test Data

Over the past several years, many new data sets have become available regarding the characterization of the Utah FORGE reservoir. These include, but are not limited to, the stimulation of Well 16A, the drilling and completion of Well 16B, and interwell circulation confirmatory testing. As part of the characterization efforts, conceptual models of the reservoir are re-examined as new data become available. As part of the planning for FORGE activities, numerical models are often used to predict the reservoir response to the planned testing. Stochastic methods are often employed to bound uncertainty and allow for evaluation of comprehensive ranges of key reservoir parameters. For the most recent interwell circulation confirmatory testing (July 2023), a priori numerical model predictions did bound the observed behavior (Xinj et al., 2023), but key deviations from expected behavior prompted the FORGE team to reevaluate our conceptual model of the reservoir. In early October 2023, key members of the development, testing, and monitoring teams met for 2 days to review newly collected data and discuss ‘interesting’ or ‘key’ observations. From these discussions, 15 Key Observations were documented, with several significant ones being that the discrete fracture network developed from the 16A stimulation data may not be appropriate and that the early time pressure data obtained during the summer 2023 reservoir testing were best described using radial solutions. In July 2023, two campaigns of interwell confirmatory testing were conducted, the first set of tests occurred on July 4-5 and the second set on July 18-19. The second set of circulation tests conducted at the Utah FORGE site between the injection well 16A(78)-32 and production well 16B(78)-32 on July 18 and 19, 2023 are used to calibrate material properties in a thermal-hydraulic-mechanical (THM) simulation of the discrete fracture network connecting the wells. The spatially and temporally varying reservoir properties are calibrated to match the time dependent pressure and production profiles from the circulation tests. In future work, this calibrated model will be coupled to the native state THM model of the FORGE reservoir to predict surface deformation and strains resulting from pumping schedules.

15 GEOTHERMAL ENERGY↗

Unfolding of Proteins: Thermal and Mechanical Unfolding

We have employed a Hamiltonian model based on a self-consistent Gaussian appoximation to examine the unfolding process of proteins in external - both mechanical and thermal - force elds. The motivation was to investigate the unfolding pathways of proteins by including only the essence of the important interactions of the native-state topology. Furthermore, if such a model can indeed correctly predict the physics of protein unfolding, it can complement more computationally expensive simulations and theoretical work. The self-consistent Gaussian approximation by Micheletti et al. has been incorporated in our model to make the model mathematically tractable by signi cantly reducing the computational cost. All thermodynamic properties and pair contact probabilities are calculated by simply evaluating the values of a series of Incomplete Gamma functions in an iterative manner. We have compared our results to previous molecular dynamics simulation and experimental data for the mechanical unfolding of the giant muscle protein Titin (1TIT). Our model, especially in light of its simplicity and excellent agreement with experiment and simulation, demonstrates the basic physical elements necessary to capture the mechanism of protein unfolding in an external force field.

Hur, Joe S.↗

Impacts of Spatial Resolution in a High-Fidelity Capacity Expansion Model: An ERCOT Case Study

Capacity expansion models are important tools in examining the evolution of the electric power sector. Embedded in these tools are many modeling choices with consequential impacts on computational burden and associated analysis. In this study, we adjust the spatial resolution of the Regional Energy Deployment System (ReEDS) to understand the implications of higher-fidelity modeling on energy system projections and model solve times. The native ReEDS regions capture the contiguous United States in 134 balancing areas whereas the regions in the higher-resolution version are defined by over 3,000 U.S. counties. Using both resolutions, we conduct a case study of the Texas Interconnection (The Electric Reliability Council of Texas [ERCOT]) to explore differences in model projections and to inform appropriate applications of high spatial resolution in a large-scale, applied capacity expansion model.

county↗

Mixed lipid bilayers enable enhanced stability and activity retention of Lipase A in low-pH environments

Lipase A (LipA) from Bacillus subtilis is a versatile and industrially relevant enzyme, but its activity is compromised under acidic conditions due to aggregation and deactivation. In this study, we investigated the use of mixed lipid bilayers to stabilize and improve activity retention of LipA at low pH. Attachment to lipid bilayers, particularly those containing cationic lipids, greatly enhanced the long-term stability of LipA in acidic conditions while also leading to improvements in activity. Tethering to cationic bilayers not only shifted the apparent pH activity profile toward more acidic conditions, but also significantly enhanced activity retention upon incubation at pH 6. Notably, this protective effect persisted even without direct tethering, indicating that reversible, non-covalent interactions with the bilayer surface are sufficient for long-term stability. Circular dichroism further revealed that the secondary structure of LipA was retained, while dynamic light scattering suggested that activity loss in solution was primarily due to aggregation of the native state. Together, these findings support a model in which lipid bilayers mitigate aggregation and stabilize LipA through transient interactions and electrostatic modulation of the local surface environment. Furthermore, this approach exemplifies a simple, tether-free strategy for enhancing enzyme performance in acidic or destabilizing conditions, with implications for biocatalysis, biosensing, and therapeutic delivery.

Biocatalysis↗

Cryo-$\mathrm{EM}$ reveals an entangled kinetic trap in the folding of a catalytic $\mathrm{RNA}$

Functional RNAs fold through complex pathways that can contain misfolded “kinetic traps.” A complete model of RNA folding requires understanding the formation of these misfolded states, but they are difficult to characterize because of their transient and potentially conformationally dynamic nature. We used cryo–electron microscopy (cryo-EM) to visualize a long-lived misfolded state in the folding pathway of the Tetrahymena thermophilagroup I intron, a paradigmatic RNA structure-function model system. The structure revealed how this state forms native-like secondary structure and tertiary contacts but contains two incorrectly crossed strands, consistent with a previous model. This incorrect topology mispositions a critical catalytic domain and cannot be resolved locally as extensive refolding is required. This work provides a structural framework for interpreting decades of biochemical and functional studies and demonstrates the power of cryo-EM for the exploration of RNA folding pathways.

59 BASIC BIOLOGICAL SCIENCES↗

Native Architecture of Wheat Straw Cell Walls: A Unified Model from X-ray Scattering and Solid-State NMR

Plant secondary cell walls constitute the dominant reservoir of renewable biomass, comprising tightly packed cellulose, hemicellulose, and lignin at the nanoscale. Recent advances in solid-state NMR spectroscopy and the availability of small-angle X-ray scattering for biomass characterization have led to an accumulation of experimental data on cell wall organization, yet no explicit structure model has simultaneously satisfied both Xray and NMR observations. Using wheat straw as a model system, we propose a structural framework consistent with current knowledge of cellulose biosynthesis, X-ray scattering data, and one- and two-dimensional 13 C solid-state NMR spectra. In this model, 18-chain elementary fibrils align in parallel and populate the cross-section at random. Arabinose-substituted xylan shows no conformational dependence for cellulose-binding in wheat, and only a minor fraction of 2-fold xylan appears in close proximity to cellulose, unlike in Arabidopsis, where xylan is more tightly attached to the cellulose surface. While NMR data cannot unambiguously resolve the internal arrangement of the 18 glucan chains, X-ray scattering profiles uniquely constrain the fibril size and exclude the possibility of tight bundling in the intact walls. The specific interaction between the matrix polymers and the cellulose elementary fibrils must be reconsidered in light of the small interfibril spaces, which bring the matrix components into spatial proximity with cellulose even in the absence of attractive interactions. These findings provide fundamental molecular-level insight into cellulose fibril architecture and matrix−polymer interactions, resolving longstanding discrepancies between spectroscopic and scattering data and advancing our understanding of biopolymer assembly into structurally and functionally versatile lignocellulosic biomaterials.

Carbohydrates↗

Nuclear Power's Future Role in a Decarbonized U.S. Electricity System

This study explores the potential future role of nuclear energy in a decarbonized U.S. electricity system through a multi-model comparison approach. We employ four state-of-the-art CEMs with native and harmonized input assumptions, layered with different policy and technology trajectories. Comparing outputs across models, technology assumptions, and policy scenarios informs model understanding, interpretation, and development decisions. Under current policies, models differ in their projections for nuclear retirements, but nuclear power plants consistently run with high capacity factors and new builds only occur in scenarios with very low nuclear costs. String power sector carbon policies drive models to align in keeping existing nuclear capacity and employing nuclear plant flexibility, but they may not be enough to bring new nuclear capacity online in the absence of significant cost declines. Therefore, significant economic deployment of new nuclear capacity requires both a stringent electric sector CO2 policy and very low cost assumptions for new nuclear. While these scenarios should not be interpreted as predictions, they are informative for understanding differing model assessments of the relative competitiveness of nuclear energy under a range of policy and technology conditions.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Evaluation of native Earth system model output with ESMValTool v2.6.0

Earth system models (ESMs) are state-of-the-art climate models that allow numerical simulations of the past, present-day, and future climate. To extend our understanding of the Earth system and improve climate change projections, the complexity of ESMs heavily increased over the last decades. As a consequence, the amount and volume of data provided by ESMs has increased considerably. Innovative tools for a comprehensive model evaluation and analysis are required to assess the performance of these increasingly complex ESMs against observations or reanalyses. One of these tools is the Earth System Model Evaluation Tool (ESMValTool), a community diagnostic and performance metrics tool for the evaluation of ESMs. Input data for ESMValTool needs to be formatted according to the CMOR (Climate Model Output Rewriter) standard, a process that is usually referred to as “CMORization”. While this is a quasi-standard for large model intercomparison projects like the Coupled Model Intercomparison Project (CMIP), this complicates the application of ESMValTool to non-CMOR-compliant climate model output. In this paper, we describe an extension of ESMValTool introduced in v2.6.0 that allows seamless reading and processing of “native” climate model output, i.e., operational output produced by running the climate model through the standard workflow of the corresponding modeling institute. This is achieved by an extension of ESMValTool's preprocessing pipeline that performs a CMOR-like reformatting of the native model output during runtime. Thus, the rich collection of diagnostics provided by ESMValTool is now fully available for these models. For models that use unstructured grids, a further preprocessing step required to apply many common diagnostics is regridding to a regular latitude–longitude grid. Extensions to ESMValTool's regridding functions described here allow for more flexible interpolation schemes that can be used on unstructured grids. Currently, ESMValTool supports nearest-neighbor, bilinear, and first-order conservative regridding from unstructured grids to regular grids. Example applications of this new native model support are the evaluation of new model setups against predecessor versions, assessing of the performance of different simulations against observations, CMORization of native model data for contributions to model intercomparison projects, and monitoring of running climate model simulations. For the latter, new general-purpose diagnostics have been added to ESMValTool that are able to plot a wide range of variable types. Currently, five climate models are supported: CESM2 (experimental; at the moment, only surface variables are available), EC-Earth3, EMAC, ICON, and IPSL-CM6. As the framework for the CMOR-like reformatting of native model output described here is implemented in a general way, support for other climate models can be easily added.

58 GEOSCIENCES↗

Purification of the small mechanosensitive channel of Escherichia coli (MscS): the subunit structure, conduction, and gating characteristics in liposomes

The small mechanosensitive channel, MscS, is a part of the turgor-driven solute efflux system that protects bacteria from lysis in the event of osmotic downshift. It has been identified in Escherichia coli as a product of the orphan yggB gene, now called mscS (Levina et al., 1999, EMBO J. 18:1730). Here I show that that the isolated 31-kDa MscS protein is sufficient to form a functional mechanosensitive channel gated directly by tension in the lipid bilayer. MscS-6His complexes purified in the presence of octylglucoside and lipids migrate in a high-resolution gel-filtration column as particles of approximately 200 kDa. Consistent with that, the protein cross-linking patterns predict a hexamer. The channel reconstituted in soybean asolectin liposomes was activated by pressures of 20-60 mm Hg and displayed the same asymmetric I-V curve and slight anionic preference as in situ. At the same time, the single-channel conductance is proportional to the buffer conductivity in a wide range of salt concentrations. The rate of channel activation in response to increasing pressure gradient across the patch was slower than the rate of closure in response to decreasing steps of pressure gradient. Therefore, the open probability curves were recorded with descending series of pressures. Determination of the curvature of patches by video imaging permitted measurements of the channel activity as a function of membrane tension (gamma). Po(gamma) curves had the midpoint at 5.5 +/- 0.1 dyne/cm and gave estimates for the energy of opening DeltaG = 11.4 +/- 0.5 kT, and the transition-related area change DeltaA = 8.4 +/- 0.4 nm(2) when fitted with a two-state Boltzmann model. The correspondence between channel properties in the native and reconstituted systems is discussed.

NASA Discipline Cell Biology↗

System Modeling of Metabolic Heat Regenerated Temperature Swing Adsorption (MTSA) Subassembly for Prototype Design

This paper describes modeling methods for the three core components of a Metabolic heat regenerated Temperature Swing Adsorption (MTSA) subassembly: the sorbent bed, a sublimation (cooling) heat exchanger (SHX), and a condensing icing (warming) heat exchanger (CIHX). The primary function of the MTSA, removing carbon dioxide from a ventilation loop, is performed via the sorbent bed. The CIHX is used to heat the sorbent bed for desorption and to remove moisture from the ventilation loop while the SHX is alternately employed to cool the sorbent bed via sublimation of a spray of water at low pressure to prepare the reconditioned bed for the next cycle. This paper describes a system level model of the MTSA as developed in Thermal Desktop and SINDA/FLUINT including assumptions on geometry and physical phenomena, modeling methodology and relevant pa ra mete rizatio ns. Several areas of particular modeling interest are discussed. In the sorbent bed, capture of the translating CO2 saturation front and associated local energy and mass balance in both adsorbing and desorbing modes is covered. The CIHX poses particular challenges for modeling in SINDA/FLUINT as accounting for solids states in fluid submodels are not a native capability. Methods for capturing phase change and latent heat of ice as well as the transport properties across a layer of low density accreted frost are developed. This extended modeling capacity is applicable to temperatures greater than 258 K. To extend applicability to the minimum device temperature of 235 K, a method for a mapped transformation of temperatures from below the limit temperatures to some value above is given along with descriptions for associated material property transformations and the resulting impacts to total heat and mass transfer. Similar considerations are shown for the SHX along with assumptions for flow mechanics and resulting model methods for sublimation in a flow.

Bower, Chad↗

Modeling Of Metabolic Heat Regenerated Temperature Swing Adsorption (MTSA) Subassembly For Prototype Design

This paper describes modeling methods for the three core components of a Metabolic heat regenerated Temperature Swing Adsorption (MTSA) subassembly: a sorbent bed, a sublimation (cooling) heat exchanger (SHX), and a condensing icing (warming) heat exchanger (CIHX). The primary function of the MTSA, removing carbon dioxide from a space suit Portable Life Support System (PLSS) ventilation loop, is performed via the sorbent bed. The CIHX is used to heat the sorbent bed for desorption and to remove moisture from the ventilation loop while the SHX is alternately employed to cool the sorbent bed via sublimation of a spray of water at low pressure to prepare the reconditioned bed for the next cycle. This paper describes subsystem heat a mass transfer modeling methodologies relevant to the description of the MTSA subassembly in Thermal Desktop and SINDA/FLUINT. Several areas of particular modeling interest are discussed. In the sorbent bed, capture of the translating carbon dioxide (CO2) front and associated local energy and mass balance in both adsorbing and desorbing modes is covered. The CIHX poses particular challenges for modeling in SINDA/FLUINT as accounting for solids states in fluid submodels are not a native capability. Methods for capturing phase change and latent heat of ice as well as the transport properties across a layer of low density accreted frost are developed. This extended modeling capacity is applicable to temperatures greater than 258 K. To extend applicability to the minimum device temperature of 235 K, a method for a mapped transformation of temperatures from below the limit temperatures to some value above is given along with descriptions for associated material property transformations and the resulting impacts to total heat and mass transfer. Similar considerations are given for the SHX along with functional relationships for areal sublimation rates as limited by flow mechanics in t1he outlet duct.

Bower, Chad E.↗

Designing Antifouling and Antimicrobial Interfaces: Structural Characterization using CryoEM, Automated Microscopy, and AI Image Segmentation

The design of functionalized surfaces for interactions with biological systems is critical across sectors such as healthcare, energy, and agriculture. Tailoring materials for specific applications, such as antifouling and antimicrobial surfaces, demands a comprehensive understanding of topology and chemistry across multiple length and time scales on both biological and materials systems. This work presents the development and characterization of nanostructured surfaces with controlled topographies and chemistries that enhance bacterial membrane disruption, reduce biofilm formation, and improve antimicrobial and antifouling capabilities. Two specific use cases will be presented - the use of cellulose nanocrystals (CNCs) for bacterial growth inhibition and the development of antifouling surfaces to prevent protein and bacterial adsorption [1-4]. By leveraging large language models (LLMs) for image segmentation and training [5], we enable automated analysis of terabyte-scale cryogenic electron microscopy (cryoEM) datasets. This analysis provides statistical insights into the biotic/abiotic interface and facilitates automated electron microscopy experiments to mitigate time and dose. The integration of cryogenic electron tomography (cryoET) and cryogenic focused ion beam (cryoFIB) milling enables high-resolution, near-native-state imaging and 3D reconstructions of bio/material interfaces [6]. Orthogonal characterization techniques and computational modeling further enhances our understanding, offering a robust platform for the design and optimization of next-generation functional surfaces [7].

Williams, Alexis [ORNL] (ORCID:0000000252835822)↗

Ventenata dubia projected to expand in the western United States despite future novel conditions

Abstract Distributions of both native and invasive species are expected to shift under future climate. Species distribution models (SDMs) are often used to explore future habitats, but sources of uncertainty including novel climate conditions may reduce the reliability of future projections. We explore the potential spread of the invasive annual grass ventenata ( Ventenata dubia ) in the western United States under both current and future climate scenarios using boosted regression tree models and 30 global climate models (GCMs). We quantify novel climate conditions, prediction variability arising from both the SDMs and GCMs, and the agreement among GCMs. Results demonstrate that currently suitable habitat is concentrated inside the invaded range of the northwest, but substantial habitat exists outside the invaded range in the Southern Rockies and southwestern US mountains. Future suitability projections vary greatly among GCMs, but GCMs commonly projected decreased suitability in the invaded range and increased suitability along higher elevations of interior mountainous areas. Climate novelty did not appear to undermine the prediction reliability in many cases where the climate–species relationship was fully represented by the occurrence data. GCM‐derived variability resulting from variation in future cool season precipitation and temperature seasonality was greatest in the Rocky Mountains. SDM‐derived variability was higher in currently suitable habitat, and few GCMs projections agreed that these areas would contain future suitable habitat. However, while prediction variability was high, many GCM projections agreed that parts of the Rocky, Wasatch, and Uinta Mountains would contain highly suitable habitat in the future. As disturbances in the interior mountains occur in coming decades, reducing some natural barriers to invasion, land managers, and conservationists will need to monitor for ventenata in post‐disturbance environments. Changes to invasion potential may not play out for several decades, but results related to current potential may have applications for early detection and rapid response planning.

Environmental Sciences & Ecology↗

Hybrid Oscillator-Qubit Quantum Processors: Instruction Set Architectures, Abstract Machine Models, and Applications

This tutorial offers a pedagogical guide to hybrid quantum processors that integrate discrete-variable (DV) qubits and continuous-variable (CV) oscillators. Aimed at computer scientists, engineers, and physicists, it provides an overview of the experimental, algorithmic, and architectural aspects of this novel and rapidly developing hardware model. Experimental realizations of this model include superconducting, trapped-ion, and neutral-atom platforms. By combining DV and CV components, hybrid oscillator-qubit processors enable a powerful new paradigm that offers complementary strengths for quantum control, error correction, computation, and simulation. Working toward the goal of a full-stack system connecting applications to CV-DV hardware, we define and formulate abstract machine models and instruction set architectures. These essential abstractions enable codesign of hardware and software, and resource estimation for exploring the potential of current and future hardware for computational and simulation tasks. Using these abstractions, we present both new and existing examples that illustrate the benefits of hybrid CV-DV processors relative to traditional DV-only hardware in computation as well as quantum simulation of physical models. Examples include algorithms for transferring states between DV and CV systems, performing the quantum Fourier transform, and simulation of lattice gauge theories. Relative to qubit-only hardware, the bosonic degrees of freedom natively available in hybrid architectures can substantially reduce the circuit complexity of simulations for physical models containing bosons. A key technique is the extension of quantum signal processing ideas to CV-DV systems. This work is intended to serve as a timely and comprehensive guide to this relatively unexplored yet promising approach to quantum computation and to provide a road map to guide future development.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Structure and identification of the native PLP synthase complex from Methanosarcina acetivorans lysate

Many protein-protein interactions behave differently in biochemically purified forms as compared to their in vivo states. As such, determining native protein structures may elucidate structural states previously unknown for even well-characterized proteins. Here, we apply the bottom-up structural proteomics method, cryoID , toward a model methanogenic archaeon. While they are keystone organisms in the global carbon cycle and active members of the human microbiome, there is a general lack of characterization of methanogen enzyme structure and function. Through the cryoID approach, we successfully reconstructed and identified the native Methanosarcina acetivorans pyridoxal 5′-phosphate (PLP) synthase (PdxS) complex directly from cryogenic electron microscopy (cryo-EM) images of fractionated cellular lysate. We found that the native PdxS complex exists as a homo-dodecamer of PdxS subunits, and the previously proposed supracomplex containing both the synthase (PdxS) and glutaminase (PdxT) was not observed in cellular lysate. Our structure shows that the native PdxS monomer fashions a single 8α/8β TIM-barrel domain, surrounded by seven additional helices to mediate solvent and interface contacts. A density is present at the active site in the cryo-EM map and is interpreted as ribose 5-phosphate. In addition to being the first reconstruction of the PdxS enzyme from a heterogeneous cellular sample, our results reveal a departure from previously published archaeal PdxS crystal structures, lacking the 37-amino-acid insertion present in these prior cases. This study demonstrates the potential of applying the cryoID workflow to capture native structural states at atomic resolution for archaeal systems, for which traditional biochemical sample preparation is nontrivial.

Methanosarcina acetivorans↗

Matrix Approach to Accelerate Spin-Up of CLM5

Numerical models have been developed to investigate and understand responses of biogeochemical cycle to global changes. Steady state, when a system is in dynamic equilibrium, is generally required to initialize these model simulations. However, the spin-up process that is used to achieve steady state pose a great burden to computational resources, limiting the efficiency of global modeling analysis on biogeochemical cycles. This study introduces a new Semi-Analytical Spin-Up (SASU) to tackle this grand challenge. We applied SASU to Community Land Model version 5 and examined its computational efficiency and accuracy. At the Brazil site, SASU is computationally 7 times more efficient than (or saved up to 86% computational cost in comparison with) the traditional native dynamics (ND) spin-up to reach the same steady state. Globally, SASU is computationally 8 times more efficient than the accelerated decomposition spin-up and 50 times more efficient than ND. In summary, SASU achieves the highest computational efficiency for spin-up on site and globally in comparison with other spin-up methods. It is generalizable to wide biogeochemical models and thus makes computationally costly studies (e.g., parameter perturbation ensemble analysis and data assimilation) possible for a better understanding of biogeochemical cycle under climate change.

54 ENVIRONMENTAL SCIENCES↗