Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic profiling”

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 883 records · Page 49

Dynamic Modeling for the Separation of Rare Earth Elements Using Solvent Extraction: Predicting Separation Performance Using Laboratory Equilibrium Data

Industrial rare earth element (REE) separation facilities utilize acidic cation exchange ligands such as 2-ethylhexylphosphonic acid mono-2-ethylhexyl ester (PC88A) for solvent extraction processes. REE separations are costly and difficult due to their chemical similarities and subsequent low separation factors. Several empirical correlations are available in the literature to predict steady state extraction equilibria for various solvent systems. However, complete solvent extraction flowsheet design for REE separations requires complex scrubbing and stripping circuits to separate and produce individual pure species. Furthermore, the sensitive process chemistry is susceptible to process upsets and may manifest as irreversible transients that complicate startup and continuous operation of a solvent extraction plant. Thus, dynamic modeling of REE separations will aid in process design, optimization, and management of process fluctuations. A time-dependent MATLAB/SIMULINK modeling framework has been developed for REE solvent extraction processes using laboratory equilibrium data. The model was used to develop a flowsheet that produced high purity neodymium from a 25 wt% praseodymium and 75 wt% neodymium feed. Laboratory mixer-settlers were used to verify and validate model performance. Results indicated that the model reasonably predicts the dynamic behavior of a counter-current REE separation process, and accurately predicts the steady state REE concentration profiles across the cascade. The modeling framework can be applied to any REE separation or solvent system provided adequate laboratory equilibrium data are available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rapidly Assessing the Quality of Targeted Proteomics Experiments through Monitoring Stable-Isotope Labeled Standards

Selected reaction monitoring (SRM), also known as multiple reaction monitoring (MRM), is a data acquisition technique used in targeted analysis of molecules, including targeted proteomic studies. It exploits the unique capability of triple-quadrupole (QQQ) mass spectrometers to monitor the predefined precursor and fragment ion pairs of specific molecules of interest throughout a liquid chromatography (LC) elution profile. Compared to shotgun proteomics, targeted proteomics using SRM has high selectivity, high sensitivity, and a wide linear dynamic range,(1-3) which makes it especially useful in the accurate and reproducible quantification of low-abundance proteins in highly complex biological samples. SRM has been widely used in the fields of biomarker discovery,(4-7) analysis of protein post-translational modifications(8) and characterization of biological protein networks.(4,9) In the recent years, multiple technical advances have greatly improved the throughput of SRM analyses, allowing for the quantification of hundreds of peptides in a single analysis.(6,7,10) For example, a single 800-plex SRM assay (e.g., 400 unlabeled and heavy-labeled peptide pairs and 2400 transitions with retention time scheduling) using ultrahigh-performance liquid chromatography (UHPLC) has been developed to quantify proteins in plasma.(11) Advanced labeling techniques utilizing in vitro proteome-assisted MRM for protein absolute quantification (iMPAQT) demonstrated the capability of SRM in genome-wide protein quantification of over 18?000 human proteins.(12) Moreover, the scan speed of QQQ instrumentation has been greatly improved in recent releases of commercial instrumentation. The newly developed TSQ Altis (released in 2017) can scan more than 600 transitions per second, which is 6 times more than a traditional QQQ scan speed of 100 transitions per second. The breadth of measurement enabled by these technological improvements to QQQ mass spectrometry has increased the feasibility and popularity of large scale targeted proteomics studies,(13) a major application of which will be in clinical studies in which up to hundreds of protein candidates need to be quantified in hundreds of clinical specimens.(14) Quantitative accuracy is a primary motivating factor for utilizing a targeted proteomics protocol. The precise and reproducible absolute quantification produced by SRM assays is essential to many clinical and laboratory experiments.(15,16) Because the abundance of an endogenous peptide is calculated from the measurement of the spiked-in reference standard, it is essential to assess the data quality of these references.(17) In early applications of targeted proteomics, when instrument speed greatly limited the number of transitions that could be monitored, much of this quality assessment was done manually. However, recent improvements in instrument performance and experimental design have enabled a dramatic increase in the number of target peptides and associated SRM transitions, which makes manual quality assessment an untenable and laborious task. A variety of computational tools assist in SRM experiment design and data analysis. The first task in creating an SRM experiment is the choice of proteins and representative peptides to monitor. Achieving a reliable protein abundance requires appropriately choosing peptides that have a strong signal and are free from interferences in the biological matrix. Numerous computational tools exist to facilitate assay design by identifying peptides and refining SRM transitions.(18-21) To help share these assays and eliminate time spent designing the same transitions at multiple institutions, community portals have begun to host well-designed and vetted assays.(22-24) Analyzing the experimental data requires significant computational effort to align files across replicates and experimental conditions, pick peaks and produce quantitative values, normalize data and perform statistical tests, etc.(25-27)

