SEARCH · Search NASA

Results for “Frozen”

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.

SourceResearch connectionFrom yearThrough year

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

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2024

The liquid-conduction, vapor-flow heat pipe model in Sockeye

A single-phase heat pipe flow model implemented in the heat pipe application Sockeye is described. This model solves one-dimensional, compressible flow equations for the vapor phase in the center of a heat pipe, which are coupled to the two-dimensional heat conduction equation for the wick, liquid, and cladding, as well as to an ordinary differential equation tracking the working fluid inventory in the evaporator section of the heat pipe. This model is demonstrated with several test problems, including comparisons to analytic solutions for the vapor flow fields, analytic curves for sonic and capillary limitations of heat pipe operation, and some experimental data. The numerical solution gives excellent agreement for verification problems and good agreement with experimental results. Also, demonstrations show that the model is very robust, allowing for full simulations of heat pipe transients, including frozen startup, sonic-limited (supersonic) flow, and heat pipe shutdown.

42 ENGINEERING↗

DOE OSTI2025

Phenomena Identification and Ranking Table (PIRT) for heat pipes

Heat pipes are advanced passive thermal management devices that utilize phase change and capillary action to achieve efficient heat transfer. However, due to the complexity of the phenomena coupled in heat pipes, including capillary, phase change, turbulence, and compressibility effects, there are high uncertainties in the predictability of their operational regimes and performance. This PIRT exercise, conducted as a collaborative effort involving the Department of Energy (DOE) Microreactor Program (MRP), the Nuclear Regulatory Commission (NRC), and university partners systematically identifies, reviews, and prioritizes critical phenomena affecting the operation of heat pipes based on their importance and knowledge levels. Additional analyses and discussion are provided for phenomena with high importance and low knowledge, such as wick de-wetting, critical heat flux, contact angles, and pressure dynamics. The discussions included the recognizing challenges and proposing future research directions for both modeling and simulation and experimental efforts. Additionally, the report addresses phenomena with medium importance and low knowledge that could impact heat pipe operation during non-normal or transient operation, including frozen startup, laminar to turbulent transition, geyser boiling, wick priming, underfilling conditions, surface roughness of the wick, NCGs trapped in the wick, and the timescales of startup and shutdown. In conclusion, this comprehensive evaluation serves as a valuable resource for guiding future research and development efforts, supporting the successful integration of heat pipes into critical applications such as nuclear reactors, and contributing to the advancement of heat pipe technologies in safety-critical industries.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

DOE OSTI2026

Modeling heat pipe startup and noncondensable gases in Sockeye

For this work, a one-dimensional gas mixture flow model was developed and implemented in the heat pipe code Sockeye to model the effects of noncondensable gases. Additionally, a startup model based on the dusty gas model was implemented to model the transition from rarefied gas dynamics to continuum flow, which occurs during the frozen startup of high-temperature heat pipes. Multiple startup and noncondensable gas models were tested against experimental data for sodium heat pipes, showing excellent agreement. Additionally, the newly developed gas mixture model for modeling noncondensable gas is further tested with a theoretical case study with arbitrary heating configurations. Finally, several recommendations and conclusions are made from the studies in this work to guide future heat pipe modeling efforts.

97 - MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Depth-specific distribution of bacterial MAGs in permafrost active layer in Ny Ålesund, Svalbard (79°N)