Gibbons, Bryson C.↗

Molecular Mechanism for Isoform-Selective Inhibition of Acyl Protein Thioesterases 1 and 2 (APT1 and APT2)

Post-translational S-palmitoylation directs the trafficking and membrane localization of hundreds of cellular proteins, often involving a coordinated palmitoylation cycle that requires both protein acyl transferases (PATs) and acyl protein thioesterases (APTs) to actively redistribute S-palmitoylated proteins toward different cellular membrane compartments. This process is necessary for the trafficking and oncogenic signaling of S-palmitoylated Ras isoforms, and potentially many peripheral membrane proteins. The depalmitoylating enzymes APT1 and APT2 are separately conserved in all vertebrates, suggesting unique functional roles for each enzyme. The recent discovery of the APT isoform-selective inhibitors ML348 and ML349 has opened new possibilities to probe the function of each enzyme, yet it remains unclear how each inhibitor achieves orthogonal inhibition. Herein, we report the high-resolution structure of human APT2 in complex with ML349 (1.64 Å), as well as the complementary structure of human APT1 bound to ML348 (1.55 Å). Although the overall peptide backbone structures are nearly identical, each inhibitor adopts a distinct conformation within each active site. In APT1, the trifluoromethyl group of ML348 is positioned above the catalytic triad, but in APT2, the sulfonyl group of ML349 forms hydrogen bonds with active site resident waters to indirectly engage the catalytic triad and oxyanion hole. Reciprocal mutagenesis and activity profiling revealed several differing residues surrounding the active site that serve as critical gatekeepers for isoform accessibility and dynamics. Structural and biochemical analysis suggests the inhibitors occupy a putative acyl-binding region, establishing the mechanism for isoform-specific inhibition, hydrolysis of acyl substrates, and structural orthogonality important for future probe development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Extraction and biomolecular analysis of dermal interstitial fluid collected with hollow microneedles

Dermal interstitial fluid (ISF) is an underutilized information-rich biofluid potentially useful in health status monitoring applications whose contents remain challenging to characterize. Here, we present a facile microneedle approach for dermal ISF extraction with minimal pain and no blistering for human subjects and rats. Extracted ISF volumes were sufficient for determining transcriptome, and proteome signatures. We noted similar profiles in ISF, serum, and plasma samples, suggesting that ISF can be a proxy for direct blood sampling. Dynamic changes in RNA-seq were recorded in ISF from induced hypoxia conditions. Finally, we report the first isolation and characterization, to our knowledge, of exosomes from dermal ISF. The ISF exosome concentration is 12–13 times more enriched when compared to plasma and serum and represents a previously unexplored biofluid for exosome isolation. This minimally invasive extraction approach can enable mechanistic studies of ISF and demonstrates the potential of ISF for real-time health monitoring applications.

59 BASIC BIOLOGICAL SCIENCES↗

An EQT-cDFT approach to determine thermodynamic properties of confined fluids

We present a continuum-based approach to predict the structure and thermodynamic properties of confined fluids at multiple length-scales, ranging from a few angstroms to macro-meters. The continuum approach is based on the empirical potential-based quasi-continuum theory (EQT) and classical density functional theory (cDFT). EQT is a simple and fast approach to predict inhomogeneous density and potential profiles of confined fluids. We use EQT potentials to construct a grand potential functional for cDFT. The EQT-cDFT-based grand potential can be used to predict various thermodynamic properties of confined fluids. In this work, we demonstrate the EQT-cDFT approach by simulating Lennard-Jones fluids, namely, methane and argon, confined inside slit-like channels of graphene. We show that the EQT-cDFT can accurately predict the structure and thermodynamic properties, such as density profiles, adsorption, local pressure tensor, surface tension, and solvation force, of confined fluids as compared to the molecular dynamics simulation results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics of gradient formation by intracellular shuttling

A number of important cellular functions rely on the formation of intracellular protein concentration gradients. Experimental studies discovered a number of mechanisms for the formation of such gradients. One of the mechanisms relies on the intracellular shuttling of a protein that interconverts between the two states with different diffusivities, under the action of two enzymes, one of which is localized to the plasma membrane, whereas the second is uniformly distributed in the cytoplasm. Recent work reported an analytical solution for the steady state gradient in this mechanism, obtained in the framework of a one-dimensional reaction-diffusion model. Here, we study the dynamics in this model and derive analytical expressions for the Laplace transforms of the time-dependent concentration profiles in terms of elementary transcendental functions. Inverting these transforms numerically, one can obtain time-dependent concentration profiles of the two forms of the protein.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulation of diatomic gas-wall interaction and accommodation coefficients for negative ion sources and accelerators

Particle-wall interactions determine in different ways the operating conditions of plasma sources, ion accelerators, and beams operating in vacuum. For instance, a contribution to gas heating is given by ion neutralization at walls; beam losses and stray particle production—detrimental for high current negative ion systems such as beam sources for fusion—are caused by collisional processes with residual gas, with the gas density profile that is determined by the scattering of neutral particles at the walls. This paper shows that Molecular Dynamics (MD) studies at the nano-scale can provide accommodation parameters for gas-wall interactions, such as the momentum accommodation coefficient and energy accommodation coefficient: in non-isothermal flows (such as the neutral gas in the accelerator, coming from the plasma source), these affect the gas density gradients and influence efficiency and losses in particular of negative ion accelerators. For ideal surfaces, the computation also provides the angular distribution of scattered particles. Classical MD method has been applied to the case of diatomic hydrogen molecules. Single collision events, against a frozen wall or a fully thermal lattice, have been simulated by using probe molecules. Different modelling approximations are compared.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The rising stellar velocity dispersion of M87 from integrated starlight

We have measured the line-of-sight velocity distribution from integrated stellar light at two points in the outer halo of M87 (NGC 4486), the second-rank galaxy in the Virgo Cluster. The data were taken at R = 480'' (∼41.5 kpc) and R = 526'' (∼45.5 kpc) along the SE major axis. The second moment for a non-parametric estimate of the full velocity distribution is 420 ± 23 km s{sup –1} and 577 ± 35 km s{sup –1}, respectively. There is intriguing evidence in the velocity profiles for two kinematically distinct stellar components at the position of our pointing. Under this assumption, we employ a two-Gaussian decomposition and find the primary Gaussian having rest velocities equal to M87 (consistent with zero rotation) and second moments of 383 ± 32 km s{sup –1} and 446 ± 43 km s{sup –1}, respectively. The asymmetry seen in the velocity profiles suggests that the stellar halo of M87 is not in a relaxed state and confuses a clean dynamical interpretation. That said, either measurement (full or two component model) shows a rising velocity dispersion at large radii, consistent with previous integrated light measurements, yet significantly higher than globular cluster measurements at comparable radial positions. These integrated light measurements at large radii, and the stark contrast they make to the measurements of other kinematic tracers, highlight the rich kinematic complexity of environments like the center of the Virgo Cluster and the need for caution when interpreting kinematic measurements from various dynamical tracers.

79 ASTRONOMY AND ASTROPHYSICS↗

Demonstrating the Potential of Adaptive LMS Filtering on FPGA-Based Qubit Control Platforms for Improved Qubit Readout in 2D and 3D Quantum Processing Units