Arctic soil microbial communities may shift with increasing temperatures and water availability from climate change. We examined temperature and volumetric liquid water content (VWC) in the upper 80 cm of permafrost-affected soil over 2 years (2018–2019) at the Bayelva monitoring station, Ny Ålesund, Svalbard. We show VWC increases with depth, whereas in situ temperature is more stable vertically, ranging from -5°C to 5°C seasonally. Prokaryotic metagenome-assembled genomes (MAGs) were obtained at 2–4 cm vertical resolution collected while frozen in April 2018 and at 10 cm vertical resolution collected while thawed in September 2019. The most abundant MAGs were Acidobacteriota, Actinomycetota, and Chloroflexota. Actinomycetota and Chloroflexota increase with depth, while Acidobacteriota classes Thermoanaerobaculia Gp7-AA8, Blastocatellia UBA7656, and Vicinamibacteria Vicinamibacterales are found above 6 cm, below 6 cm, and below 20cm, respectively. All MAGs have diverse carbon-degrading genes, and Actinomycetota and Chloroflexota have autotrophic genes. Further, genes encoding β -glucosidase, N-acetyl-β-D-glucosaminidase, and xylosidase increase with depth, indicating a greater potential for organic matter degradation with higher VWC. Acidobacteriota dominate the top 6 cm with their classes segregating by depth, whereas Actinomycetota and Chloroflexota dominate below ~6 cm. This suggests that Acidobacteriota classes adapt to lower VWC at the surface, while Actinomycetota and Chloroflexota persist below 6cm with higher VWC. This indicates that VWC may be as important as temperature in microbial climate change responses in Arctic mineral soils. Here we describe MAG-based Seqcode type species in the Acidobacteriota, Onstottus arcticum, Onstottus frigus, and Gilichinskyi gelida and in the Actinobacteriota, Mayfieldus profundus.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Activity-induced migration of viscous droplets on a solid substrate

Active matter exploits motion to induce changes in shape and conformation via external input. Here, in this paper, we establish theoretically that viscous liquid droplets containing magnetic nanoparticles with frozen-in magnetic moments, sitting on a solid substrate and surrounded by an ambient gas phase, can deform and migrate under the influence of a magnetic torque. The effect arises because the collective rotation of the magnetic nanoparticles at the liquid–gas interface tilts the droplet away from a symmetric configuration, breaks the reflection symmetry with respect to the centre axis, and leads to a left–right asymmetry of the contact angles. A sufficiently strong magnetic torque leads the contact angles to overcome hysteresis effects leading the droplet to migrate. We develop a general framework to explain how symmetry-breaking affects droplet migration. Thus previous results of droplet spreading and migration can be recovered as special cases. Such droplets can be employed as agents in active surfaces and can move against gravity, chemical and thermal gradients, providing a mechanism that could be utilized by both industry and medicine.