Advancements in quantum computing underscore the critical need for sophisticated qubit readout techniques to accurately discern quantum states. This abstract presents our research intended for optimizing readout pulse fidelity for 2D and 3D Quantum Processing Units (QPUs), the latter coupled with Superconducting Radio Frequency (SRF) cavities. Focusing specifically on the application of the Least Mean Squares (LMS) adaptive filtering algorithm, we explore its integration into the FPGA-based control systems to enhance the accuracy and efficiency of qubit state detection by improving Signal-to-Noise Ratio (SNR). Implementing the LMS algorithm on the Zynq UltraScale+ RFSoC Gen 3 devices (RFSoC 4x2 FPGA and ZCU216 FPGA) using the Quantum Instrumentation Control Kit (QICK) open-source platform, we aim to dynamically test and adjust the filtering parameters in real-time to characterize and adapt to the noise profile presented in quantum computing readout signals. Our preliminary results demonstrate the LMS filter's capability to maintain high readout accuracy while efficiently managing FPGA resources. These findings are expected to contribute to developing more reliable and scalable quantum computing architectures, highlighting the pivotal role of adaptive signal processing in quantum technology advancements.

Johnson, Hans↗

Additive Manufacture of Multi-Functional, Large-Scale Components

The poster summarizes the main results for a LDRD focused on developing a arc-based additive manufacturing (AM) system capable of fabricating large components with a composition gradient otherwise known as functionally graded materials (FGMs). FGMs are capable of addressing local extreme conditions in an optimal and cost-effective manner. The AM system (referred to as multi-Wire Arc Additive Manufacturing or m-WAAM) consisted of modifying a standard automated gas tungsten arc welding system to deliver multiple welding wires to the weld pool at controlled rates to deposit desired alloy compositions at critical locations. Numerous shapes and composition gradients were fabricated with the system and the effect of various deposition parameters on the nature of the composition profile were evaluated. Additionally, heat treatment of the FGM was performed to develop graded microstructures which should enhance dynamic impact performance of inexpensive steels. The project highlights the potential of m-WAAM to efficiently fabricate FGMs. However, this project only explored the science behind the technique and additional development work will be required to develop a fully-functional system that can build more complex components.

36 MATERIALS SCIENCE↗

Offshore Wind to Hydrogen - Modeling, Analysis, Testing and International Collaboration Work

This project explores electrolytic hydrogen production from an offshore wind turbine by: Modeling: Simulate an offshore wind turbine and generate power output profiles; Analysis: Analyze offshore conditions and determine operational modes of the stack; Testing: Perform hardware testing to evaluate dynamic characteristics; and International Collaboration Work: Share learnings and publish. Potential Outcomes: Accelerate development of an integrated, in-turbine offshore wind hydrogen system; and Support two DOE Energy Earthshots - Hydrogen and Floating Offshore Wind.

electrolyzer↗

Demonstrating the Potential of Adaptive LMS Filtering on FPGA-Based Qubit Control Platforms for Improved Qubit Readout in 2D and 3D Quantum Processing Units

Advancements in quantum computing underscore the critical need for sophisticated qubit readout techniques to accurately discern quantum states. This abstract presents our research intended for optimizing readout pulse fidelity for 2D and 3D Quantum Processing Units (QPUs), the latter coupled with Superconducting Radio Frequency (SRF) cavities. Focusing specifically on the application of the Least Mean Squares (LMS) adaptive filtering algorithm, we explore its integration into the FPGA-based control systems to enhance the accuracy and efficiency of qubit state detection by improving Signal-to-Noise Ratio (SNR). Implementing the LMS algorithm on the Zynq UltraScale+ RFSoC Gen 3 devices (RFSoC 4x2 FPGA and ZCU216 FPGA) using the Quantum Instrumentation Control Kit (QICK) open-source platform, we aim to dynamically test and adjust the filtering parameters in real-time to characterize and adapt to the noise profile presented in quantum computing readout signals. Our preliminary results demonstrate the LMS filter's capability to maintain high readout accuracy while efficiently managing FPGA resources. These findings are expected to contribute to developing more reliable and scalable quantum computing architectures, highlighting the pivotal role of adaptive signal processing in quantum technology advancements.

Johnson, Hans↗

SU-G-TeP4-15: The Roucoulette: A Set of Quality Control Tests for Dynamic Trajectory (4Pi) Treatment Delivery Techniques