Aggarwal, A. [Northwestern Univ., Evanston, IL (Un↗

DOE OSTI2024

Site Disorder Drives Cyanide Dynamics and Fast Ion Transport in Li 6 PS 5 CN

Halide argyrodite solid-state electrolytes of the general formula Li 6 PS 5 X exhibit complex static and dynamic disorder that plays a crucial role in ion transport processes. Here, we unravel the rich interplay between site disorder and dynamics in the plastic crystal argyrodite Li 6 PS 5 CN and the impact on ion diffusion processes through a suite of experimental and computational methodologies, including temperature-dependent synchrotron powder X-ray diffraction, AC electrochemical impedance spectroscopy, 7 Li solid-state NMR, and machine learning-assisted molecular dynamics simulations. Sulfide and (pseudo)halide site disorder between the two anion sublattices unilaterally improves long-range lithium diffusion irrespective of the (pseudo)halide identity, which demonstrates the importance of site disorder in dictating bulk ionic conductivity in the argyrodite family. Furthermore, we find that anion site disorder modulates the presence and time scales of cyanide rotational dynamics. Ordered configurations of anions enable fast, quasi-free rotations of cyanides that occur on time scales of 10 11 Hz at T = 300 K. In contrast, we find that cyanide dynamics are slow or frozen in Li 6 PS 5 CN when site disorder between the cyanide and sulfide sublattices is present at T = 300 K. We rationalize the observed differences in cyanide dynamics in the context of elastic dipole interactions between neighboring cyanide anions and local strain induced by the configurations of site disorder that may impact the energetic landscape for cyanide rotational dynamics. Through this study, we find that anion disorder plays a decisive role in dictating the extent and time scales of both lithium ion and cyanide dynamics in Li 6 PS 5 CN.

36 MATERIALS SCIENCE↗

DOE OSTI2024

High-Throughput Screening for Boride Superconductors

A high-throughput screening using density functional calculations is performed to search for stable boride superconductors from the existing materials database. The workflow employs the fast frozen-phonon method as the descriptor to evaluate the superconducting properties quickly. Twenty-three stable candidates were identified during the screening. The superconductivity was obtained earlier experimentally or computationally for almost all found binary compounds. Previous studies on ternary borides are very limited. Here our extensive search among ternary systems confirmed superconductivity in known systems and found several new compounds. Among these discovered superconducting ternary borides, TaMo 2 B 2 shows the highest superconducting temperature of ∼12 K. Most predicted compounds were synthesized previously; therefore, our predictions can be examined experimentally. Our work also demonstrates that the boride systems can have diverse structural motifs that lead to superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Magnetic and EPR Spectroscopic Studies of Thiolate Bridged Divalent Ni, Pd, and Pt Ions Capped with VO(N 2 S 2 ) Metalloligands

Reactions of the metallodithiolate complex VO(bme-dach) (hereafter abbreviated as V, where bme-dach = N,N ′-bis(2-mercaptoethyl)-1,4-diazacycloheptane) with [Pd II (CH 3 CN) 4 ](BF 4 ) 2 and [Pt II (CH 3 CN) 4 ](BF 4 ) 2 yield the V–M–V trimetallic compounds [VPdV](BF 4 ) 2 (2) and [VPtV](BF 4 ) 2 (3). Reaction of a similar metalloligand, VO(bme-daco) (hereafter abbreviated as V′ where bme-daco = N,N′-bis(2-mercaptoethyl)-1,5-diazacyclooctane) with [Ni II (CH 3 CN) 6 ](BF 4 ) 2 afforded the related salt [V′NiV′](BF 4 ) 2 (1). X-ray structural analyses revealed that cations in 1, 2, and 3 adopt a stairstep C 2h structure consisting of two terminal VO(N 2 S 2 ) moieties bridged via thiolate sulfur to the group 10 metal ions. Weak ferromagnetic superexchange coupling (J = 0.282 cm –1 for 1, 0.954 cm –1 for 2, and 1.372 cm –1 for 3) was observed between the two S = 1/2 V IV centers separated by distances in the range of 5.9–6.3 Å, with J values varying following the order Ni < Pd < Pt. Frozen-solution EPR spectra measured on the more soluble [BArF 24 ] − (BArF 24 – = tetrakis((3,5-trifluoromethyl)phenyl)borate) analogues revealed that the [VPtV] 2+ cation exhibits a 15-line hyperfine splitting of 225 MHz at g = 4 in parallel mode, confirming exchange coupling between the two 51 V, I = 7/2 nuclei. Density-functional theory (DFT) calculations indicate an S = 1 ground state for 1–3. These results demonstrate that the choice of paramagnetic metallodithiolate ligand and diamagnetic bridge in such trimetallic species influences the sign and magnitude of magnetic interactions.

anions↗

DOE OSTI2025

Gold-Standard Chemical Database 137 (GSCDB137): A Diverse Set of Accurate Energy Differences for Assessing and Developing Density Functionals

We present GSCDB137, a rigorously curated benchmark library of 137 data sets (8377 entries) covering main-group and transition-metal reaction energies and barrier heights, (intra- and intermolecular) noncovalent interactions, dipole moments, polarizabilities, electric-field response energies, and vibrational frequencies. Legacy data from GMTKN55 and MGCDB84 have been updated to today's best reference values; redundant or low-quality points were removed, and many new, property-focused sets were added. Testing 29 popular density functional approximations (DFAs) confirms the expected Jacob's-ladder hierarchy overall but also reveals notable exceptions: functional performance for frequencies and electric-field properties correlates poorly with that for other ground-state energetics. ωB97M-V and ωB97X-V are the most balanced hybrid meta-GGA and hybrid GGA, respectively; B97M-V and revPBE-D4 lead the meta-GGA and GGA classes. Double hybrids lower mean errors by about 30% versus their hybrid analogues but demand careful frozen-core, basis set, and spin contamination treatment. GSCDB137 offers a comprehensive, openly documented platform for rigorous validation of DFA and universal machine learning potentials, and training of the next generation of exchange-correlation functionals.

Liang, Jiashu [University of California, Berkeley,↗

DOE OSTI2026

Reduced-Cost Four-Component Relativistic Double Ionization Potential Equation-of-Motion Coupled-Cluster Approaches with 4-Hole–2-Particle Excitations and Three-Body Clusters

The double ionization potential (DIP) equation-ofmotion (EOM) coupled-cluster (CC) method with 4-hole−2- particle (4h-2p) excitations on top of the CC with singles, doubles, and triples calculation, abbreviated as DIP-EOMCCSDT(4h-2p), along with its perturbative DIP-EOMCCSD(T)(a)(4h-2p) approximation, are extended to a relativistic four-component (4c) framework. In addition, we introduce and test a new computationally practical DIP-EOMCC approach, which we call DIPEOMCCSD( T)(ã)(4h-2p), that approximates the treatment of 4h- 2p correlations within the DIP-EOMCCSD(T)(a)(4h-2p) method and reduces the $\mathcal{N}$ 8 scaling characterizing DIP-EOMCCSDT(4h- 2p) and DIP-EOMCCSD(T)(a)(4h-2p) to $\mathcal{N}$ 7 with the system size $\mathcal{N}$. Further improvements in computational efficiency are obtained using the frozen natural spinor (FNS) approximation to reduce the numbers of unoccupied spinors entering the correlated steps of the DIP-EOMCC calculations according to a well-defined occupation-number-based threshold. The resulting 4c-FNS-DIPEOMCC approaches are used to compute DIPs for the series of inert gas atoms from argon to radon as well as the vertical DIPs in Cl 2 , Br 2 , HBr, and HI, which have been experimentally examined in the past. We demonstrate that, when using complete basis set extrapolations and FNS truncation threshold of 10 −4.5 , the 4c-FNS-DIP-EOMCCSD(T)(ã)(4h-2p) calculations are capable of predicting DIPs in agreement with experimental data, improving upon their nonrelativistic and spin-free scalar-relativistic counterparts, particularly when examining DIPs characterized by stronger spin−orbit coupling effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Defect Stability in CdTe Based on Formation Energies and Migration Barriers

Native point defects are thought to play a key role in CdTe, either as compensation centers in intentionally doped material, as a source of conductivity in nominally undoped material, or as electron–hole recombination centers. Here, the discussion of their concentration and impact has often centered only on formation energies and transition levels. Using hybrid density functional calculations, including the effects of spin–orbit coupling (SOC), we discuss the stability of native point defects in CdTe based on their formation energies and migration barriers. We show that although Cd interstitials are the lowest energy donor defects, they are unstable at room temperature due to a low migration barrier. They are important for maintaining charge neutrality during growth or annealing at high temperatures, but once the material is brought to room temperature, they are not frozen in as often assumed and are expected to anneal out, leaving the other more stable defects to determine the conductivity. Taking this into account in the solution of the charge neutrality equation, we are able to predict the conductivity type and carrier concentrations that are in good agreement with experimental observations.

14 SOLAR ENERGY↗

DOE OSTI2024

Electrified Operando -Freezing of Electrocatalytic CO 2 Reduction Cells for Cryogenic Electron Microscopy

The ability to freeze and stabilize reaction intermediates in their metastable states and obtain their structural and chemical information with high spatial resolution would be very powerful to unravel the fundamentals in many important materials technologies such as catalysis and batteries. Here, we develop an electrified operando-freezing methodology for the first time to preserve these metastable states under electrochemical reaction conditions for cryogenic electron microscopy (cryo-EM) imaging and spectroscopy. Using Cu catalysts for CO 2 reduction as a model system, we observe restructuring of the Cu catalyst in a CO 2 atmosphere while the same catalyst remains intact in an air atmosphere at the nanometer scale. Furthermore, we discover the existence of single valance Cu (1+) state and C-O bonding at the electrified liquid-solid interface of the operando-frozen samples, which are key reaction intermediates that traditional ex situ measurements fail to detect. Finally, this work highlights our novel technique to study the local structure and chemistry of electrified liquid-solid interfaces, which has broad impact for many electrochemical reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Kinetic Understanding of Field-Induced Phase Transition from Tetragonal to Ferroelectric Orthorhombic Phase in Ferroelectric CeO 2 –HfO 2 –ZrO 2 Films

The ferroelectric properties and structural phase transition behaviors of fluorite-type CeO 2 −HfO 2 −ZrO 2 films were investigated. The epitaxial films on indium tin oxide (ITO) (111)/yttria-stabilized zirconia (YSZ) (111) substrates were grown through pulsed laser deposition at room temperature and subsequently heat-treated at 1000 °C under a N 2 gas flow. The crystalline phases and Curie temperatures of the films were investigated by X-ray diffraction. An increase in the Ce or Zr content in the films led to a higher crystallographic symmetry, such as orthorhombic or tetragonal. In addition, electrical characterization revealed that the orthorhombic films and some of the tetragonal films displayed ferroelectricity. This was due to the field-induced phase transition from the tetragonal to ferroelectric orthorhombic phase in the films, where the Curie temperatures were relatively low. The tetragonal metastable phase was kinetically frozen and could not change into the stable orthorhombic phase at such a low temperature. The critical electric field where the field-induced phase transition occurred was below 0.8 MV/cm, which was sufficiently small compared to the coercive field. These results evidence the kinetic driving force that causes a field-induced phase transition from the paraelectric tetragonal phase to the ferroelectric orthorhombic phase in HfO 2 -based ferroelectrics. They also enhance our understanding of the thermodynamic phase stabilities of HfO 2 -based material polymorphs.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Bayesian Optimization of Catalysis with In-Context Learning

Large language models (LLMs) can perform accurate classification with zero or few examples through in-context learning (ICL), allowing the model to observe query-relevant examples at inference time and eliminating the need for additional weight updates to generalize beyond its original training data. We extend this capability to regression with uncertainty estimation using frozen LLMs (e.g., GPT-4o, Gemini), enabling Bayesian optimization (BO) in natural language without explicit model training or feature engineering. We apply this to materials discovery by representing materials as synthesis and testing procedures for use in natural language prompts. This Bayesian, design-first approach prioritizes optimization toward target material properties before detailed characterization, in contrast to conventional experimental workflows that often emphasize characterization of suboptimal materials. On benchmarks like aqueous solubility and oxidative coupling of methane (OCM), BO-ICL matches or outperforms Gaussian processes. In live experiments on the reverse water–gas shift (RWGS) reaction, BO-ICL identifies multimetallic catalysts that approach equilibrium CO yield within 6 and 10 iterations from a pool of 3,700 and 360,000 candidates, respectively. Our method redefines materials representation and accelerates discovery, with broad applications across catalysis, materials science, and AI.

Calibration↗

DOE OSTI2024

Catalytic Deconstruction of Ethylene Vinyl Acetate Copolymer and Polyethylene Mixtures via Hydroconversion: Challenges and Solutions

We explore hydrogenolysis over ruthenium supported on zirconia (Ru/ZrO 2 ) and hydrocracking over platinum (Pt) supported on zeolites as an effective end-of-life strategy for ethylene vinyl acetate (EVA)–a widely used performance heat sealant in hard-to-recycle multilayer packaging. For Ru/ZrO 2 hydrogenolysis, EVA reacts slower than low-density polyethylene (LDPE) and the catalyst deactivates due to carbonaceous deposits originating from polyenes generated in situ during EVA thermal degradation. High H 2 pressures and temperatures can overcome catalyst deactivation; however, CH 4 yields are excessive due to cascade hydrogenolysis stemming from strong C=C/metal interactions. Polyene hydrogenation allows chains anchored by C=C to desorb from Ru, shifting product selectivity from CH 4 to higher-value liquids. Hydrogenolysis of mixed EVA and linear low-density polyethylene (LLDPE), mirroring typical frozen food packaging formulations, results in comparable catalyst activity and CH 4 yield as the pure EVA resin. For Pt/zeolite hydrocracking, pure EVA and EVA:LLDPE mixtures are deconstructed to propane or light naphtha with minimal CH 4 production. Among catalysts tested, Pt/HY gives the highest liquid productivity (g C5+products /g cat ·h). Furthermore, these findings showcase the recalcitrant nature of EVA and its associated mixtures for Ru/ZrO 2 hydrogenolysis, highlighting that hydrocracking catalysts may be superior for complex packaging waste.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Immersion Freezing in Particle-Based Aerosol-Cloud Microphysics: A Probabilistic Perspective on Singular and Time-Dependent Models

Cloud droplets containing immersed ice-nucleating particles (INPs) may freeze at temperatures above the homogeneous freezing threshold temperature in a process referred to as immersion freezing. In modeling studies, immersion freezing is often described using either so-called “singular” or “time-dependent” parameterizations. Here, we compare both approaches and discuss them in the context of probabilistic particle-based (super-droplet) cloud microphysics modeling. First, using a box model, we contrast how both parameterizations respond to idealized ambient cooling rate profiles and quantify the impact of the polydispersity of the immersed surface spectrum on the frozen fraction evolution. Presented simulations highlight that the singular approach, constituting a time-integrated form of a more general time-dependent approach, is only accurate under a limited range of ambient cooling rates. The time-dependent approach is free from this limitation. Second, using a prescribed-flow two-dimensional cloud model, we illustrate the macroscopic differences in the evolution in time of ice particle concentrations in simulations with flow regimes relevant to ambient cloud conditions. The flow-coupled aerosol-budget-resolving simulations highlight the benefits and challenges of modeling cloud condensation nuclei activation and immersion freezing on insoluble ice nuclei with super-particle methods. The challenges stem, on the one hand, from heterogeneous ice nucleation being contingent on the presence of relatively sparse immersed INPs, and on the other hand, from the need to represent a vast population of particles with relatively few so-called super particles (each representing a multiplicity of real particles). We discuss the critical role of the sampling strategy for particle attributes, including the INP size, the freezing temperature (for singular scheme) and the multiplicity.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Exposing and Reducing Biases of Simulating Mixed-Phase Clouds in the Convection-Permitting E3SM Atmosphere Model: Lessons From an Arctic Cold-Air Outbreak

Mixed-phase clouds modulate the water and energy cycles of high-latitude regions, yet their liquid-ice phase partitioning has long been poorly simulated in climate models. Here, simulations of Arctic mixed-phase clouds by the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) are assessed against large-eddy simulations, satellite data, and ground-based observations during the Cold-Air Outbreaks in the Marine Boundary Layer Experiment field campaign. SCREAM simulates nearly completely frozen clouds, which is attributed largely to the unreasonably strong Wegener–Bergeron–Findeisen (WBF) process that converts liquid to ice excessively and partly to the early over-abundant ice production at cold temperatures from a temperature-deterministic deposition ice nucleation scheme. Assuming no subgrid variation for the WBF process in the original formulation particularly conflicts with the instantaneous saturation adjustment assumption in the condensation scheme that assumes subgrid variability, leading to exaggerated WBF process rates. A proposed simple physically-based improvement on the treatment of subgrid cloud overlap substantially increases supercooled liquid water content and notably improves cloud-top phase partitioning, aligning better with observations. Improvement of supercooled liquid water content also converges with increasing horizontal resolution. The deposition ice nucleation scheme is found responsible for a falsely-produced ice cloud aloft that is not observed, biasing the simulated cloud radiative effects and top-of-atmosphere radiative fluxes. This study identifies key deficiencies in cloud parameterizations that continue to challenge convection-permitting models.

Geosciences↗

DOE OSTI2024

Whole-soil warming leads to substantial soil carbon emission in an alpine grassland

The sensitivity of soil organic carbon (SOC) decomposition in seasonally frozen soils, such as alpine ecosystems, to climate warming is a major uncertainty in global carbon cycling. Here we measure soil CO 2 emission during four years (2018–2021) from the whole-soil warming experiment (4 °C for the top 1 m) in an alpine grassland ecosystem. We find that whole-soil warming stimulates total and SOC-derived CO 2 efflux by 26% and 37%, respectively, but has a minor effect on root-derived CO 2 efflux. Moreover, experimental warming only promotes total soil CO 2 efflux by 7-8% on average in the meta-analysis across all grasslands or alpine grasslands globally (none of these experiments were whole-soil warming). We show that whole-soil warming has a much stronger effect on soil carbon emission in the alpine grassland ecosystem than what was reported in previous warming experiments, most of which only heat surface soils.

54 ENVIRONMENTAL SCIENCES↗