Purpose: To present and validate a set of quality control tests for trajectory treatment delivery using synchronized dynamic couch (translation and rotation), MLC and collimator motion. Methods: The quality control tests are based on the Picket fence test, which consist of 5 narrow band 2mm width spaced at 2.5cm intervals, and adds progressively synchronized dynamic motions. The tests were exposed on GafChromic EBT3 films. The first test is a regular (no motion and MLC static while beam is on) Picket Fence test used as baseline. The second test includes simultaneous collimator and couch rotation, each stripe corresponding to a different rotation speed. Errors in these tests were introduced (0.5 degree and 1 degree error in rotation synchronization) to assess the error sensitivity of this test. The second test is similar to the regular Picket Fence but now including dynamic MLC motion and couch translation (including acceleration during delivery) while the beam is on. Finally in the third test, which is a combination of the first and second test, the Picket Fence pattern is delivered using synchronized collimator and couch rotation and synchronized dynamic MLC and couch translation including acceleration. Films were analyzed with FilmQA Pro. Results: The distance between the peaks in the dose profile where measured (18.5cm away from the isocentre in the inplane direction where non synchronized rotation would have the largest effect) and compared to the regular Picket Fence tests. For well synchronized motions distances between peaks where between 24.9–25.4 mm identical to the regular Picket Fence test. This range increased to 24.4–26.4mm and 23.4–26.4mm for 0.5 degree and 1 degree error respectively. The amplitude also decreased up to 15% when errors are introduced. Conclusion: We demonstrated that the Roucoulette tests can be used as a quality control tests for trajectory treatment delivery using synchronized dynamic motion.

60 APPLIED LIFE SCIENCES↗

The Onset of Nucleate Boiling Prediction in a Mini Channel with the Eulerian Multiphase Flow

The prediction of boiling incipience is a critical issue for a reactor. Vapor lock causes operating instability and undesirable temperature rise, especially for small channels. This study uses the Eulerian multiphase flow boiling model to predict the onset of nucleate boiling in a millimeter scale rectangular channel. Two wall boiling models, namely the RPI (Rensselaer Polytechnic Institute) and the non-equilibrium sub cooled model are applied. Computational fluid dynamics is employed for this study to investigate the impact of operating conditions like liquid flow rate, operating pressure, heat flux, velocity profile, turbulence model, and inlet liquid temperature. Here, the effect of fundamental boiling parameters on the surface wall temperature and vapor volume fraction is also studied. The parameters under consideration are the bubble departure diameter, bubble departure frequency, nucleate site density, quenching time period, and interface heat transfer coefficients. A series of CFD calculations is conducted by varying the considered variables systematically in a wide range of flow conditions covering laminar, transition, and turbulent flows. The sub-component heat fluxes like the liquid convective, the vapor convective, the evaporative, and the quenching terms are monitored to infer the boiling dynamics in the transition region from the single-phase flow to the nucleate boiling zone. It is found that the CFD approach to detecting the boiling incipience point agrees reasonably with available experimental data. However, its limitations, like inaccuracy for lower flow rates and early transition to boiling flow, are also noticed due to the nature of the considered RPI model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Liquid-based gating mechanism with tunable multiphase selectivity and antifouling behaviour

Living organisms make extensive use of micro- and nanometre-sized pores as gatekeepers for controlling the movement of fluids, vapours and solids between complex environments. The ability of such pores to coordinate multiphase transport, in a highly selective and subtly triggered fashion and without clogging, has inspired interest in synthetic gated pores for applications ranging from fluid processing to 3D printing and lab-on-chip systems(1-10). But although specific gating and transport behaviours have been realized by precisely tailoring pore surface chemistries and pore geometries(6,11-17), a single system capable of controlling complex, selective multiphase transport has remained a distant prospect, and fouling is nearly inevitable(11,12). Here we introduce a gating mechanism that uses a capillary-stabilized liquid as a reversible, reconfigurable gate that fills and seals pores in the closed state, and creates a non-fouling, liquid-lined pore in the open state. Theoretical modelling and experiments demonstrate that for each transport substance, the gating threshold-the pressure needed to open the pores-can be rationally tuned over a wide pressure range. This enables us to realize in one system differential response profiles for a variety of liquids and gases, even letting liquids flow through the pore while preventing gas from escaping. These capabilities allow us to dynamically modulate gas-liquid sorting in a microfluidic flow and to separate a three-phase air-water-oil mixture, with the liquid lining ensuring sustained antifouling behaviour. Because the liquid gating strategy enables efficient long-term operation and can be applied to a variety of pore structures and membrane materials, and to micro- as well as macroscale fluid systems, we expect it to prove useful in a wide range of applications.

Hou, X↗

Separating Dynamic and Thermodynamic Impacts of Climate Change on Daytime Convective Development over Land

Climate change affects the dynamics and thermodynamics of moist convection. Changes in the dynamics concern, for instance, an increase of convection strength due to increases of convective available potential energy (CAPE). Thermodynamics involve increases in water vapor that the warmer atmosphere can hold and convection can work with. Small-scale simulations are conducted to separate these two components for daytime development of unorganized convection over land. The simulations apply a novel modeling technique referred to as the piggybacking (or master–slave) approach and consider the global climate model (GCM)-predicted change of atmospheric temperature and moisture profiles in the Amazon region at the end of the century under a business-as-usual scenario. The simulations show that the dynamic impact dominates because changes in cloudiness and rainfall come from cloud dynamics considerations, such as the change in CAPE and convective inhibition (CIN) combined with the impact of environmental relative humidity (RH) on deep convection. The small RH reduction between the current and future climate significantly affects the mean surface rain accumulation as it changes from a small reduction to a small increase when the RH decrease is eliminated. The thermodynamic impact on cloudiness and precipitation is generally small, with the extreme rainfall intensifying much less than expected from an atmospheric moisture increase. These results are discussed in the context of previous studies concerning climate change–induced modifications of moist convection. Future research directions applying the piggybacking method are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Halted-pendulum Relaxation: Application to White Dwarf Binary Initial Data

Studying compact-star binaries and their mergers is integral to determining progenitors for observable transients. Today, compact-star mergers are typically studied via state-of-the-art computational fluid dynamics codes. One such numerical technique, smoothed particle hydrodynamics (SPH), is frequently chosen for its excellent mass, energy, and momentum conservation. The natural treatment of vacuum and the ability to represent highly irregular morphologies make SPH an excellent tool for the study of compact-star binaries and mergers. For many scenarios, including binary systems, the outcome of simulations is only as accurate as the initial conditions. For SPH, it is essential to ensure that the particles are distributed regularly, representing the initial density profile but without long-range correlations. Particle noise in the form of high-frequency local motion and low-frequency global dynamics must be damped out. Damping the latter can be as computationally intensive as the actual simulation. We discuss a new and straightforward relaxation method, halted-pendulum relaxation (HPR), to remove global oscillation modes of SPH particle configurations. In combination with effective external potentials representing gravitational and orbital forces, we show that HPR has an excellent performance in efficiently relaxing SPH particles to the desired density distribution and removing global oscillation modes. We compare the method to frequently used relaxation approaches and test it on a white dwarf binary model at its Roche-lobe overflow limit. We highlight the importance of our method in achieving accurate initial conditions and its effect on achieving circular orbits and realistic accretion rates when compared with other general relaxation methods.

79 ASTRONOMY AND ASTROPHYSICS↗

Radiocarbon constraints imply reduced carbon uptake by soils during the 21st century

Soil is the largest terrestrial carbon reservoir and may influence the sign and magnitude of carbon cycle-climate feedbacks. Many Earth system models (ESMs) estimate a significant soil carbon sink by 2100, yet the underlying carbon dynamics determining this response have not been systematically tested against observations. For this research, we used 14 C data from 157 globally distributed soil profiles sampled to 1-meter depth to show that ESMs underestimated the mean age of soil carbon by a factor of more than six (430 ± 50 years versus 3100 ± 1800 years). Consequently, ESMs overestimated the arbon sequestration potential of soils by a factor of nearly two (40 ± 27%). Finally, these inconsistencies suggest that ESMs must better represent carbon stabilization processes and the turnover time of slow and passive reservoirs when simulating future atmospheric carbon dioxide dynamics.

54 ENVIRONMENTAL SCIENCES↗